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We present a detailed report of a measurement of the neutron /3-asymmetry parameter Ao, the 
parity-violating angular correlation between the neutron spin and the decay electron momentum, 
performed with polarized ultracold neutrons (UCN). UCN were extracted from a pulsed spallation 
solid deuterium source and polarized via transport through a 7-T magnetic field. The polarized UCN 
were then transported through an adiabatic-fast-passage spin-flipper field region, prior to storage in 
a cylindrical decay volume situated within a 1-T 2 x 2tt solenoidal spectrometer. The asymmetry 
was extracted from measurements of the decay electrons in multiwire proportional chamber and 
plastic scintillator detector packages located on both ends of the spectrometer. From an analysis 
of data acquired during runs in 2008 and 2009, we report A = -0.11966 ± 0.00089 tHoilo, from 
which we extract a value for the ratio of the weak axial-vector and vector coupling constants of the 
nucleon, A = gA/gv ~ —1.27590 ± 0.00239 !o'oo377- Complete details of the analysis are presented. 

PACS numbers: 12.15.Ff, 12.15.Hh, 13.30.Ce, 14.20.Dh, 23.40.Bw 



I. INTRODUCTION 



Precise measurements of neutron /3-decay observables 
determine fundamental parameters of the weak interac- 
tion and contribute to tests of the Standard Model [1- 
6]. Because the momentum transfer in neutron /3-decay 
(n — > p + er + 77 e + 781.5 keV) is small compared to the 
W~ mass, the decay can be modeled as a four-fermion 
contact interaction with an amplitude under the Stan- 
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where q is the four-momentum transfer; M is the nucleon 
mass; and g v {q 2 ), 3wm(<? 2 ), 3s(<? 2 ), 3a(<? 2 ), 3t(<? 2 ), and 



dard Model given by 

M= GM* w>)Liii (1) 

where Gf is the Fermi weak coupling constant, V u d is 
the weak-quark- mixing CKM matrix element, and = 
Tie7 J u(l — 75)1^77 is the leptonic weak vector and axial- 
vector current. In its most general form, the hadronic 
weak vector and axial-vector (V — A) current includes 
six form factors [7, 8] , 



u n , (2) 

I 

gp(q 2 ) are the vector, weak magnetism, induced scalar, 
axial vector, induced tensor, and induced pseudoscalar 
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form factors, respectively. In the limit of q 2 — > 0, the 
hadronic weak current is dominated by the weak vector 
and axial vector coupling constants of the nucleon, de- 
fined to be the values of the vector and axial vector form 
factors at q 2 = 0, gA = <?a(<Z 2 = 0) and gv = gv(l 2 = 0). 
Under the Conserved Vector Current (CVC) hypothesis 
of the Standard Model and the assumption of isospin 
symmetry, the vector coupling constant is gv = 1 (in- 
dependent of the nuclear medium). Isospin-symmetry- 
breaking effects to the value of gv in neutron /3-decay 
have been calculated in chiral perturbation theory, with 
the correction to gv found to be at a negligible —4 x 10 -5 
level [9]. Also per the CVC hypothesis, the weak mag- 
netism coupling constant, gwM = gwAi(q 2 = 0), which 
appears at recoil order in the vector current, is related 
to the proton and neutron anomalous magnetic moments 
by gwM — k p — K n . 

In contrast to the vector current, the axial- vector cur- 
rent is renormalized by the strong interaction such that 
the value of gA must be determined experimentally and 
also by lattice Quantum Chromodynamics (QCD) calcu- 
lations. Any contribution from the induced pseudoscalar 
coupling constant, gp, to neutron /3-decay observables 
is expected to be negligibly small, with the contribution 
of gp to the energy spectrum calculated to be of order 
m 2 e /ME e ~ 10~ 4 [10]. 

The two remaining terms, the induced scalar, gs(q 2 ), 
in the vector current, and the induced tensor, gr(q 2 ), 
in the axial- vector current, are termed second-class cur- 
rents due to their transformation properties under G- 
parity. Under the requirement of G-parity symmetry, 
both gs(q 2 ) — St(<Z 2 ) = 0. However, G-parity symmetry 
is violated within the Standard Model due to differences 
in the u and d quarks' charges and masses (i.e., isospin 
symmetry breaking effects). An estimate for <?t includ- 
ing SU(3) breaking effects suggested \gr\ in neutron (3- 
decay to be < 0.03 [11], and an evaluation of gx/gA using 
QCD sum rules found g T /gA = -0.0152(53) [12]. Fi- 
nally, recent lattice QCD studies of SU(3) breaking in 
semi-leptonic decays find small, O(0. 1), values for both 
gs{g 2 ) and gr(q 2 ) in neutron /3-decay, but the results 
are statistically limited and consistent with zero at 1-2 
standard deviations [13]. However, despite these hints 
for non-zero values of these second-class currents, their 
contributions to neutron /3-dccay observables are again 
expected to be negligibly small, as they also appear at 
order m 2 /ME e ~ 10~ 4 in the energy spectrum [10]. 

Therefore, under the assumption that any such contri- 
butions from gp, g$, and <?t are negligibly small relative 
to the current level of experimental precision, it is clear 
that a description of neutron /3-decay under the CVC hy- 
pothesis of the Standard Model requires the specification 
of only two parameters, V u d and gA, given the high pre- 
cision results for Gf achieved in muon decay [14]. Both 
V u d and gA can be accessed via measurements of two dif- 
ferent types of neutron /3-decay observables: the lifetime, 
and angular correlation coefficients in polarized and un- 
polarized /3-decay. The first of these, the lifetime, as 



calculated from the amplitude and integration over the 
allowed phase space, is of the form [15] 



1 



2tt 3 



■K 2 d(l + 3A 2 )/(l + RC) 



(3) 



where m e is the electron mass and the parameter A is 
defined to be the ratio of the axial vector and vector cou- 
pling constants, A = gA/gv- The numerical value for 
the phase space factor of / = 1.6887 [15] includes cor- 
rections for the Fermi function, the finite nucleon mass, 
the finite nucleon radius, and the effect of recoil on the 
Fermi function. The factor (1 + RC) denotes the total 
effect of all electroweak radiative corrections, including 
the 0(a) outer (long-distance loop and bremsstrahlung 
effects) and inner (short distance, including axial- vector- 
current, loop effects) radiative corrections; an 0(a 2 ) cor- 
rection resulting from factorization of the Fermi function; 
and 0(a 2 ) leading- log and next-to-leading- log correc- 
tions (for lepton and quark loop insertions in the photon 
propagator) [15]. The total electroweak radiative correc- 
tion has been calculated to be (1 + RC) = 1.0390±0.0004 
[16] where the ±0.0004 uncertainty was reduced by a 
factor of two (from its previous value of ±0.0008 [15]) 
after the development of a new method for calculat- 
ing hadronic effects in the matching of long- and short- 
distance contributions to axial- vector current loop effects 
(primarily from the box diagram). 

The second type of observable, angular correlation co- 
efficients, parametrize the angular correlations between 
the momenta of the decay products and the spin of the 
initial-state neutron. In general, the directional distribu- 
tion of the electron and antineutrino momenta and the 
electron energy in polarized /3-decay is of the form [17] 
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where E e (E v ) and p e (p v ) denote, respectively, the elec- 
tron's (antineutrino's) total energy and momentum; Eq 
(= 781.5 keV + m e ) is the electron endpoint energy; and 
(a n ) is the neutron polarization. The angular correlation 
coefficients a (e-j7 e -asymmetry), A (/3-asymmctry), and 
B (77 e -asymmetry) are, to lowest order, functions only of 
A where, under a A < sign convention, 



ao = 



1-A 2 
1 + 3A 2 ' 



A = -2 



A(A±1) 
1 + 3A 2 1 



B = 2 



A(A-l) 
1 + 3A 2 ' 
(5) 



The contributions of terms in Eq. (4) proportional 
to the Fierz interference term b and the time-reversal- 
odd triplc-correlation-cocfficient D are at recoil order for 
Standard Model interactions [10, 18-20], and are negligi- 
ble at the current level of experimental precision. Note 



3 



that to our knowledge, there arc no published direct mea- 
surements of b in neutron /3-decay. 

As already noted, at lowest order ao, Aq, and Bq are 
functions only of A. However, recoil-order corrections, 
including the effects of weak magnetism and gy-gA in- 
terference, introduce energy- dependent corrections to the 
asymmetry, and are of 0(1%) for a and A [10, 21, 22] and 
0(0.1%) for B [23]. For A, the recoil-order corrections 
are of the explicit functional form [10, 21, 22] 



A = An 



Rx 



A 2 R + A 3 Rx, 



(6) 



where R = E /M, e = (m e /M) 2 , x = E e /E , and the A, 
(z = 1, 2, 3) coefficients arc functions only of A and gwu 
(assuming gg = g? = 0, and negligible contributions 
from g P ). Under these assumptions, A\ = —0.3054, A 2 = 
0.7454, and A 3 = -3.0395. Note that the q 2 dependence 
of the form factors does not appear until next-to-leading 
recoil order [22]. 

In addition to the above recoil-order corrections, there 
is a small energy-dependent radiative correction (for vir- 
tual and bremsstrahlung processes) to polarized asym- 
metries, resulting in a O(0.1%) correction to A [24, 25]. 
After application of these recoil-order and radiative cor- 
rections, a value for A can be extracted from ao, Aq, and 
Bq. Note, however, that for a given (relative) statistical 
precision, the sensitivity of A to A is slightly higher than 
that of ao, and a factor of ~ 8 greater than that of Bo, 
where at leading order the relative uncertainties compare 
as 
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Thus, measurements of the angular correlation coeffi- 
cients determine a value for A (or gA, assuming gy = 1 
under the CVC hypothesis), a fundamental parameter in 
the nucleon weak current. A precise value for gA is also 
important in many other contexts. In hadronic physics 
studies of the spin structure of the nucleon [26, 27], the 
Bjorken sum rule relates the difference in the first mo- 
ments of the proton and neutron spin-dependent g\ struc- 
ture functions (i.e., isovector channel), as probed in po- 
larized deep inelastic electron scattering, to gA- In QCD, 
the assumption of a partially conserved axial- vector cur- 
rent (PC AC), valid in the limit of a massless pion (iden- 
tified as the Goldstone boson of the spontaneously bro- 
ken chiral symmetry), leads to the Goldberger-Treiman 
relation [7], relating the value of gA to the pion decay 
constant /„•, the weak pion- nucleon- nucleon coupling con- 
stant g n NNi and the nucleon mass. The value of gA is also 
important in astrophysical processes, including calcula- 
tions of solar fusion cross sections and rates, in particular, 
of the pp fusion reaction, impacting the solar neutrino 
flux for this process [28] . High-precision experimental re- 
sults for gA also serve as an important benchmark for 
theoretical calculations of gA, both in fundamental lat- 
tice QCD calculations [29] and in rclativistic constituent 
quark model calculations [30]. A precise value for gA is 



also important as a phcnomcnological input parameter 
(together with other low energy constants, such as the 
pion decay constant f w , the nucleon mass, etc.) to ef- 
fective field theory calculations involving the axial vector 
current [31]. 

Although not a fundamental weak interaction param- 
eter by itself, a precise value for the lifetime is impor- 
tant for Big Bang Nucleosynthesis calculations, impact- 
ing the neutron-to-proton ratio and hence the primordial 
4 He abundance at the time of freeze-out, when the weak 
reaction rate became less than the Hubble expansion rate 
[32]. The value of the lifetime is also important for the 
interpretation of data from neutrino oscillation experi- 
ments employing antincutrinos from reactors [33], which 
typically search for the reaction 77 e +p — > e + +n in detec- 
tors. The cross section for this reaction is inversely pro- 
portional to the neutron lifetime; therefore, an accurate 
and precise experimental value for the lifetime is needed 
for an interpretation of measured detector antincutrino 
reaction rates in terms of the underlying neutrino oscil- 
lation physics. 

Measurements of the lifetime and a value for A from 
measurements of angular correlation coefficients permit 
the extraction of a value for V u d solely from neutron /3- 
decay observables according to Eq. (3). Although a value 
for V u d from neutron /3-decay [34] is not yet competitive 
with the definitive value deduced from measurements of 
ft values in superallowed + — > + nuclear /3-dccay [35], 
the appeal of such an extraction is that it does not require 
corrections for isospin-symmetry-brcaking and nuclear- 
structure effects. Ultimately, when the precision on a 
neutron-based value for V u d approaches the precision of 
the + — > + result, the two values must agree in the ab- 
sence of new physics. However, given that the neutron- 
based value for V u d is not yet competitive, one can treat 
the + — > + value for V u d as a fixed input parameter, 
and instead perform a robust test of the consistency of 
the various measured neutron /3-dccay observables un- 
der the Standard Model. In particular, results for gA 
extracted from correlation coefficient measurements can 
then be directly compared with results from measure- 
ments of the neutron lifetime r„. 

Finally, measurements of the angular correlation coef- 
ficients themselves are sensitive to beyond-the-Standard- 
Model physics, such as scalar and tensor interactions 
[18, 36, 37]. With the projected improvements to the 
experimental precision in future years, neutron /3-decay 
measurements will be sensitive to any such sources of new 
physics at energy scales rivaling those probed directly at 
the Large Hadron Collider [18]. 

The remainder of this article is organized as follows. In 
Section II, we summarize the current status of measure- 
ments of Aq. We then outline the experimental motiva- 
tion for a measurement of Aq with ultracold neutrons in 
Section III, and then present a detailed description of the 
UCNA ("Ultracold Neutron Asymmetry") Experiment 
[38, 39] at the Los Alamos National Laboratory. Our 
measurement procedures and experimental geometrical 
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TABLE I. Summary of published measurements of the neutron /3-asymmetry parameter Ao. The error on the average has been 
increased by a factor of \fx 2 / v = 2.47. 



Experiment 


Years Published 


Type 


Polarization 


Ao Result 


Notes 


PERKEO [40] 


1986 


cold neutron beam 


0.974 ± 0.005 


-0.1146 ±0.0019 


a 


PNPI [41, 42] 


1991, 1997 


cold neutron beam 


0.770 ± 0.007 


-0.1135 ±0.0014 


b 


ILL-TPC [43, 44] 


1995, 1997 


cold neutron beam 


0.981 ±0.003 


-0.1160 ±0.0015 


c 


PERKEO II [45, 46] 


1997, 2002 


cold neutron beam 


0.989 ± 0.003 


-0.1189 ±0.0007 


d 


UCNA [38, 39], this work 


2009, 2010 


stored ultracold neutrons 


1 +0 

i - u -0.O052 


-0.11966 ± 0.00089lo;ooi4o 


c 


PERKEO II [47] 


2012 


cold neutron beam 


0.997 ±0.001 


-0.11996 ±0.00058 


r 




Current Average 


Value: A = -0.11846 ± 0.00104 (x^/v = 


24.35/4) 





a Included a ~ 10% correction to the asymmetry for magnetic mirror effects. 

b The result reported in [42] superseded that reported in [41] of Ao = —0.1116(14), on the basis of a revised value for the polarization. 
c The final result reported in [44] was identical to a first result reported in [43]. 

d The final result of A = -0.1189(7) was the combined result of -0.1189(12) reported in [45] and -0.1189(8) reported in [46]. 
° The result of Ao = — 0.1138(46)(21) reported in [38] was from a proof-of-principle measurement and was not included in the result 
reported in [39]. 

f Accounting for correlated systematic errors in [46, 47], the combined PERKEO II result is Aq = —0.11951 ± 0.00050. 



configurations are reported in Section IV. Results from 
our calibration and analysis procedures are discussed in 
Section V. Details of our procedure for the extraction of 
asymmetries are presented in Section VI, and the correc- 
tions to measured asymmetries for various systematic ef- 
fects are discussed in Section VII. Systematic uncertain- 
ties are summarized in Section VIII, and our final results 
for Aq are then reported in Section IX. We then conclude 
with a brief summary of the physics impact of our work 
in Section X. The data presented here were obtained dur- 
ing data-taking runs in 2008-2009 and published rapidly 
in 2010 [39]; in this article we provide a more detailed 
account of the experiment and analysis procedures. 
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II. STATUS OF MEASUREMENTS OF A 



FIG. 1. Results from measurements of the /3-asymmetry pa- 
rameter Ao [38-46]. The band (±lo-) indicates the current 
world average value of A Q = -0.11846 ± 0.00104. 



The current status of published results [38-47] for the 
neutron /3-asymmctry parameter Aq is summarized in Ta- 
ble I and shown in Fig. 1. Other than the UCNA Exper- 
iment, all of the experiments have been performed with 
beams of polarized cold neutrons, with reported values 
for the polarization ranging from 0.770 ± 0.007 [41, 42] 
to 0.997 ± 0.001 [47]. Magnetic solenoidal spectrome- 
ters providing 2 x 27r solid angle acceptance for detection 
of the decay electrons were employed in the PERKEO 
[40] and PERKEO II [45-47] experiments at the Institut 
Laue-Langevin (ILL). In contrast, the solid angle was 
defined by the geometric acceptance in an experiment at 
the Petersburg Nuclear Physics Institute (PNPI) [41, 42] 
in which the decay electrons and protons were detected 
in coincidence in detectors surrounding the beam decay 
region, and in an experiment at the ILL [43, 44] which 
utilized a time projection chamber for reconstruction of 
the electron track. 

The current world average value for Aq — — 0.1 1846 ± 



0.00104 includes the most recent PERKEO II result 1 [47] , 
but excludes the UCNA proof-of-principle result [38]. 
Note that the current err or ba r of ±0.00104 includes the 
Particle Data Group's \/x 2 / v scaling [34]. The need for 
this expanded error bar suggests an incomplete assess- 
ment of the systematic errors in one or more of the cold- 
neutron-based experiments. 



1 Note that in computing the world average, we employed the com- 
bined PERKEO II result of -0.11951 ± 0.00050 reported in [47] 
which accounted for correlations of systematic errors in the two 
separately published PERKEO II results [46, 47]. 
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III. UCNA EXPERIMENT 
A. Overview of Experiment 

The UCNA Experiment, installed in Area B of the 
Los Alamos Neutron Science Center (LANSCE) at the 
Los Alamos National Laboratory (LANL), was designed 
to perform the first-ever measurement of the neutron j3- 
asymmetry parameter with ultracold neutrons (UCN), 
and to-date is the only experimental measurement of 
any neutron /3-decay angular correlation coefficient per- 
formed with ultracold neutrons (UCN). UCN are de- 
fined to be neutrons with kinetic energies sufficiently low 
(< 335 neV, corresponding to speeds < 8 m s _1 ) such 
that they undergo total external reflection at any angle 
of incidence from an effective potential barrier (a vol- 
ume average of Fermi potentials Upcrmi) at the surfaces 
of certain materials [48]. Thus, UCN can be stored in 
material-walled vessels, whereas cold neutrons (kinetic 
energies 0.05-25 meV, speeds 100-2200 m s _1 ) must be 
transported along neutron guides at reflection angles less 
than the guide critical angle, resulting in short residency 
times in an apparatus. 

A schematic diagram of the UCNA Experiment is 
shown in Fig. 2, and the basic principle of the experi- 
ment is as follows. Spallation neutrons resulting from the 
interaction of a pulsed (typically 0.2 Hz) 800 MeV pro- 
ton beam with a tungsten target were moderated in cold 
polyethylene to the cold neutron regime, and then down- 
scattered to the UCN regime in a solid deuterium (SD 2 ) 
crystal. The UCN were then transported along a series 
of UCN guides through a 7.0-Tesla solenoidal polarizing 
magnet, where the spin-dependent — /x- B potential (±60 
neV T^ 1 ) served as a spin-state selector for magnetic 
moments pi oriented parallel to the direction of the longi- 
tudinal magnetic field B. The polarized UCN were then 
transported along non-magnetic UCN guides through an 
adiabatic- fast-passage (AFP) spin-flipper 1.0-Tcsla field 
region, used to prepare UCN with spins either parallel or 
anti-parallel to the magnetic field. The UCN were then 
directed to the center of a 12.4-cm diameter, 3-m long 
cylidrical decay storage volume located within the warm 
bore of a 1.0-T solenoidal spectrometer. Emitted /3-decay 
electrons then spiraled (with a maximum Larmor diam- 
eter of 7.76 mm for 782 keV endpoint electrons emitted 
perpendicular to the 1.0-T field) along the magnetic field 
lines towards one of two electron detectors located on 
both ends of the spectrometer. 

In principle, the /3-asymmetry A can be extracted from 
measurement of the W(9) oc (1 + PnA/3 cos 9) angular 
distribution by forming an energy- dependent "measured 
asymmetry", A meas (E e ), of the detectors' (background- 
subtracted) count rates. 

A meas (Ee) = ri \f\l r2 [f\ = P n AP(cose), (8) 
ri(E e ) + r 2 {E e ) 

where ri^)(E e ) denote the energy-dependent count rates 
observed in the two detectors, P n denotes the neutron 



polarization, /3 denotes the electron velocity in units of 
c, and (cos 9} is the average value of cos 9 integrated 
over the detectors' angular acceptance for that partic- 
ular value of E e . Note that for nominal values of P n ~ 1 , 
A ~ -0.12, /3 ~ 0.75, and (cos6>) ~ 1/2, the experimen- 
tal measured asymmetry is of order |A meas | ~ 0.04. 

In practice, the asymmetry is extracted from ratios of 
the two detectors' energy- dependent count rates for the 
two neutron spin states with polarizations oriented par- 




FIG. 2. Schematic diagram (not to scale) of the primary 
components of the UCNA Experiment's /3-asymmetry mea- 
surement. The top panel depicts the layout of the 7.0-Tesla 
polarizing magnet and AFP spin-flipper, the 1.0-Tesla spec- 
trometer, the decay trap, and the electron detectors. The 
bottom panel depicts the UCN source and shielding, the lay- 
out of the UCN transport guides, the location of the UCN 
gate valve, the location of the UCN switcher, the locations 
of all the magnets, and the locations of all the UCN moni- 
tor detectors. All of these components are described in detail 
later in the text. [Note that during our data-taking runs in 
2009, the (bare) Cu UCN guides located between the spin- 
flipper and the entrance to the decay trap (as depicted in the 
top panel) were replaced with diamondlike carbon-coated Cu 
guides; see Section IV B for details.] 
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allel and anti-parallel to the magnetic field via a "super- 
ratio" technique. Here, the super-ratio, i?, is defined in 
terms of the measured energy-dependent detector count 



rates for the two spin states, 



„+(-) 

1(2) 



{E e ), to be 



R 



r^(E e )-r+(E e ) 



(9) 



r+(£ e )-r 2 -(£y 
with the energy- dependent measured asymmetry. 

A 

ing to 



A meas (E e ), then calculated from the super-ratio accord- 



i - Vr 
i + Vr 



P n A(i (cos 9). 



(10) 



The merit of this super-ratio technique is that effects due 
to differences in the two detectors' efficiencies and spin- 
dependent differences in the efficiencies for transport of 
the two UCN spin states into the spectrometer cancel to 
first order. In a binned analysis, energy- dependent de- 
tection efficiencies also largely cancel in the super-ratio, 
and are negligible for the energy bin sizes used in this 
work. 

The motivation for the development of the UCNA Ex- 
periment was several- fold. First, the use of UCN in a 
neutron /3-asymmetry experiment controls key neutron- 
related systematic corrections and uncertainties, includ- 
ing the neutron polarization and neutron-generated back- 
grounds. As discussed in detail later in this article, the 
polarization has been demonstrated to be > 99.48% at 
the 68% C.L., with the precision, at present, limited only 
by statistics. Further, neutron-generated backgrounds 
have been constrained to be negligible, a direct result of 
the relatively small number of neutrons present in the 
apparatus at any time, the small probability for their 
capture and subsequent generation of accompanying irre- 
ducible gamma ray backgrounds, and the fact that nearly 
all of the neutrons present in the apparatus are located 
within the spectrometer's decay volume. In the UCNA 
Experiment, a relatively large fraction, ~ 1/40, of the 
UCN stored in the decay volume contribute to the mea- 
sured decay rate, whereas in cold neutron beam exper- 
iments typically only ~ 1/10 7 of the neutrons passing 
through the apparatus contribute to the decay rate [45]. 
Therefore, control of neutron-generated backgrounds is 
expected to be intrinsically more challenging in cold neu- 
tron beam experiments. 

Second, as described in detail elsewhere [49, 50], the 
electron detector system developed for the UCNA Ex- 
periment, consisting of a low-pressure multiwire propor- 
tional chamber (MWPC) backed by a plastic scintillator, 
provides position sensitivity, suppresses ambient gamma- 
ray backgrounds, and permits the reconstruction of low- 
energy-deposition electron backscattering events. The 
UCNA Experiment is the first neutron /3-asymmetry ex- 
periment to employ a MWPC, providing the experiment 
with two critical advantages. First, the position infor- 
mation permits the definition of a fiducial volume on an 



event- by- event basis. Second, the position information 
also permits for an event- by- event correction for the scin- 
tillator's position-dependent energy response. 

We now provide a more detailed description of the pri- 
mary components of the UCNA Experiment. 



B. UCN Source and Guide Transport System 

A detailed description of the design principles and per- 
formance of the LANL SD2 UCN source is given else- 
where [51-53]; therefore, we provide only a brief descrip- 
tion here. Protons from the 800 MeV LANSCE accel- 
erator were delivered in pulsed mode 2 at a repetition 
rate of 0.2 Hz to a tungsten spallation target, which 
was surrounded by a room-temperature beryllium reflec- 
tor. With the spallation target operated in this pulsed 
mode, prompt beam related backgrounds can be elimi- 
nated with simple timing cuts, with negligible loss of duty 
factor for the /3-decay measurements performed with the 
UCN stored in the electron spectrometer. The spalla- 
tion neutrons were moderated in cold-helium-gas-cooled 
polyethylene (maintained at a temperature of ~ 150 K for 
time-averaged proton beam currents of 5.8 /iA) located 
between the tungsten target and the beryllium reflec- 
tor. The moderated cold neutrons were then downscat- 
tered to the UCN regime in a ~ 2 L cylindrical volume 
of ortho-state SD2 [54, 55] embedded at the bottom of 
a vertically-oriented cylindrical liquid-helium-cooled alu- 
minum cryostat coated with 58 Ni, presenting a nominal 
effective potential of 342 neV to the emerging UCN flux. 
The SD2 was maintained at temperatures < 10 K dur- 
ing the proton beam pulses on the spallation target. A 
butterfly-style "flapper" valve coated with 58 Ni was lo- 
cated immediately above the SD2 volume. This "flap- 
per" valve opened and subsequently closed (with open- 
ing and closing response times of about 0.1 s) with each 
proton pulse, in order to increase the storage lifetime of 
the UCN in the volume of the source above the SD2 vol- 
ume. A typical UCN lifetime with the flapper open was 
9.6 ± 0.2 s, whereas the lifetime with the flapper closed 
was 39.4 ± 0.1 s. The flapper leads to a corresponding 
increase in the UCN density. 

The UCN were then extracted from the source along 
horizontally-oriented 10.16-cm diameter stainless steel 
guides (presenting a nominal potential of 189 neV) 
through the biological shielding surrounding the source 
and out into the experimental area. As shown in Fig. 2, 
this system of guides through the biological shielding in- 
cluded two 45° bends to eliminate neutrons with kinetic 
energies above the stainless steel guide potential. The 
maximum UCN density at the biological shield exit that 
we have obtained is 52 ± 9 cm~ 3 [53], but for this work 



2 Each proton beam pulse consisted of five 625 £ts beam bursts 
separated by 0.05 s, with 5.2 s between each pulse's leading edge 
burst. 
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(runs in 2008-2009) typical densities were ~ 35 ±6 cm -3 . 
After exiting the biological shield, the UCN were trans- 
ported along stainless steel guides through a gate valve, 
which served to separate the UCN source from the ex- 
periment, thus permitting measurements of backgrounds 
in the electron spectrometer detectors with the proton 
beam still operating in its normal pulsed mode, but no 
accompanying UCN transport to the spectrometer. A 
6.0-T superconducting solenoidal prc-polarizing magnet 
(PPM) was located immediately downstream of this gate 
valve. The PPM was included in the experiment design 
in order to minimize UCN transport losses through a 
thin (0.0254-mm thick) Zr foil which served to separate 
the vacuum in the SD2 source from the downstream vac- 
uum in the remainder of the experiment. Note that the 
UCN population was polarized after transport through 
the PPM's longitudinal magnetic field. 

To preserve this initial polarization, 10.16-cm diameter 
electropolishcd Cu guides (nominal potential of 168 neV) 
were installed downstream of the PPM magnet. The 
UCN were then transported along these guides through 
a "switcher" valve, which allowed the downstream guides 
comprising the /3-asymmetry measurement to be con- 
nected to either the upstream guides from the UCN 
source, or to a 3 Hc UCN detector [56] used, as de- 
scribed later, for measurements of the depolarized pop- 
ulation. These electropolishcd Cu guides then trans- 
ported the UCN through the primary 7.0-T polarizing 
magnet (called the AFP magnet). A 100-cm long quartz 
guide section coated with diamondlike carbon (DLC) [57] 
passed through the center of a resonant (1.0-T) "bird- 
cage" r.f. cavity [58], used for adiabatic fast passage 
(AFP) spin-flipping of the UCN. 

Downstream of this DLC-coated quartz guide, another 
section of 10.16-cm diameter Cu guide transitioned to a 
7-cm (vertical) x 4-cm (horizontal) rectangular Cu guide, 
which transported the UCN through a horizontal pene- 
tration in the 1.0-T solenoidal electron spectrometer coil 
into the decay trap. Permanent magnets were attached 
to the outer surfaces of the rectangular guide in order 
to suppress Majorana spin-reorientations [59] of neutrons 
passing through "field zeros" in the 1.0-T solenoidal spec- 
trometer's field. 

The UCN rate along the transport guide system was 
monitored with 3 He UCN detectors [56] at two key lo- 
cations: at the gate valve (for monitoring of the SD2 
source performance), and slightly downstream of the 
APF spin flipper (for monitoring of the AFP spin- flipping 
efficiency). These 3 He UCN detectors were coupled to the 
guide system via small (0.64-cm diameter) holes in the 
bottom of the guides. Note that the UCN density in the 
spectrometer for the spin-flipped state was smaller than 
that for the unflipped state, because of losses (after the 
2.0-T /i- B energy boost associated with the spin flip) in 
the transport guides located between the AFP spin-flip 
region and the electron spectrometer. The measured /3- 
decay rates for the spin-flipped state were ~ 25% smaller 
than those for the non-spin-flipped state. 



The maximum neutron /3-decay rates measured in the 
spectrometer during the 2008-2009 runs correspond to a 
stored density of approximately 1 cm~ 3 in the decay trap. 
Major sources of loss are transport through the high field 
regions in the PPM and in the AFP magnet. The trans- 
port though the PPM is about 25%. Approximately half 
of the loss is due to polarization of the neutrons, and 
the other half is due to UCN absorption in the Zr foil 
and non-specular scattering on the UCN guide walls in 
the high field region of the magnet which leads to en- 
hanced wall losses. There is an approximate 15% loss in 
UCN density in the transition from the stainless steel to 
the copper guides because of the lower Fermi potential 
of the copper. Transmission through the AFP magnet is 
about 60%, again due to non-specular scattering in the 
high field region. There is a 50% loss in density in load- 
ing the decay trap because the loading time (which is 
determined by the aperture of the above-described 7-cm 
x 4-cm rectangular guide) and decay trap lifetime are 
nearly the same. Finally, there is an approximate fac- 
tor of two loss in the transport from the biological shield 
exit to the decay trap due to guide losses. (The typical 
loss per bounce in the guide system is 3 x 10~ 4 which is 
dominated by gaps in the guide couplings.) Thus, all of 
these factors combined account for the reduction in the 
UCN density from its intial value of ~ 35 cm -3 at the 
biological shield exist to ~ 1 cm -3 in the spectrometer 
decay trap during runs in 2008-2009. 



C. Decay Trap Geometry 

The decay trap consisted of a 300-cm long, 12. 4-cm di- 
ameter electropolished Cu tube situated along the warm 
bore axis of the 1.0-T solenoidal electron spectrometer. 
The vacuum pressure in the UCN guides downstream of 
the Zr foil in the PPM and in the decay trap was typically 
~ 10 ~ 5 Torr. The ends of the decay trap were closed off 
with variable thickness mylar end-cap foils, whose inside 
surfaces were coated with 300 nm of Be (nominal 252 
neV potential) which served to increase the UCN storage 
time in the decay trap (and, hence, the /3-decay rate). 
An additional important feature of the end-cap foils is 
that they eliminated the possibility for neutron /3-decay 
events in the region of the spectrometer where the field is 
expanded from 1.0-T to 0.6-T (discussed later in Section 
HIE 1). Collimators with inner radii of 5.84 cm mounted 
on the two ends of the decay trap suppressed events orig- 
inating near the decay trap walls and also functioned as 
mounts for the end-cap foils. As discussed in more de- 
tail later, the thicknesses of the mylar end-cap foils were 
varied from 0.7 fim to 13.2 /im to study key systemat- 
ics related to electron energy loss in, and backscattering 
from, these foils. 

The UCN density in the decay trap was monitored by a 
3 He UCN detector coupled to a small (0.64-cm diameter) 
hole in the bottom of the decay trap center, as indicated 
schematically in Fig. 2. 



8 



D. Polarization and Spin Flipping System 

A detailed description of the 7.0 T polarizing mag- 
net and the AFP spin- flip system is given elsewhere [58] ; 
therefore, we provide only a brief description of this sys- 
tem here. 

The solenoidal superconducting magnet which serves 
as the primary UCN polarizer (the AFP magnet) and 
provides the requisite environment for an adiabatic fast- 
passage (AFP) spin flipper was designed by American 
Magnetics using a cryostat supplied by Ability Engineer- 
ing. It possesses a 194.9 cm long, 12.7 cm diameter warm 
bore and provides both a peak field of 7.0-T near the en- 
trance as well as a 44.5-cm long precision gradient spin- 
flip region with an average field of 1.0 T, chosen to re- 
duce neutron spectral differences between the flipped and 
unflipped spin states in the electron spectrometer decay 
trap volume. When energized to 96.45 A, the main coil 
of this magnet produces both the maximum polarizing 
field as well as the 1 T field, with an average gradient of 
6 x 10~ 5 T cm -1 through the latter. Ten superconduct- 
ing shim coils centered on the uniform field region and 
spaced every 5.1 cm provide the ability to further tailor 
the uniform field in order to optimize performance of the 
spin flipper. 

Due to the high field in the spin flip region, the spin 
flipper was constructed using an efficient high-pass bird- 
cage resonant cavity geometry [58]. For the UCNA Ex- 
periment, this configuration was realized with eight Cu 
tubes (rungs) arranged in a cylindrical geometry and con- 
nected at the top by 820 pF American Technical Ceram- 
ics chip capacitors and at the bottom by a sliding Cu 
tuning ring whose position determined the inductance 
presented by the Cu tubes. When excited by an r.f. sig- 
nal such a geometry is resonant, with the fundamental 
mode corresponding to a discrctized sinusoidal distribu- 
tion of current in the rungs. This current distribution 
provides a transverse r.f. field, one rotating component 
of which is utilized for AFP spin flipping, providing ef- 
ficient spin reversal over a wide band of neutron speeds 
[58]. The specific operation frequency was adjusted by 
moving the tuning ring, and the cavity formed after tun- 
ing the spin flipper to operate at ~ 28 MHz was ~ 15 cm 
long (8.74 cm diameter), coaxial with a 7 cm diameter 
DLC-coated quartz UCN guide. 

The UCNA spin flipper was typically operated with 
40 W of input power, which necessitated an impedance- 
matching system comprised of a calculated length of drive 
line and three Jennings vacuum variable capacitors. Wa- 
ter cooling was also required, and was accomplished by 
flowing chilled, filtered tap water serially through the 
rungs. In order to provide for stable electrical opera- 
tion and to prevent r.f. radiation from inducing noise 
elsewhere in the experiment, the birdcage cavity was 
driven in a balanced mode and electrically shielded by 
a grounded Al cylindrical enclosure which also provided 
a vacuum seal around the DLC-coated quartz guide. The 
interior of this enclosure connected through four bellows 



to the outside of the AFP magnet so that the actual r.f. 
cavity remained at atmosphere while the AFP magnet 
bore and the guide system were under vacuum. This ar- 
rangement also provided feedthroughs for the r.f. drive 
line, water cooling lines, an RTD temperature sensor, 
and an r.f. field sensor loop. 

Initial characterization of the spin flipper was per- 
formed in a crossed polarizer analyzer geometry as de- 
scribed in [58], which determined the average spin flip 
efficiency to be 0.9985 ± 0.0004. During the actual run- 
ning of the UCNA experiment during the years 2008- 
2009, tuning of the spin flipper, as well as run-to-run 
monitoring of its performance, was accomplished using a 
3 He UCN monitor located just downstream of the AFP 
magnet, ~ 1 m below a ~ 0.64 cm hole in the bottom of 
the guide (the location of this UCN monitor is indicated 
schematically in Fig. 2). This detector had a magnetized 
Fe foil covering the detector acceptance which provided 
for spin state selection. 



E. Electron Spectrometer System 

The electron spectrometer system, consisting of a 1.0- 
T superconducting solenoidal magnet and a multiwire 
proportional chamber (MWPC) and plastic scintillator 
electron detector package, is described in detail elsewhere 
[49, 50]. Nevertheless, for completeness, we discuss the 
primary components of this system here. As described 
below, the electron spectrometer system was designed 
both to suppress the total electron backscattering frac- 
tion and to reconstruct low-energy-deposition backscat- 
tering events. The primary components of the two identi- 
cal MWPC and plastic scintillator detector packages are 
shown in Fig. 3. 



1. 1.0-T Superconducting Solenoidal Magnet 

The spectrometer magnet [50] is a warm-bore 35-cm di- 
ameter, 4.5-m long superconducting solenoid (hereafter, 
SCS magnet). The coil, which was designed and fabri- 
cated by American Magnetics, Inc., consists of a main coil 
winding with a single persistence heater switch, 28 shim 
coil windings (each with individual persistence heater 
switches) , and three rectangular 7-cm x 4-cm radial pen- 
etrations (two providing horizontal access, and one pro- 
viding vertical access, to the warm bore). These penetra- 
tions are located at the center of the coil. The magnet's 
1600-L capacity liquid helium cryostat was designed and 
fabricated by Meyer Tool and Manufacturing, Inc. The 
magnet's full energized field strength of 1.0 T requires 
a current of 124 A in the main coil winding. Note that 
the magnet's shim coils were not energized during the 
data-taking runs reported in this article. 

An important feature of the SCS magnet design was 
that the field is expanded, as indicated schematically in 
Fig. 2, from 1.0 T in the decay trap region to 0.6 T at the 
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FIG. 3. Schematic diagram (not to scale) of the MWPC and plastic scintillator electron detector package. 



location of the MWPC and plastic scintillator electron 
detectors, to suppress large-pitch-angle backscattering. 
In particular, the field-expansion ratio of 0.6 maps pitch 
angles of 90° in the 1.0- T region to pitch angles of 51° 
in the 0.6-T region. Another important feature of the 
magnet design concerned the field uniformity in the decay 
trap region. Electrons emitted with momentum po = 
(Pj_ o +^|o) 1 ^ 2 ' w ith P±.o (P||,o) the initial transverse 
(longitudinal) momentum component, in some local field 
Bq will be reflected from field regions B if B > B cr a = 
(Po/j?l o)-^Oj thus contributing to a false asymmetry. By 
this same process, electrons emitted at large pitch angles 
in the vicinity of a local field minimum will be trapped. 

The SCS field profile measured during the data-taking 
period reported in this article is shown in Fig. 4. As can 
seen there, the field was uniform to the level of 5 x 10~ 3 
over the decay trap length, but included a ~ 0.005 T 
"dip" near the center of the decay trap. Note that the 
field uniformity shown here was degraded from that orig- 
inally published in [50]; this was the result of damage 
to the shim coil persistence heater switches during mul- 
tiple magnet quenches. The impact of this field non- 
uniformity on the measured asymmetry is discussed later 
in this article. 



2. Multiwire Proportional Chamber 

Some of the most important features of the MWPC 
[49] arc as follows. First, because an MWPC is relatively 
insensitive to gamma rays, requiring a coincidence be- 
tween the MWPC and the scintillator greatly suppressed 
gamma-ray backgrounds, a dominant background source 
in previous experiments (see, e.g., background spectrum 
in [45]). 

Second, the MWPC permitted reconstruction of an 
event's transverse (x,y) position. This permitted the 
definition of a fiducial volume and the subsequent rejec- 
tion of events occurring near the decay trap walls. Such 
electrons can scatter from the decay trap walls, leading 
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FIG. 4. Measured SCS field profile during 2008 (circles) and 
2009 (squares). The two data sets have been offset for pur- 
poses of clarity. The coordinate system is such that the 300- 
cm long decay trap is centered at 2 = 0. 



to a distortion in the energy spectrum and/or a bias to 
the asymmetry. The (x, y) position information from the 
MWPC also permitted a characterization of the scintilla- 
tor's position-dependent response, as the scintillator was 
viewed by four photomultiplier tubes (discussed in the 
next section). The 64-wire anode plane was strung with 
10-/im diameter gold-plated tungsten wires, and the two 
cathode planes (oriented at 90° relative to each other) 
were each strung with 64 50-^m diameter gold-plated 
aluminum wires. The wire spacing on both the anode 
and cathode planes was 2.54 mm, yielding an active area 
of 16.3 x 16.3 cm 2 . This area in the 0.6-T field-expansion 
region mapped to a \/0.6(16.3 x 16.3) cm 2 = 12.6 x 12.6 
cm 2 square in the 1.0- T region, thereby providing full 
coverage of the 12.4-cm diameter decay trap volume. As 
demonstrated previously [49, 50], the center of the event 
(i.e., the center of the charge cloud resulting from the 
electron's Larmor spiral in the MWPC gas) could be re- 
constructed with an accuracy of better than 2 mm, suf- 
ficient for the definition of a fiducial volume. 
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Third, to suppress "missed backscattering events" (i.e., 
those events depositing no energy above threshold in 
any detector element along the electron's trajectory prior 
to backscattering), the entrance window separating the 
MWPC fill gas from the spectrometer vacuum was de- 
signed to be as thin as possible. Fourth, because of this 
thin entrance window requirement, the fill gas pressure 
was required to be as low as possible. The chosen fill 
gas, C5H12 (2,2-Dimcthylpropane, or "ncopentane" ) , a 
low-Z heavy hydrocarbon, was shown to yield sufficient 
gain at a pressure of 100 Torr and a bias voltage of 2700 
V. At this pressure, the minimum window thickness (over 
the MWPC's 15-cm diameter entrance and exit windows) 
shown to withstand this 100 Torr pressure differential 
with minimal leaks from pinholes was 6 /jm of aluminizcd 
mylar. Note that the front window was further reinforced 
by Kcvlar fibers. 



3. Plastic Scintillator Detector 

The plastic scintillator detector was a 15-cm diameter, 
3.5-mm thick disk of Eljen Technology EJ-204 scintilla- 
tor. This 15-cm diameter mapped to a 11.6-cm diame- 
ter disc in the 1.0-T decay trap region, providing nearly 
full coverage of the decay trap volume. The range of 
an endpoint energy electron in the plastic was 3.1 mm; 
therefore, the 3.5-mm thickness was sufficient for a mea- 
surement of the full /3-decay energy spectrum, while min- 
imizing the ambient gamma ray background rate. 

With the scintillator located in the 0.6-T field- 
expansion region at a distance of 2.2 m from the center 
of the SCS magnet, light from the disc was transported 
over a distance of ~ 1 m along a scries of UVT light 
guides to photomultiplier tubes which were mounted in a 
region where the magnetic field was ~ 0.03 T. The light 
guide system, shown schematically in Fig. 3, consisted 
of twelve rectangular strips (39-mm wide x 10-mm thick 
UVT) coupled to the edge of the scintillator disc with 
optical grease. These twelve rectangular strips were then 
bent through 90° over a 35-mm radius, transported over 
a distance of ~ 1 m away from the scintillator, and then 
adiabatically transformed into four 39 x 30 mm 2 rectan- 
gular clusters, with 5.08-cm diameter Burle 8850 photo- 
multiplier tubes (PMTs) glued to each of these four rect- 
angular clusters. Therefore, each PMT effectively viewed 
one it/2 quadrant of the scintillator face. 

The magnetic shielding for each of the PMTs consisted 
of an array of active and passive components, including 
(moving from the outside to inside) steel and medium- 
carbon-steel shields, a bucking solcnoidal coil wound on 
the surface of a thin /x-metal foil, and a /i-metal cylinder. 
Magnetic end caps were not required. 

The vacuum housing enclosing the scintillator, light 
guides, and PMTs was maintained at ~ 95 Torr of nitro- 
gen, and was separated from the MWPC volume (with its 
100 Torr of neopentane gas) by the MWPC exit window. 
The nitrogen volume pressure was maintained at a some- 



what lower pressure than the MWPC pressure to ensure 
that the MWPC exit window bowed out, or away, from 
the MWPC interior, to avoid contact with the MWPC 
wire planes. 



F. Scintillator Calibration and PMT Gain 
Monitoring 

The scintillators were calibrated periodically with 
conversion electron sources, including commercially- 
available 109 Cd (63 keV, 84 keV), 139 Ce (127 keV, 160 
keV), 113 Sn (364 keV, 388 keV), 85 Sr (499 keV), and 
207 Bi (481 keV, 975 keV, 1047 keV) conversion-electron 
sources, and a custom-prepared 114m In (162 keV, 186 
keV, 189 keV, 190 keV) conversion-electron source (via 
implantation of 113 In onto an Al substrate and subse- 
quent irradiation in a reactor [60]). These calibrations 
were conducted in-situ using a vacuum load-lock source 
insertion system which permitted insertion and removal 
of calibration sources with the electron spectrometer un- 
der vacuum. The insertion point for these sources was 
through one of the superconducting solenoid magnet's 
horizontal rectangular penetrations at the center of the 
coil. Note that this source insertion system permitted the 
sources to be positioned only along the horizontal axis of 
the decay trap's circular geometry; however, as described 
later, the position dependence of the energy calibration 
over the full circular geometry was achieved by compar- 
ing the reconstructed neutron /3-decay endpoint in a large 
number of binned positions over the scintillator face. 

The PMT gains were monitored on an approximate 
daily basis with a 113 Sn source using this source insertion 
system. Fits to the minimum-ionizing peak of cosmic-ray 
muons served as a run-to-run gain monitor. 

G. Cosmic- Ray Muon Veto System 

The electron spectrometer was surrounded with a 
cosmic-ray muon veto system which consisted of the fol- 
lowing components. First, as shown in Fig. 3, a 15-cm 
diameter, 25-mm thick plastic scintillator (the "back- 
ing veto") was located immediately behind each of the 
spectrometer scintillators. Second, a large-scale plastic 
scintillator and sealed drift tube veto counters [61] sur- 
rounded the electron spectrometer magnet. 

H. Electronics and Data Acquisition 

The frontend electronics for the experiment consisted 
of a VME-bascd system for the event trigger logic (via 
discriminators and programmable logic units (PLUs)) 
and for the readout of scalers, analog-to-digital conver- 
ters (ADCs) and time-to-digital convertors (TDCs). A 
NIM-based system coupled to the VME system was em- 
ployed for the implementation of a "busy logic" , which 
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FIG. 5. Simplified schematic diagram of the data acquisition 
event trigger logic. 



merit of the total charge measured in each PMT. The ana- 
log signals from the cosmic-ray muon backing vetos were 
also read out by these QADC modules. Peak-sensing 
CAEN V785 ADC (PADC) modules, triggered for read- 
out by a ~ 12 [is gate from a CAEN V462 gate generator, 
digitized the MWPC anode and cathode-plane signals. 
Note that the anode signal that was read out was the 
summation (i.e., single channel per anode plane) of the 
signals on all 64 of the wires comprising the anode plane. 
The 64 wires on each of the two cathode planes were 
read-out in groups of four (i.e., 16 channels per cath- 
ode plane); hereafter, we will simply call each of these 
four- wire groups a "wire". Analog signals from the 3 He 
UCN monitors and the drift tube cosmic-ray muon veto 
counters were also read out with PADCs. 



served to veto event triggers arriving during ADC/TDC 
conversion times (i.e., during these modules' busy states). 
This busy logic also prevented re-triggering by correlated 
scintillator afterpulscs (mostly occuring over a ~ 1 [is 
window [62]), and was implemented with a LeCroy 222 
gate generator in latch mode. 

A simplified schematic diagram of the trigger logic is 
shown in Fig. 5. For each detector package, a trigger was 
defined by a two-fold PMT coincidence trigger above the 
discriminator threshold for each PMT (nominally, set at 
0.5 photoelectrons) . The resulting two-fold PMT trigger 
rate in each scintillator was ~ 50 s _1 (primarily from 
low-energy background gamma rays); the singles rates 
in each PMT as determined by counting in scalers were 
typically ~ 500-1000 s _1 (from both dark noise and low- 
energy backgrounds). The main event trigger was then 
defined to be the OR of the two detectors' two-fold PMT 
triggers and other experiment triggers (e.g., from the 
3 He UCN monitor detectors). The logic for the two-fold 
PMT coincidence triggers and the main event trigger was 
performed with CAEN V495 Dual PLUs. Those main 
event triggers not vetoed by the busy logic then trig- 
gered gate/delay generators for the readout of the ADCs, 
TDCs, and scaler modules. The total number of two-fold 
PMT coincidence triggers were counted in scalers as a 
monitor of the DAQ dead time. 

CAEN V775 TDC modules were used for the relative 
measurement of the time-of-flight between the two de- 
tectors' two-fold PMT coincidence triggers. This relative 
timing information provided for the identification of the 
detector with the earlier arriving trigger, important, as 
discussed later, for the assignment of the initial direction 
of incidence for electron backscattering events triggering 
both scintillators. These TDCs were also used to record 
the timing information from the plastic-scintillator-based 
muon veto detectors. A global event-by-event time stamp 
was defined by the counting of a 1 MHz clock in a CAEN 
V830 scaler. 

CAEN V792 charge-integrating ADC (QADC) mod- 
ules, triggered for readout by a ~ 140 ns gate from a 
CAEN V486 gate/delay generator, provided a measure- 



The data acquisition (DAQ) system was based on the 
MIDAS package [63], with a dedicated Linux- based work- 
station for implementation of the frontend electronics 
acquisition code and a separate dedicated Linux-based 
workstation for run control and online analysis. The 
frontend acquisition code accessed the VME crate via a 
Struck PCI/VME interface. The MIDAS raw data banks 
were subsequently decoded into CERNLIB PAW [64] and 
ROOT [65] file formats for data analysis. 

A separate data acquisition system, based on the PCDAQ 
software package [66] was used to monitor the proton 
beam charge incident upon the UCN source's tungsten 
spallation target and to asynchronously monitor environ- 
mental variables in the experimental area. The incident 
proton flux was measured using an integrating current 
toroid mounted around the proton beam line 8 m up- 
stream of the tungsten target, just before the proton 
beam entered the biological shield. As noted earlier in 
Section IIIB, each proton beam pulse consisted of five 
625 [is beam bursts separated by 0.05 s, with 5.2 s be- 
tween each pulse's leading edge burst. The proton charge 
integrating system measured only the charge of the first 
of the five beam bursts in each pulse; the resulting value 
was then scaled by five to yield the total proton charge 
delivered during each of the 0.2 Hz beam pulses. 

The environmental monitoring system asynchronously 
read and stored up to 96 variables, on a typical time scale 
of 0.2 s to 1.0 s between readings. Variables measured in- 
cluded cryogenic temperatures in the UCN source (read 
by Lakeshore 218 temperature monitors), pressures in 
the different segments of the UCN guide system (read by 
capacitance manometer, thermocouple, and cold-cathode 
ion vacuum gauges), ambient temperature in the exper- 
imental area, and liquid helium levels and gas pressures 
throughout the cryogenic systems. The environmental 
data were time-stamped for later comparison to the /3- 
decay data acquired with the main data acquisition sys- 
tem. 
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TABLE II. Run structure for the octet data taking sequence, 
consisting of A- and B-type quartets. See text for details. 



Al 


A2 


A3 


A4 


A5 


A6 


A7 


A8 


A9 


A10 


All 


A12 


B" 


P~ 


D" 


B+ 


P + 


D+ 


P + 


D+ 


B+ 


P~ 


D" 


B" 


Bl 


B2 


B3 


B4 


B5 


B6 


B7 


B8 


B9 


B10 


Bll 


B12 


B+ 


P + 


D+ 


B~ 


P~ 


D" 




D~ 


B~ 


P + 


D+ 


B+ 



IV. MEASUREMENTS, EXPERIMENTAL 
GEOMETRIES, AND POLARIZATION 

In this section we provide a detailed description of our 
measurement procedures for /3-decay and ambient back- 
ground runs; the various geometrical configurations of 
the experiment during our /3-decay runs; and our proce- 
dures for, and results from, measurements of the neutron 
polarization. 



A. /3-Decay Run Cycle 

1. Octet Data-Taking Structure 

The data taking during normal /3-decay production 
running was organized into octets, each consisting of 
A- and B-type quartet run sequences. The structure of 
these quartet and octet run sequences, shown in Table 
II, was such that the neutron spin state (hereafter des- 
ignated + or — , with +(— ) corresponding to the loading 
of UCN with AFP-spin-flippcr-on (-off) spin states into 
the electron spectrometer) was toggled according to a 

— + H 1 h spin-sequence (for octets in which A-type 

runs preceded B-type runs) or a H 1 h H — (i-c, 

complement) spin-sequence, with the order of /3-decay 
and ambient background run pairs toggled for a particu- 
lar spin state within each A-type or B-type run sequence. 
Within each octet, the decision for whether the A-type 
runs would precede or follow the B-type runs was made 
randomly. The notation in Table II is such that B + (~) 
and /3 + ' - ) denote, respectively, ambient background and 
/3-decay runs for the two spin states. The notation for 
depolarization runs is such that D + , for example, denotes 
a measurement of the depolarized spin-state population 
for which the spin-state was polarized in the + spin-state 
during the preceding /3-decay run. 

As described in detail later, the /3-decay yields were 
ultimately obtained from background subtraction. Al- 
though such a procedure is potentially subject to sys- 
tematic bias from time- varying backgrounds, the merit of 
this octet data-taking structure is that linear background 
drifts cancel to all orders (provided that the durations of 
the background and /3-decay runs do not change during 
the octet) in the definition of asymmetries based on com- 
plete octet-structure data sets. Linear drifts in detector 
efficiency which might affect background subtraction also 
cancel under the octet structure. 



2. Run Cycle Procedure 

As noted previously in Section IIIB and shown in 
Fig. 2, a gate valve separated the UCN source from 
the ^-asymmetry experiment. Measurements of the am- 
bient backgrounds (runs Al/Bl, A4/B4, A9/B9, and 
A12/B12) were performed with this gate valve closed 
(i.e., with no UCN in the decay trap), but with the pro- 
ton beam still operating in its normal pulsed mode and 
the AFP spin-flipper in its appropriate run-paired state 
(i.e., so as to properly account for beam-related back- 
grounds and any noise/backgrounds associated with the 
operation of the AFP spin-flipper). These background 
runs were nominally 0.2 hours in duration. The /3-decay 
runs (runs A2/B2, A5/B5, A7/B7, A10/B10) were per- 
formed with the gate valve open and the AFP spin-flipper 
in its appropriate state for the entire duration of the run, 
nominally 1.0 hour in duration. 

During the /3-decay runs, an equilibrium density of 
both correctly polarized and incorrectly polarized UCN 
developed in the decay trap. With the spin flipper off, the 
incorrectly polarized population was dominated by depo- 
larization due to material interactions between the UCN 
and the walls of the decay trap and guides. When the 
spin flipper was active, this incorrectly polarized popula- 
tion was increased as a result of spin flipper inefficiency. 
In the spin flipper off case, the lifetime of correctly polar- 
ized UCN in the decay trap, dominated by the decay trap 
exit aperture, was ~ 21 s, and the lifetime of incorrectly 
polarized UCN trapped in the experimental geometry by 
the 7 T polarizing field was ~ 31 s, dominated by losses in 
the low-field region between the AFP magnet and SCS. 
In the spin flipper on state, the lifetime of correctly po- 
larized UCN was ~ 17 s, and the lifetime of incorrectly 
polarized UCN was ~ 44 s. 

At the immediate conclusion of a /3-decay run, a de- 
polarization run (runs A3/B3, A6/B6, A8/B8, All/Bll) 
was conducted to measure the depolarized fraction of the 
UCN population in the decay trap via the following pro- 
cedure. First, the gate valve was closed, the proton beam 
was gated off, and the switcher valve was re-configured 
such that the guides downstream of this valve (i.e., from 
the switcher valve all the way through the decay trap) 
were connected to a 3 He UCN detector (see the discussion 
in Section IIIB). The state of the AFP spin- flipper was 
unchanged from its state during the immediately preced- 
ing /3-decay run. At this point, UCN of the "correct" spin 
state in the experimental volume could exit the geometry 
through the 7-T polarizing field, which now served as a 
spin-state analyzer, and were counted in the UCN detec- 
tor located at the switcher valve. This cleaning phase 
lasted 25 s, and the number of counts recorded in the 
UCN detector during this time interval was proportional 
to the number of correctly polarized UCN present in the 
experimental geometry 

Following this cleaning phase, the state of the spin- 
flipper was changed. This then permitted those UCN of 
the "wrong" spin state located downstream of the spin- 
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TABLE III. Foil thicknesses for the different decay trap end-cap window and MWPC window Geometries and the number of 
/3-decay events collected in each Geometry passing all analysis cuts. 



Geometry (Year) 


Decay Trap End-Cap 
Window Thickness [/an] 


MWPC Window 
Thickness [/im] 


Number of 
/3-Decay Events 


A (2008) 


0.7 (mylar) + 0.3 (Be) 


25 


5.2 x 10 b 


B (2008) 


13.2 (mylar) + 0.3 (Be) 


25 


5.3 x 10 6 


C (2008) 


0.7 (mylar) + 0.3 (Be) 


6 


2.4 x 10 6 


D (2009) 


0.7 (mylar) + 0.3 (Be) 


6 


1.8 x 10 6 



flipper, which until now had been trapped within this 
volume by the 7-T polarizing field, to exit this volume 
through the 7-T field, and to be counted in the UCN 
detector. This counting during the unloading phase was 
performed for ~ 200 s, which provided for a measurement 
of both the number of wrong spin-state neutrons as well 
as a measurement of the UCN detector background on a 
depolarization run-by-run basis. 

B. Experiment Geometries 

During the data-taking runs in 2008-2009 for the re- 
sults reported in this article, the experiment was op- 
erated in four different geometries with different decay 
trap end-cap and MWPC entrance and exit window foil 
thicknesses, to study key systematic corrections and un- 
certainties related to energy loss in, and backscattering 
from, these foils. The foil thicknesses for these four dif- 
ferent experimental geometries, termed Geometries A, B, 
C, and D, are given in Table III. The number of /3-decay 
events collected in each Geometry passing all of the anal- 
ysis cuts detailed later in this article are also listed there. 

Note that although the foil thicknesses for Geometries 
C and D were identical, wc defined separate geometries 
for these data-taking periods because the UCN transport 
guides in the region between the APF spin-flip region and 
the decay trap (i.e., the circular and rectangular guides, 
see Section III B) were upgraded from (bare) electropol- 
ished Cu in Geometry C to DLC-coated electropolishcd 
Cu in Geometry D. This change to a guide system with a 
higher effective UCN potential in the region downstream 
of the AFP spin-flip region resulted in a different veloc- 
ity spectrum for those UCN stored in the decay trap, the 
details of which were important for the interpretation of 
measurements of the UCN polarization, described below. 



C. Polarization Measurements 

A pair of depolarization measurements for each spin 
state, i.e. a D~ run (following the loading of spin-fiipper- 
off spin states during the preceding /3-decay run) and a 
D + run (following the loading of spin-flipper-on states) 
as described in Section IV A 2, provide, in principle, an 
in situ measurement of the UCN polarization at the end 
of the associated /3-decay interval. This pair of mea- 



surements automatically incorporates all depolarization 
mechanisms including, in the case of flipper-on loading, 
spin flipper inefficiency. Fig. 6 depicts the arrival time 
spectra in the switcher UCN monitor detector and the 
decay trap UCN monitor detector (hereafter, SCS moni- 
tor detector) characteristic of a depolarization measure- 
ment during each of the "load", "clean", and "unload" 
intervals. The states of the gate valve, switcher, and the 
spin flipper during each of these intervals are indicated 
schematically there. 

Determination of the equilibrium polarization at the 
end of the "load" interval (corresponding to a time 
t = 200 s in Fig. 6) was accomplished by using the 
switcher UCN detector to count both the number of cor- 
rectly polarized UCN in the decay trap, & p (t = 200 s), 
during the "clean" interval, and then by changing the 
state of the spin nipper to count the number of incor- 
rectly polarized UCN in the decay trap, &a{t = 225 s), 
during the "unload" interval. From these signals, it is 
possible to extrapolate back to a depolarized fraction 
£ ~ g d ft=200 s) f rom which detection efficiencies cancel 
to first order, and which provides the equilibrium polar- 
ization via P = 1 — 2£. The extrapolation procedure 
requires knowledge of the appropriate storage lifetimes 
for correctly and incorrectly polarized neutrons in the 
system, obtained from the SCS monitor detector. 

Models of the UCN transport confirm the intuitive ex- 
pectation that the time dependence of the "clean" and 
"unload" switcher detector signals is characterized by 
double exponential behavior: the shorter time constant 
is associated with emptying the guide system between 
the switcher detector and the narrow rectangular guide 
to the decay trap, while the longer time constant is as- 
sociated with emptying the decay trap through this rect- 
angular guide. Analysis of the arrival time spectra gen- 
erated by the D± measurements which formed part of 
the beta asymmetry run cycle was accomplished by first 
fitting the switcher detector timing spectrum during the 
clean interval to a double exponential plus background 
(where the background was determined from the last 100 
s of the unload interval) . This established the two ampli- 
tudes and associated time constants A^^~\ t[ + ^' and 

A% \t~2 + ^ which characterize the population of cor- 
rectly polarized UCN in the system at the end of the /3- 
decay (loading) interval, where — corresponds to flippcr- 
off loading and + corresponds to flippcr-on loading. Sim- 
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DT Measurement 



SCS Monitor Detector 



CLOSED/OFF 




FIG. 6. (Color online) Arrival time spectra from the switcher 
UCN monitor detector (two primary figures) and the SCS 
monitor detector (top inset figure) characteristic of depolar- 
ization trapping (DT) and reloaded background (R) measure- 
ments. (The start time t = corresponds to the start of the 
load interval or, conceptually, 200 s prior to the end of a /3- 
decay run.) The states of the gate valve, switcher, and the 
spin flipper during the load, clean, and unload intervals are 
indicated schematically. (With the switcher in the "down" 
state, UCN located between the switcher and the decay trap 
are transported to the switcher UCN monitor detector; in the 
opposite state, UCN are transported from the SD2 source to 
the electron spectrometer.) The timing and states depicted 
here are for the case of spin-flipper-off loading. Flipper-on 
loading measurements reverse the state of the spin flipper. 



ilarly, fitting the unload interval determined the ampli- 
tude ' associated with the smaller time constant 



_(+/- 



and the amplitude A 



(+/-) 



associated with the 



larger time constant t> ' , which ideally characterize 
the population of incorrectly polarized UCN present at 
the end of the cleaning interval. In order to extrapolate 
this population back to the end of the /3-decay measure- 



ment interval, the storage lifetimes r + ; of UCN trapped 
in the system due to their spin state relative the state of 
the spin flipper must also be determined, where here + 
(— ) corresponds to the storage lifetime of UCN whose 
spins are parallel (anti-parallel) to the local magnetic 
field and are thus trapped downstream of the spin flip- 
per when it is off (on) . This was accomplished by fitting 
the unload interval of the SCS monitor timing spectrum 
with a single exponential plus background. Note that 
t + is determined during a flippcr-on loading (D + ) de- 
polarization measurement while r_ is determined during 
a flippcr-off loading depolarization (D - ) measurement. 
Monte Carlo studies indicated that using these storage 
lifetimes to capture the average behavior of the depolar- 
ized population during the clean interval introduced no 
significant bias (at the current level of precision) to the 
extrapolation of this depolarized population back to the 
/3-decay measurement. 

In an ideal depolarization measurement, the cleaning 
interval is made of sufficient length that contributions to 
the signal observed during the unloading interval from 
correctly polarized UCN which are not trapped when the 
spin flipper changes state are negligible. If the number of 
free correctly polarized UCN is large and the depolarized 
signal sufficiently small, however, waiting long enough for 
adequate cleaning can reduce the incorrectly depolarized 
signal to levels below the measurement threshold. Since 
this was the case for the UCNA geometries utilized in 
the 2008-2009 run period, the cleaning time A was set 
to 25 s, just long enough to resolve T\ and Ti. This en- 
hanced the depolarized signal but necessitated separate 
measurements to determine the correctly polarized back- 
ground in the unload timing spectrum, which for this 
clean interval was on the same order as the depolarized 
signal. In particular, depolarized UCN coming from the 
decay volume are expected to appear as part of the A < 
component, but the short cleaning time created a non- 
negligible population of correctly polarized UCN in the 
guides between the spin flipper and the polarizing field 
which are not trapped by the spin flipper and which en- 
ter the decay trap before being detected in the switcher 
detector, causing them to appear as part of A K . 

In order to correct for this reloaded (R) background, 
ex situ measurements, denoted R^~^~\ were performed. 
In these measurements, whose characteristic switcher de- 
tector timing spectrum is shown in Fig. 6, thirteen sec- 
onds prior to the start of the unloading phase the spin 
flipper state was changed for three seconds in order to 
trap an additional reloaded population, which then con- 



tributed to the amplitude A 



(+/-) 



determined from the 



unload phase of the corresponding reload measurement. 
With this additional observable, the reload-corrected po- 
larizations were determined via 
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P~ = 1 - 2 



e A/r +(r -_ r -) {cr(i-t-) 



P+ = 1 - 2- 



A< - j¥LA< 



^ A/T - (T+-T+) (Cl + (l-t + ) 



A+ - ^ + A+ 



a: 



(flipper-off loading) 



(fiippcr-on loading), 



(11) 



where r is a Monte Carlo calculated parameter on the 
order of 0.60 needed to account for the presence of an 
extra population between the spin flipper and the 7.0-T 
region trapped by the three second flipper cycle, 



£(+/-) = e -(A 1 +A 2 )/T 2 (+/ - ) e A 2 /r ( _ / + ) 



(12) 



(where Ai is the length of the flipper cycle and A2 is the 
interval between the end of the flipper cycle and the start 
of the unload phase) is a scaling factor which accounts for 
the evolution of the reloaded population trapped during 
the flipper cycle and corrects for the larger population of 
correctly polarized UCN present to be reloaded during 
the flipper cycle, 3> p is the total number of background- 
subtracted counts recorded during the clean interval, (2 
is a factor which uses Ax, t\, A2, and T2 to extrapolate 
3> v to the number of counts which would be observed 
for an infinitely long cleaning period, and JV is a nor- 
malization factor. Values of P + and P~ were obtained 
separately for the 2008 and 2009 data sets by summing 
all corresponding D and R runs and applying Eq. (11). 
Since there was no statistically significant difference at 
the la level between any of the four measurements, a 
single rcload-corrected value P for the polarization was 
obtained by performing a weighted average over the four 
measurements . 

Spin flipper inefficiency decreases the UCN polariza- 
tion for flipper-on loading, resulting in the expectation 
that P~ > P + . The resulting decreased polarization for 
the case of flipper-on loading due to the spin flipper inef- 
ficiency is the actual polarization of the UCN population 
stored in the decay trap, and no further correction to 
the P + value is required. However, the spin flipper in- 
efficiency also leads to a (smaller) increase in P~ since 
correctly polarized UCN which should remain trapped 
during the unloading phase are freed when they are not 
flipped, adding to the observed unload signal. Since this 
population is generated after the /3-decay measurement 
interval it requires a correction, which will decrease the 
value of P~ . The accumulated data limited this correc- 
tion to be no larger than ~0.15% of the total polariza- 
tion, and error bars on P were expanded accordingly. It is 
also possible to have depolarized UCN populations whose 
storage lifetime in the system is much shorter than the de- 
polarization measurement time, and which therefore have 
a low efficiency for detection in a depolarization mea- 
surement. Neutrons with sufficient energy to surmount 
the potential barrier presented by the 7.0-T polarizing 
field (which may therefore enter the experiment in the 



I 

wrong spin state) and initially polarized UCN with ener- 
gies higher than the material potential of the decay trap 
walls (which can survive in the system when confined to 
trajectories that sample the walls at sufficiently oblique 
angles) are examples of such populations. Monte Carlo 
calculations estimating the effect of these populations on 
the neutron polarization indicated a negligible contribu- 
tion at the current level of precision, due largely to the 
short residency times that such UCN posses. Expanding 
the error bars on P to account also for variations in the 
polarization due to pulsed loading (estimated to be on 
the order of 0.04% of P), a la lower limit of P > 0.9948 
was determined. 



V. CALIBRATION AND RECONSTRUCTION 

We now turn to a discussion of our data analysis and 
energy calibration procedures. We begin by defining the 
various possible event types in the experiment, the selec- 
tion rules for the observable event types, and a descip- 
tion of our position recontruction algorithm using the 
MWPC signals. We also discuss our data "blinding" 
procedure, which ultimately resulted in "blinded" asym- 
metries which were scaled by a randomly chosen scaling 
factor at the O(0-5%) level. Next, we discuss our en- 
ergy calibration procedures for the scintillator and the 
MWPC, and then compare for the different event types 
our reconstructed energy spectra with simulated Monte 
Carlo spectra. Finally, we conclude this section with a 
discussion of our procedure for the assignment of the ini- 
tial energy of the electron. 

Hereafter, in our discussions of the data analysis of 
the detector signals, we will refer to the two electron 
detectors as the "East" and the "West" detectors, corre- 
sponding to their actual physical locations in the UCNA 
Experiment. 



A. Event Type Definitions 

Measurement of the /3-asymmetry requires an accu- 
rate determination of the decay electron's initial direc- 
tion of incidence. This determination is complicated by 
backscattcring effects, some of which are not detectable. 
We define the various classes of event types, shown in 
Fig. 7. 
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FIG. 7. Schematic diagram of the various event types defined 
in the text. 

• No backscattering events: Events in which an elec- 
tron, incident initially on one of the detectors, does 
not backscatter from any element of that detector, 
and then generates a two-fold PMT trigger in that 
side's scintillator. 

• Type 1 backscattering events: Events in which an 
electron, incident initially on one of the detectors, 
generates a two-fold PMT trigger in that side's 
scintillator, backscatters from that scintillator, and 
then generates a two-fold PMT trigger in the op- 
posite side's scintillator. Note that a measurement 
of the relative time-of-flight between the two scin- 
tillators' two-fold PMT triggers determines the ini- 
tial direction of incidence for Type 1 backscattering 
events. 

• Type 2 backscattering events: Events in which 
an electron, incident initially on one of the de- 
tectors, deposits energy above threshold in that 
side's MWPC, backscatters from some element of 
that side's MWPC (e.g., gas, wire planes, or the 
exit window) or the scintillator without generating 
a two-fold PMT trigger (e.g., backscattering from 
the scintillator's dead layer, or triggering only one 
PMT), and is then detected in the opposite side's 
scintillator. As can be seen in Fig. 7, the initial 
direction of incidence of a Type 2 backscattering 
event would be misidentified using only scintillator 
two-fold PMT trigger information. 

• Type 3 backscattering events: Events in which an 
electron, incident initially on one of the detectors, 
generates a two-fold PMT trigger in that side's scin- 
tillator, backscatters from that side's scintillator, 
deposits energy above threshold in the opposite 
side's MWPC, and is then stopped in some ele- 
ment of the opposite side's MWPC or scintillator 
without generating a two-fold PMT trigger (i.e., in 
the dead layer, or triggering only one PMT). Note 
that Type 2 and Type 3 backscattering events can- 
not be distinguished using only scintillator two-fold 
PMT trigger information and a threshold cut on 



the MWPC response (i.e., the Type 2 and Type 3 
events depicted in Fig. 7 would not be distinguish- 
able) . 

• Missed backscattering events: Events in which an 
electron, incident initially on one of the detectors, 
backscatters from cither the decay trap end-cap foil 
or that side's MWPC without depositing energy 
above threshold (e.g., from the entrance window, 
or from the gas in the region between the cathode 
plane and the entrance window), and is then de- 
tected in the opposite side's MWPC and scintilla- 
tor. Note that missed backscattering events cannot 
be identified experimentally. 

• Lost events: Events in which an electron, incident 
initially on one of the detectors, deposits signifi- 
cant energy in a decay trap end-cap foil and/or the 
MWPC, and does not generate a two-fold PMT 
trigger in cither of the scintillators. Note that be- 
cause these events do not generate a DAQ event 
trigger, they cannot be identified experimentally. 

No backscattering events and Missed backscattering 
events cannot be distinguished experimentally, and are, 
hereafter, termed "Type 0" events. Based on scintilla- 
tor information alone, Type 2 backscattering events can- 
not be distinguished from Type 3 backscattering events. 
Thus, we will refer to these types of events as "Type 2/3" 
events. Later in Section VI A 2 we will discuss the sepa- 
ration of Type 2/3 events using MWPC information and 
simulation input. Finally, Lost events cannot, of course, 
be reconstructed and can only be corrected for in simu- 
lation. 



B. Run Selection 

Proton beam delivery constraints and other experi- 
mental issues prevented on occasion the accumulation of 
complete octet data sets during normal /3-decay produc- 
tion running. In the absence of a complete octet, runs 
forming a quartet (i.e., runs A1-A12 or B1-B12 in Ta- 
ble II) or spin-pair (i.e., A1-A6, A7-A12, B1-B6, or B7- 
B12) were retained for analysis. Runs with clear detector 
issues (e.g., noisy channels associated with the MWPC 
cathode planes) were discarded. 

C. Data Blinding 

We performed a blinded analysis of our asymmetry 
data by applying separate spin-dependent randomly cho- 
sen scaling factors to the two detectors' count rates, thus 
effectively adding an unknown scaling factor to the mea- 
sured asymmetry. This was implemented via the fol- 
lowing procedure. First, note that the detector count 
rates were based on a global event-by-event clock time 
which was defined, as described earlier in Section III H, 
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by counting a 1 MHz clock in a scaler. Second, we gen- 
erated two random scale factors, f\ and which were 
constrained to be between 1.00 ± (0.05 x 0.04), where 
0.04 represents the approximate value of the measured 
asymmetry. 

For runs with the AFP spin- flipper on (+ spin state), 
we then scaled the east and west detector clock times, t\ 
and tyy, by these scale factors according to 

t + E = t-h, 4r = t-f*> ( 13 ) 

where t denotes the true global clock time. Similarly, for 
those runs with the AFP spin- flipper off (— spin state), 

t E =t-f 2l tw=t-h. (14) 

In calculating the super-ratio for a spin-state run pair 
according to Eq. (10), the resulting blinded measured 
asymmetry, Abimd, is then a function of the blinded and 
true super-ratios, i?biind and i?true, according to 

^ _ 1 — V -Rblind _ 1 — /v -Rtruc 

1 + V -Rblind 1 + fyRtrue 

where / = / 2 //i- 

All of our analysis was performed with these blinded 
asymmetries. Note that this did not impact our assess- 
ment of our energy reconstruction algorithms (which do 
not, of course, depend on the asymmetry), or our assess- 
ment of our systematic corrections for backscattcring and 
the cos ^-dependence of the acceptance, as these were cal- 
culated (and subsequently benchmarked) in units of A . 



D. Data Quality Cuts and Live Time Definition 

We subjected each run to a number of so-called 
"global" and "event-by-event" data quality cuts, result- 
ing in the removal of either a consecutive range of events 
or a single event, respectively. First, we note that there 
were sporadic corruptions to the data stream resulting 
from malfunctioning DAQ electronics modules. These 
electronics problems resulted in cither the corruption of 
all subsequent events following the occurence of the prob- 
lem, or the corruption of only a single isolated event. 
Electronics problems resulting in the corruption of all 
subsequent events included misalignments of the VME 
data banks (e.g., of the QADC data bank relative to the 
PADC data bank) and sudden shifts in the TDC channel 
peak positions of the two detectors' two-fold PMT self- 
trigger timing peaks. After the identification of either of 
these problems, a global data quality cut was applied, re- 
sulting in the removal of all subsequent events in that run. 
Electronics problems which resulted in the corruption of 
only a single event included corruptions to the headers 
and/or footers of the PADC, QADC, or TDC event banks 
(e.g., an event lacking a header or footer) and corruptions 
to the TDC bank event counter relative to the MIDAS data 




Beam Burst Scintillator Trigger Time [s] 

FIG. 8. Typical scintillator two-fold PMT trigger rates during 
the 0.2 Hz proton pulse repetition cycle. For the clock timing 
shown here, each of the five beam bursts comprising each 
proton pulse occur at t — 0, as the clock was reset to t = 
with each beam burst; thus, under this reset scheme, all of the 
events collected during the five beam bursts appear between 
and 0.05 s (as noted in Section IIIB, the time between beam 
bursts was 0.05 s). The dashed line indicates the cut at 0.05 s. 



acquisition event counter. An event- by- event data qual- 
ity cut was then applied to those events found to have 
cither of these latter two types of electronics problems. 

A significant fraction of the data acquired during the 
Geometries A and B running was discarded due to the 
above-described electronics problems. In particular, of 
the data acquired during the second half of the Geometry 
A running and then during the entire Geometry B run- 
ning, ~ 5% of the events lacked a header and/or footer 
and up to ~ 30% of the events suffered from the TDC 
bank event counter problem. In contrast, the fraction of 
events suffering from electronics problems acquired dur- 
ing the Geometries C and D running was small, with 
< 5% of the data affected by these problems. 

Other global data quality cuts included the removal 
of all events between 0.00 s to 0.05 s after each proton 
beam burst. Typical scintillator two- fold PMT trigger 
rates during the 0.2 Hz proton beam beam pulse repeti- 
tion cycle are shown in Fig. 8. As can be seen there, the 
peak scintillator trigger rates were up to a factor of ~ 80 
higher during the proton beam bursts. This figure illus- 
trates one of the merits of a pulsed-spallation-source of 
UCN, namely, that the experiment can be performed in a 
low background environment (i.e., ambient backgrounds 
only) during the time between the proton beam bursts. 
Another global data quality cut included the removal of 
events from /3-decay runs occuring during time periods 
when the rate on the 3 He UCN monitor detector located 
near the gate valve dropped below some threshold were 
vetoed, so as not to degrade the signal-to-background 
ratio. 

After application of the above-described global data 
quality cuts, we then computed on a run- by-run basis a 
"live time" for each detector, defined to be the sum of the 
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(blinded) clock times of the run segments surviving the 
above-described global data-quality cuts. Specifically, if 
a run segment between event M and event N > M sur- 
vived these global data-quality cuts, the corresponding 
east and west detector live times, Tg and Tyy, for this 
run segment were calculated as 



rp J.N + 



E > 



Tw — ty/ — t 



,M 

<>W> 



(16) 



where (tyy) an d t% i^w) denote the blinded time 
stamps for the east (west) two-fold PMT trigger for 
events M and TV, respectively, defined previously in Eqs. 
(13) and (14). 

We note, however, an exception to the above-described 
procedure under which we applied an event- by- event data 
quality cut to events with a TDC event counter problem 
(with no subsequent correction to the live time). In par- 
ticular, as was already noted, up to ~ 30% of the events 
collected during the Geometry B running suffered from 
this problem. Further, the fraction of corrupted events 
recorded by each detector differed for the two neutron 
spin states, thus biasing the extracted asymmetry. Thus, 
it was necessary to correct the Geometry B live times. 
As discussed later in Section VIII F, we corrected for this 
Geometry B live time problem using background gamma- 
ray events, which were uncorrelated with the neutron [3- 
decay events. The correction factors to the live times 
were then defined for each detector on a run-by-run ba- 
sis to be the ratio of the number of gamma-ray events 
surviving the cvent-by-event TDC event counter cut to 
the total number of recorded gamma-ray events. 



E. Event Reconstruction and Identification 

Those events surviving the above data-quality checks 
were then reconstructed and identified (i.e., tagged as a 
Type 0, Type 1, or Type 2/3 event) based on the TDC 
measurement of the two detectors' two-fold PMT coinci- 
dence trigger time-of- flight, the pulse height (PADC) in 
the MWPCs, and the pulse height (PADC and QADC) 
and timing (TDC) information from the various muon- 
veto detectors. The available detector information is de- 
scribed in more detail below. 



1. Scintillator Timing Information 
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FIG. 9. Typical two-fold coincidence trigger relative time- 
of-flight spectrum for one of the scintillator detectors from a 
/3-decay run. See text for details. 



main event triggers generated by a East detector two-fold 
trigger, with no later arrival of a West detector two-fold 
trigger, appear at the peak at ns (i.e., a "time-out" 
for that TDC channel). Those events forming the broad 
secondary peak from ~ 20-140 ns, separated from the 
self-trigger peak by the dashed line, correspond to main 
event triggers generated by a East detector two-fold trig- 
ger, with the later arrival of a West detector two-fold trig- 
ger. That is, Type 1 backscattcring events incident ini- 
tially on the (opposite-side) East detector comprise this 
broad peak. The backscattering cut line (the dashed line 
in Fig. 9 at 140 ns) placed approximately 15 ns before 
the self-timing peak is justified by the fact that the min- 
imum time-of- flight (i.e., straight-line trajectory) for an 
800 keV electron (/3 = 0.919) traveling the 4.4 m distance 
between the two scintillators is 15.9 ns. 

Note that for the TDC dynamic range setting of ~ 140 
ns shown in Fig. 9, a "true" Type 1 event with a coin- 
cidence time-of-flight greater than ~ 140 ns would not 
appear in the opposite-side detector's coincidence timing 
spectrum, but would instead appear in its "time-out" 
peak; therefore, such a "true" Type 1 event would be 
misidentified in data analysis as a Type 2/3 event. How- 
ever, as calculated in simulation, the fraction of Type 
1 events with a time-of-flight greater than ~ 140 ns is 
small, ~ 2%, and results in a negligible effect on the 
asymmetry. 



The logic outputs from the PLUs defining the two-fold 
PMT triggers for the two detectors were routed to TDC 
channels, forming the STARTS for their respective chan- 
nels. A copy of the main event trigger was delayed by 
~ 155 ns and was used to generate the COMMON STOP for 
the TDC. A typical two- fold coincidence trigger time-of- 
flight spectrum for one of the detectors (for example, of 
the West detector) is shown in Fig. 9. Main event trig- 
gers generated by a West detector two-fold trigger ap- 
pear at the self-trigger delay time of ~ 155 ns, whereas 



2. MWPC Spectra and Particle Identification 

As noted earlier, the MWPC anode and cathode pulse 
height signals were read out on PADC channels during a 
~ 12 [is window after each event trigger. Typical anode 
and cathode (summed over all of the individual cathode 
channels) pulse height spectra during a /3-decay run are 
shown in Fig. 10. For scintillator event triggers, the ma- 
jority of the events appeared in the pedestal, and were 
tagged as gamma ray events. Those events satisfying a 
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FIG. 10. Typical MWPC anode and summed cathode pulse 
height spectra from a /3-decay run. The dashed lines indicate 
the cut positions. The peak in the anode spectrum at channel 
~ 2250 is overflow. Because the cathode sum is obtained via 
summation over all of the individual cathode channels, chan- 
nel overflows are distributed throughout the summed spec- 
trum. 
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FIG. 11. Typical scintillator spectrum for a background run, 
without (light shaded) and with (dark shaded) the require- 
ment of a MWPC coincidence cut. See the text for a descrip- 
tion of the features in the resulting concidence spectrum. 



TABLE IV. Detector selection criteria (scintillator two-fold 
PMT TDC trigger timing and MWPC pulse height) for the 
identification of Type 0, Type 1, and Type 2/3 events. A *J 
indicates the requirement of a valid detector hit, whereas a x 
indicates the requirement of no detector hit. The notation (*) 
indicates the requirement of an earlier detector hit relative to 
the opposite-side detector hit. 



Event Type 


East Detector 


West Detector 


TDC 


MWPC 


TDC 


MWPC 


Type East Incidence 
Type West Incidence 


V 

X 


V 

X 


X 

V 


X 

V 


Type 1 East Incidence 
Type 1 West Incidence 


V (*) 
V 


V 


V 


V 
V 


Type 2/3 East Trigger 
Type 2/3 West Trigger 


V 

X 


V 


X 

V 


V 



cut on either the anode PADC channel number or the 
summed cathode PADC channel number, indicated by 
the dashed lines there, were identified as charged particles 
(electrons or muons). Electron hits in a particular scin- 
tillator were then further separated from muon hits by 
the requirement of no coincident hits in any of the same- 
side muon veto detectors. Note that the coincidence win- 
dow for the plastic scintillator (drift tube) muon veto de- 
tectors was defined by the ~ 140 ns (~ 12 fis) QADC 
(PADC) data acquisition gate. 

After removal of the events occuring during the pulsed 
proton beam bursts, the requirement of a MWPC- 
scintillator coincidence cut further reduced the total in- 
tegrated ambient background rate in each scintillator by 
a factor of ~ 40 from ~ 50 s _1 to ~ 1.2 s _1 over the 
range of QADC channels corresponding to the neutron 
/3-decay energy spectrum 3 . This is illustrated in Fig. 11, 



where the (uncalibrated) scintillator spectrum in QADC 
channels for a typical background run is shown with and 
without the requirement of a MWPC-scintillator coin- 
cidence for a typical background run. The resulting 
broad peak centered approximately at channel 5000 is 
from minimum-ionizing cosmic-ray muons. The origin 
of the surviving low-energy tail (channels < 2000) was 
demonstrated in Monte Carlo studies to be the result of 
cosmic-ray muon interactions with spectrometer mate- 
rials (e.g., in the SCS magnet's infrastructure) located 
near the MWPC and scintillator detector, such as 5- 
ray production and Compton scattering of muon-induced 
gamma rays. 



3 The total background rate without the MWPC coincidence cut, 



and thus the background suppression factor under the coinci- 
dence cut, varies strongly with the scintillator threshold. The 
numbers quoted above are from a typical background run. 
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3. Event Type Classification 

Electron events were classified as Type 0, Type 1, or 
Type 2/3 events according to the detector selection cri- 
teria listed in Table IV. Electron events not satisfying 
any of these criteria were discarded. 

For Type events, the assignment of the initial direc- 
tion of incidence (i.e., on the East or West detector) was, 
of course, trivial. For Type 1 events, the initial direction 
was determined from the TDC two-fold trigger time-of- 
flight spectra, as discussed earlier in Section VE1. At 
this stage of the analysis, Type 2/3 events were iden- 
tified as such but were only tagged with the triggering 
scintillator side (i.e., not yet assigned an initial direction 
of incidence) . 

F. MWPC Position Reconstruction 

1. Algorithm 

The transverse (x, y) position of an event was recon- 
structed from the charge distribution on the MWPC 
cathode plane wires. The position reconstruction algo- 
rithm that was employed for the calculation of the po- 
sition from the signals that were digitized by the PADC 
modules began by identifying, on each of the x- and y- 
planes, the wire with the maximum PADC readout value 
above pedestal. The position was then reconstructed 
from the PADC reading of this "maximum wire" and 
those of the immediately adjacent wires, with all other 
wires ignored, by fitting these three wires' PADC read- 
ings as a function of their x- or y-positions to a Gaussian 
shape 4 . The positions for ~ 90% of the neutron /3-decay 
events could be reconstructed with this algorithm. 

The positions for the remaining ~ 10% of the events 
were reconstructed under various special circumstances. 
For example, if there was only one wire with a PADC 
value above pedestal, the position was assigned to that 
particular wire's coordinate. Or if only one of the wires 
directly adjacent to the maximum wire recorded a PADC 
value above pedestal, or if the maximum wire was located 
on one of the cathode plane edges and the adjacent wire 
recorded a PADC value above pedestal, the position was 



4 Although the distribution of charge across the cathode planes is 
technically not Gaussian (for the idealized case of a charge lo- 
cated above an infinite grounded conducting cathode plane, the 
shape of the induced charge distribution on the cathode plane 
would be equivalent to that obtained via the method of images; 
a more realistic model for the finite-length wires on the finite- 
size cathode plane would require a finite-element analysis cal- 
culation), employing a Gaussian fit was taken to be sufficient 
given that: (1) small-Larmor-radius calibration source spots 
(e.g., 139 Ce) reconstructed to correctly-sized spots, and were 
highly repeatable; and (2) the distribution of /3-decay events, 
as can be seen in Fig. 12, is fairly uniform (as expected), which 
then sets a limit on any positioning reconstruction errors. 



determined under the assumption that the width of the 
charge distribution was 0.75 times the wire spacing (as 
determined from the ~ 90% of the events reconstructed 
with the three- wire Gaussian fitting algorithm). In the 
event two or more wires recorded overflow PADC values, 
the positions were reconstructed to sharp peaks halfway 
between the two overflow wires' coordinates. 

Finally, if the x- or y-position could not be recon- 
structed (e.g., if none of the wires on one or both the 
cathode planes recorded a PADC value above pedestal), 
the x- or y-position was defined to be 0.0. This was po- 
tentially the source of a systematic bias, because as is 
discussed in detail later in Section VH, the scintillator 
energy reconstruction was position dependent. However, 
the fraction of events identified as charged particles fail- 
ing either the x- or ^-reconstruction was small, < 10 -3 , 
and is consistent with our later estimate (Section VIII I) 
of the MWPC efficiency The origin of these events may 
have been electronic noise on the MWPC anode channel 
correlated with a scintillator trigger, thus resulting in the 
identification of the event as a charged particle with no 
corresponding signal on the cathode PADC channels. 

Note that all of the (x, y) position spectra shown here- 
after have been projected back to the spectrometer's 1.0- 
T field region from the 0.6-T field-expansion region. 



2. Position Cuts, Fiducial Cut, and Coordinate Systems 

The reconstructed positions were used to define a fidu- 
cial volume in order to reject events originating near 
the decay trap collimator (inner radius of 58.4 mm). 
Examples of such problematic events include both j3- 
decay electrons (subject to cos ^-dependent acceptance 
and backscattcring effects) and background electrons 
from the Compton scattering of gamma rays in the decay 
trap material. 

For Type events, we required the reconstructed ra- 
dius on the triggering scintillator side to satisfy a con- 
servative rigger < 45 mm cut. The maximum Larmor 
diameter for an endpoint energy electron in the 1.0-T 
spectrometer field was 7.76 mm; thus, this 45 mm cut was 
safely located ~ 2 Larmor diameters from the decay trap 
wall. For events identified as Type 1 or 2/3 backscatter- 
ing events, in addition to this rt r igg er < 45 mm trigger- 
side cut, we also made a further cut on the two sides' 
i£ - x*w vertex difference to eliminate accidental back- 
grounds, with this cut defined to be \xe — xw\ < 25 mm. 

To study possible systematic effects associated with 
the position reconstruction, we considered four differ- 
ent coordinate systems. (1) Coordinates defined by 
their nominal positions reconstructed from the just- 
described Gaussian fits to the wires' PADC readings: 
xe = {xetVe) and xw = {xw,Vw)- (2) West-side co- 
ordinates defined by their nominal positions, xw, but 
transformed East-side coordinates x' E = {x' E ,y' E ) based 
on fits to Type 1 backscattering data which provided 
for a detector-to-detector coordinate system based on 
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FIG. 12. Distributions of (x, y) position distributions [for the coordinate system choice (4) described in the text] for Type 
neutron /3-decay events (left panel: East detector, right panel: West detector). The dashed circles denote the nominal 45 mm 
fiducial volume radius cut. The regularly-spaced "spikes" are an artifact of the position reconstruction algorithm: for events 
in which two adjacent cathode channels recorded overflow values, the algorithm reconstructed the positions to sharp peaks 
halfway between the two wires' positions (as discussed in Section VF1). 



these backscattering measurements of the detector-pixel- 
to-detector-pixel magnetic field map. These fits were of 
the form x' E = IZx E + c, where 1Z denotes a rotation ma- 
trix, and c = (c x , c y ) is a constant offset vector. (3) Coor- 
dinates defined by the complement of (2) (i.e., East-side 
coordinates defined by their nominal positions, but West- 
side coordinates transformed according to the Type 1 
backscattering detector-pixel-to-detector-pixel magnetic 
field map) . (4) A "split" coordinate system, in which the 
rotation matrix 1Z was applied to one set of coordinates, 
with the constant offset vector c split between the two de- 
tectors, such that x' E = lZx E + c/2, and x' w = —c/2. The 
latter was chosen as the default coordinate system in the 
analysis presented hereafter. We assess the systematic 
uncertainty associated with our choice of a coordinate 
system in Section VIII D. 

Sample (x, y) position distributions [calculated for the 
coordinate system choice (4) described in the above para- 
graph] extracted from the East and West MWPC for 
Type neutron /3-decay events (after background sub- 
traction) are shown in Fig. 12. The distribution of back- 
ground events was previously shown to be uniform over 
the MWPCs' sensitive areas [50]. 



G. Monte Carlo Simulation Programs 

We developed two independent Monte Carlo simu- 
lations of the experimental acceptance, based on the 
GEANT4 (version 9.2) [67] and PENELOPE (version 3) [68] 
simulation codes. These were employed extensively as 
input to our energy calibration procedures (described 
next) and for calculations of the systematic corrections 
for backscattering and the cos ^-dependence of the accep- 
tance (e.g., from suppression of the acceptance at large 
angles), discussed in detail in Section VII. The perfor- 



mance of these simulation programs was benchmarked 
previously in a series of measurements of backscattering 
from beryllium, plastic scintillator, and silicon targets 
and found to be accurate to within 30% in their pre- 
dictions for the angular and energy distributions of the 
backscattered electrons [62, 69-71]. 

Both of these simulations included detailed geometric 
models for: (1) the 3-m long decay trap and its end-cap 
foil geometry; (2) the MWPCs, including their entrance 
and exit windows, the Kevlar fiber support for the en- 
trance windows, the anode and cathode wire planes, and 
the 100 Torr neopentane fill gas; (3) the two dead regions 
in the MWPC (i.e., the two regions between the cathode 
planes and the entrance/exit windows); (4) the plastic 
scintillator disc; (5) the measured scintillator dead layer 
of 3.0- (im thickness [62]; and (6) the magnetic field in 
the decay trap and field expansion regions. Note that 
depending on the study of interest, the magnetic field 
was modeled either as a perfectly uniform 1.0-T field in 
the decay trap region analytically connected (subject to 
the V • B = requirement) to the 0.6-T field in the field- 
expansion region, or via bicubic spline interpolation of a 
three-dimensional grid of the field profile (e.g., such as 
for the measured field profile shown in Fig. 4). In both 
the GEANT4 and PENELOPE simulations, charged particles 
were transported through the magnetic field via Rungc- 
Kutta integration of the equations of motion. 

For simulation studies of calibration sources, events 
were generated isotropically into 4-7T from a fixed point, 
with a model for the source foil covers (assumed to be 
3.6-/im mylar). For neutron /3-decay simulations, events 
were generated uniformly over the decay volume. The 
GEANT4 events were emitted isotropically into An and 
then weighted with a W(9) oc (1 + Aft cos 9) weight fac- 
tor, whereas the PENELOPE events were sampled from the 
full phase-space distribution of Eq. (4). 
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H. Visible Energy Calibration 

1. Overview 

The overarching goal of the energy calibration pro- 
cedure was to calibrate the quantity of scintillation 
light produced by an electron which deposited a certain 
amount of "visible energy" in the scintillator and to cal- 
ibrate the electron energy deposition in the MWPC fill 
gas. We define the scintillator visible energy, E v i S , to be 
the total energy loss in the scintillator active region (i.e., 
beyond the dead layer) measured by the photomultiplier 
tubes. The scintillator visible energy is, of course, not 
equal to the initial energy of the emitted /3-decay elec- 
tron, due to reconstructable energy loss in one or more of 
the MWPCs' active regions, and non-reconstructable en- 
ergy loss in one or more of the decay trap end-cap foils, 
one or more of the MWPCs' non-active elements (e.g., 
windows, wire planes, and gas region between the cath- 
ode planes and the windows), and in one or more of the 
scintillators' dead layers. The relationship between the 
measured visible energy in the scintillator, the measured 
energy deposition in the MWPC, and the reconstructed 
"true" initial energy of the /3-decay electron, denoted 
-Brecon j was determined via comparison of conversion elec- 
tron source measurements to Monte Carlo studies (as dis- 
cussed in Section VN). 

As described earlier in Section HIE 3, each PMT effec- 
tively viewed one 7r/2 quadrant of the scintillator; there- 
fore, the response of each PMT was expected to vary as 
a function of the (x, y) position of the event. Each PMT 
then produced a signal which was a nonlinear function of 
the scintillation light it viewed, with the nonlinearity the 
result of known physics (i.e., quenching interactions in 
the scintillator [72]) and also hardware response issues. 
The scintillation light was then ultimately digitized by 
the QADC data acquisition modules. Accounting for a 
time-dependent gain of the PMT/QADC system, the to- 
tal response of the system to PMT i in terms of digitized 
QADC channels was then modeled as 

QADC i = g i (t)-f i (r H (x,y)-E Q ), (17) 

where Eq denotes the light produced in the scintilla- 
tor, rji(x 7 y) is the fraction of that light reaching PMT 
i from the position (x,y) of the event, fi is the (possibly 
non-linear) response function of PMT i, gi(t) denotes the 
time-dependent gain fluctuation of PMT i, and QADCj 
is the digitized readout for the event. Note that we use 
the notation Eq for the produced scintillation light be- 
cause the produced light should be proportional to the 
"quenched energy" , which we define to be the true energy 
deposition reduced by a Monte-Carlo-calculated quench- 
ing factor (based on studies with a low-energy electron 
gun [62, 73]). The QADC readout signal that is measured 
also includes statistical fluctuations due to photoclectron 
(PE) counting statistics, with the number of PEs propor- 
tional to the scintillation light transported to the PMT, 



Ni oc r)i(x,y)EQ, with the fluctuations in iVj expected to 
follow a Poisson distribution. 

The scintillator energy calibration procedure was thus 
divided into three primary tasks: (1) a determination of 
the linearity function /, for each PMT; (2) a determi- 
nation of the light-transport efficiency r]i{x,y) position 
response map for each PMT; and (3) a determination of 
the time-dependent gain gi(t) for each PMT. Detailed de- 
scriptions of each of these tasks follow below. However, 
note that the overall calibration procedure was highly it- 
erative, whereby previous approximate results from the 
other two tasks were used to produce new, more refined 
results for the task in question. Therefore, our discus- 
sions below of the individual tasks reference input from 
the other two tasks. 



2. PMT Linearity Functions 

The fi response functions for each of the PMTs were 
determined each time a new set of calibrations was per- 
formed at multiple (x, y) positions with the conversion 
electron sources using the remote source insertion system 
described previously in Section IIIF. Most of the useful 
calibration data were obtained with the 139 Ce, 113 Sn, and 
207 Bi sources; however, note that the 139 Ce source was 
not available during data taking for Geometries A and B. 
The 109 Cd source was visible only in Geometries C and 
D (configuration with 6-^m MWPC windows). 

The QADC spectra for each PMT were then recorded 
for each (x, y) source position. The expected Eq spec- 
trum for each PMT for each (x, y) source position was 
then determined from the Monte Carlo simulation pro- 
grams, with the Monte Carlo spectrum further smeared 
by the PE counting statistics. The smeared Monte Carlo 
spectrum was then fit to a Gaussian (or, two Gaus- 
sians for the case of 207 Bi) to determine the Eq of the 
source conversion electron peaks. The measured QADC 
spectrum was also fit according to the same procedure, 
thereby providing the measured QADC peak location for 
each conversion electron peak. Note, however, in order to 
prevent PMT nonlinearity from shifting the fit position, 
the QADC spectrum was first linearized using a linear- 
ity function fi determined from an earlier iteration of 
this procedure. The linearized spectrum of /r (QADC i ) 
was then fit, and the peak positions were converted back 
to nonlinear QADC values via the fi functions. Figure 
13 shows sample data compared with Monte Carlo spec- 
tra of the visible energy E v - ls for three of the calibration 
sources ( 139 Ce, 113 Sn, and 207 Bi). 

Using the position response map rji(x, y) for each PMT, 
the expected light reaching each PMT rji(x,y) ■ Eq was 
plotted against the observed QADC values of the peaks 
for the various sources at the various (a;, y) positions. The 
points on these plots thus trace out the linearity curve 
for each PMT, with a fit to the rji(x, y) ■ Eq values as a 
function of the QADC values providing a parametriza- 
tion of the linearity function /r for each PMT. Ex- 
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FIG. 13. Calibrated E v [ s visible energy spectra (data points with error bars), compared with Monte Carlo calculations (solid 
lines), for three of the calibration sources: 139 Ce (left panel), 113 Sn (middle panel), and 207 Bi (right panel). 
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FIG. 14. Two examples of PMT linearity curves, for a PMT exhibiting a nearly linear response (left panel) and a PMT 
exhibiting a highly nonlinear response (right panel). The curves are fits to the measured QADC values (horizontal axis) of the 
peaks for the various calibration sources at different (x, y) positions versus the expected light reaching each PMT. The large 
nonlincarity observed in some PMTs was traced to damaged bases and subsequently repaired for runs after 2009. 



amples of such linearity curves are shown in Fig. 14 for 
two of the PMTs, with one of the PMTs exhibiting a 
nearly linear response, and the other a highly nonlin- 
ear response 5 . The linearity curves were parametrized 
with a fit function which was purely linear below some 
manually-determined transition point xq. Above this 
Xo, the nonlinearities were parametrized in the form 
y(x) = exp[Ci + In a; + C 2 [m(x/x )] 2 + C 3 [ln(x/xo)] 3 ], 
with C\, C2, and C3 fitted constants. 



a The source of these nonlinearities was later (in 2010) found to 
likely be the result of problems with the PMT bases. In par- 
ticular, it was found that several of the interstage capacitors 
were electrically shorted. After this problem was discovered, the 
PMTs and bases were replaced. 



3. PMT Position Response Maps 

The overall principle for the determination of the light 
transport efficiency to each PMT as a function of posi- 
tion, r)i(x,y), was to employ the measured neutron f3- 
decay spectrum endpoint as a "standard candle" pro- 
viding coverage over the entire detector fiducial volume. 
Assuming the linearity and time-dependent gain func- 
tions, fi and gi(t), are known, the measured data provide 
a measure of the amount of light reaching each PMT, 
Li(x,y) = TH(x,y) ■ E Q = fr^QADCJg^t)). The spec- 
trum of Li(x,y) has the same shape as the spectrum 
of Eq , but is linearly "stretched" by the light-transport 
factor r]i(x,y). Determining the absolute light-transport 
efficiency would be quite difficult, but, fortunately, is un- 
necessary. Instead, measuring the relative efficiency be- 
tween different (x, y) locations is sufficient, normalizing 
r\i by convention to 37^ (0, 0) = 1, and leaving the un- 
known constant factor between the rji and the absolute 
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(unknown) light-transport efficiencies to be absorbed into 
the /,-. 

In principle, the relative rj between two different lo- 
cations can be determined by seeing how much the Li 
spectrum at one point needs to be "stretched" to line 
up with the spectrum at another point. In practice, the 
stretching factor r]i(x,y) was determined so that a Kuric 
plot of r]i(x, y) ■ Li(x, y) was aligned with a Kurie plot of 
a smeared Monte Carlo visible energy spectrum. To find 
the necessary rji(x,y), an iterative Kurie plotting proce- 
dure was used. Starting from an initial guess r] i (x,y) 
for each rji(x,y), a Kurie plot was then made from the 
spectrum of r}ifl{x,y) ■ Li(x,y). This plot was fit with 
a straight line over a visible energy range from 250-700 
keV, which yielded some intercept E- iat . The position 
of this fitted intercept relative to the expected intercept 
(from the Monte Carlo visible energy spectrum smeared 
by the PE counting statisics), £"mCj then provided an im- 
proved estimate for the stretching factor, 77^ = ^ MC ??i,o. 
This procedure was iterated several times (including it- 
erations to the Monte Carlo spectrum, to account for 
changes to the PE counting statistics resulting from im- 
proved estimates for the light-transport efficiencies at 
each (x,y) position with each iteration), until the inter- 
cepts from the Kurie plots for all of the points over the 
detector fiducial volume were aligned with the smeared 
Monte Carlo spectrum intercept. Note that the Monte 
Carlo spectrum of the visible energy varied with the par- 
ticular Geometry (i.e., A, B, or C/D, depending on the 
decay trap end-cap foil and MWPC window thicknesses). 
Thus, even though the visible energy spectra varied with 
Geometry, we emphasize that the Kurie fits to the visible 
energy were employed for relative point-to-point rji(x,y) 
calibration purposes, and not for absolute energy-scale 
calibration purposes. 

The above-described procedure for construction of 
these PMT position response maps was implemented by 
combining /3-decay data from nearly the entire 2008- 
2009 datasct. The detector face was divided into 
180 approximately-equal-sized pixels, and a background- 
subtracted Li(x,y) spectrum for each PMT was gen- 
erated for each of these pixels. The stretching factor 
r]i(x,y) was then determined for each of these pixels. 
Having determined the rji(x,y) for each of these discrete 
pixels, a continuous map of the light-transport efficiency 
was then produced via bicubic spline interpolation in po- 
lar coordinates. Fig. 15 shows an example of an interpo- 
lated position response map for one PMT, exhibiting the 
expected strong dependence with position. 

The primary source of uncertainty in the construction 
of these PMT position response maps was the statisti- 
cal uncertainty in the Kuric fitting procedure, which was 
±1.0% at each of the 180 pixels and independent of the 
rji(x,y) value determined for each pixel. On an event- 
by-event basis, there was an additional uncertainty from 
the coupling between the varying (x, y) position response, 
and the uncertainty in the reconstruction of the position 
by the MWPC. This was determined to be ~ ±1.5% per 
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FIG. 15. Example of a r/(x,y) light-transport efficiency map 
for a single PMT, nominally viewing the [x < 0, y > 0) quad- 
rant of the scintillator disc. By convention, 77(0, 0) = 1. 

RMS mm uncertainty in position reconstruction. Note 
that the RMS uncertainty in the position reconstruction 
is difficult to estimate, because on an event-by-event ba- 
sis, there will be a true variation due to the Larmor spiral 
radius. This RMS uncertainty is actually probably sig- 
nificantly less than 2 mm, as the width observed with 
the conversion electron calibration sources was mostly 
explained by their Larmor radii. Because the position re- 
sponse maps between the four PMTs for each detector are 
significantly correlated, the uncertainty of the combined 
four PMT response will be less than the ~ 0.75%/mm 
expected under the assumption that the four PMT un- 
certainties arc independent. 

4- Gain Stabilization 

The gain of the PMT/QADC system drifted over time, 
typically on the scale of a few percent over several hours, 
due to ambient temperature fluctuations, etc. Periodic 
113 Sn source calibrations conducted every few days dur- 
ing production /3-decay running provided the first layer 
of gain stabilization. For such a run, the (x, y) posi- 
tion of the source and its QADC spectrum for each PMT 
was measured. The observed QADC spectrum and a 
Monte Carlo spectrum for that (x, y) position were then 
fit with the same procedure employed for the determi- 
nation of the PMT linearity functions, yielding a QADC 
channel value for the 113 Sn peak observed in each PMT 
i, QADC Sni , and a Monte Carlo expected light yield, 
EQ,Sn,i- The PMT gain factor, gi, was then set so that 
the calibration curve placed the measured QADC value 
at its expected Eq, 

QADC Sn j = 9i ■ fi(rn(x,y) ■ L SM ). (18) 
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Note that the resolution of each PMT was also deter- 
mined from the fitted width of the 113 Sn peak spectra. 
As with the energy calibration, the resolution was ex- 
tracted from a comparison to Monte Carlo spectra, in 
order to disentagle the effects of peak broadening due to 
PE counting statistics from the multiple K, L, etc. con- 
version electron lines. 

All other runs were then further gain stabilized to 
match these 113 Sn calibration runs by comparing, on a 
run-by-run basis, the shape of the measured QADC spec- 
trum for events tagged as background cosmic-ray muons 
by the backing veto. Typically ~ 5000 such events were 
identified in the ~ 1-hour long /3-decay runs, and the 
measured QADC spectrum was fitted to a Landau dis- 
tribution. If the QADC peak position for a muon event 
in PMT i was QADC jUi (0) during the 113 Sn source cal- 
ibration run and then QADC^^t) at some later time t, 
this gain shift was then corrected for by setting the time- 
dependent gain correction factor to be 



9i(t) 



Mmfav) • E Q ,sn,i) QADC^(O) 



QADC Sn , 



QADC^tt) 



(19) 



These time-dependent gain corrections were typically on 
the order of ~ 5%. 

The uncertainty in the 113 Sn gain stabilization was 
dominated by the position response uncertainty; since 
the source data were very localized in position, they were 
subject to the statistics-limited localized position map 
fluctuations of ~ ±1% for each PMT. The fit statistics 
for the cosmic-ray muon peaks were typically ~ ±1.5% 
for each PMT on a run-by-run basis. These two un- 
certainties must then be combined in quadrature, with 
that from the muon peak uncertainty contributing twice, 
since calculation of the gi(i) factors requires results from 
the muon peak fits for both the 113 Sn calibration run 
and the run being calibrated. Therefore, the total gain 
stabilization uncertainty for each PMT was estimated to 
be ~ 2.3%. With the four PMTs contributing approxi- 
mately equally to a combined reconstruction of the visi- 
ble energy (discussed next), the gain stabilization fluctu- 
ations were estimated to be ~ 1.2%. 



I. Visible Energy Reconstruction and Resolution 

Thus far, the discussion has primarily focused on the 
calibration of the individual PMTs. The calibrated sig- 
nals from the four PMTs for each scintillator then pro- 
vided, on an event-by-event basis, four independent mea- 
surements of the visible energy. These were then statis- 
tically combined, with weighting according to their PE 
counting statistics, to produce a single, more accurate 
result for the measured visible energy. Note that the 
dominant contribution to the uncertainty, Cj, in each in- 
dividual PMT's measurement of the visible energy was 
from Poisson counting statistics, as the previously dis- 
cussed position response map (~ 1.5%) gain stabilization 



(~ 2.3%) uncertainties are small compared to the order 
~ 10%- level individual-PMT PE counting statistics. 

The motivation for extracting the event energy from 
a statistically weighted average of the individual PMT 
energy measurements (as opposed to a non-statistically- 
weighted sum of the individual PMT responses) is as fol- 
lows. Consider an event with quenched energy Eq at 
position (x,y). According to our model some fraction 
of the produced scintillation light Lj = rji(x,y)EQ will 
reach each PMT, where it will be converted to Ni ± y/Ni 
PEs according to the quantum efficiency, Cj, of the PMT, 
such that Ni = C,L;. The QADC signal for PMT i is 
then converted to an individual-PMT estimate E i: with 
estimated error <?i = Ei/y/Nl. Combining the four PMT 
estimates, with their respective 1/of statistical weights, 
and assuming the individual PMT measurements are 
such that Ei w Eq, we find 
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which proves that the statistically weighted mean yields 
an estimate for the energy which is the product of the 
sum of the total number of photoclectrons, N to t = Ei -^ii 
and the position-dependent photoelectron-to-cnergy con- 
version factor, 1/ Ei QVifei u)- Note that this form sig- 
nificantly protects against errors in the reconstructed 
(x, y) position, because the position dependence of the 
individual PMT responses appears only in the summed 
combination Ei^ 7 ?^ which is a smoother function of 
(x, y) than the individual maps. 

The energy resolution of the detector was extracted 
from fits to the measured 113 Sn peak positions and 
widths, after accounting in Monte Carlo for peak broad- 
ening from the K, L, etc. conversion electron lines. Fig- 
ure 16 shows plots of the extracted number of photoclec- 
trons as a function of position. Averaged over the fidu- 
cial volume, the detector resolution was such that ~ 400 
PEs/MeV were observed in the East detector, and ~ 500 
PEs/MeV in the West detector, translating to a resolu- 
tion of ±9% at the 113 Sn peak and ±5% at the neutron 
/3-decay endpoint energy. 



J. Visible Energy Reconstruction Uncertainty 

The quality of the visible energy calibration procedure 
was checked by comparing how closely the reconstructed 
visible energy spectra for the conversion electron calibra- 
tion sources were aligned with the Monte Carlo predicted 
spectra. The primary findings of this exercise were: 
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FIG. 16. Extracted energy response for one of the detectors in photoelectrons per MeV as a function of (x, y) position extracted 
from one PMT (left panel) and the statistically weighted combined response from all four PMTs (right panel), demonstrating 
the combined response exhibits a smoother response as a function of the {x,y) position. 
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FIG. 17. (Color online) The vertical axes appearing at the locations of the expected (i.e., as calculated in Monte Carlo) visible 
energies (on the horizontal axis) for the various calibration sources show histograms of errors (in keV) in the measured visible 
energies. The dashed lines then show our conservative estimate of an "error envelope" of the total energy calibration uncertainty 
based on the widths of the error distributions. 



• For each run, there was an overall ±1.2% energy 
uncertainty from gain stabilization. 

• At low visible energies (50-100 keV), there was a 
±5% linearity uncertainty, as deduced from a com- 
parison of measured and simulated 109 Cd and 139 Ce 
spectra. 

• The linearity uncertainty was ~ at the 113 Sn peak 
energy, since this served as an anchor point for the 
calibrations. At this visible energy range, the vis- 
ible energy reconstruction error, ±1.7%, was pri- 
marily due to the errors in the position response 
maps and the run-to-run gain stabilization uncer- 
tainty. 



• Residual nonlinearity induces a 1.3% (~ 10 keV) 
error at the /3-decay endpoint energy, based on 
fits over the visible energy range of 300-700 keV. 
This is then interpreted as the linearity uncertainty 
around 500 keV. The linearity uncertainty then in- 
creases to ±2.5% at the upper end of the visible 
energy range, ~ 900 keV, based on the 207 Bi upper 
peak. 

Figure 17 shows a histogram of the reconstruction er- 
rors in the source peak energies from calibrations across 
all of the Geometries (A, B, C, and D), together with an 
"error envelope" of the total calibration uncertainty from 
the above described sources. Note that this error enve- 
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lope is a conservative estimate for the uncertainty, based 
on the widths of the distributions (and is especially con- 
servative at energies below 200 keV, where the envelope 
is wider than the plotted 109 Cd, 139 Ce, and 114m I n distri- 
butions), as opposed to the uncertainties in the means. 
The impact of this uncertainty in the visible energy cali- 
bration on the extraction of the /3-asymmetry is discussed 
later in Section VIII C. 



K. Scintillator Trigger Efficiency Functions 

An extraction of the scintillator two-fold PMT trig- 
ger efficiencies as a function of visible energy was impor- 
tant for comparisons of Monte Carlo calculations with 
data (especially for backscattering events which deposit 
small amounts of energy in the scintillator), because the 
trigger efficiency ultimately impacts the reconstruction 
of the event type. We extracted these two-fold PMT 
trigger efficiencies from the measured /3-decay data as 
a function of the measured visible energy according to 
the following procedure. First, we employed only elec- 
tron events satisfying the MWPC position cuts described 
in Section VF2. Second, for every possible pair of 
two PMTs that both triggered (i.e., generated a signal 
above pedestal), we then incremented an E v { s "total" 
histogram for the other two PMTs. For each of the 
other two PMTs, if it triggered, we then incremented 
an E vis "trigger" histogram. Third, the trigger effi- 
ciency histogram for each PMT, ej, was then obtained 
by dividing the "trigger" histogram by the "total" his- 
togram. Finally, the total two-fold PMT trigger efficiency 
for each scintillator as a function of the visible energy 
was calculated as [1 — (probability no PMTs trigger) — 
(probability only one PMT triggers)], with the appropri- 
ate combinatorics for each of the terms in this expression 
in terms of the individual PMT e; efficiencies. 

The results from such an analysis of the individual 
PMT efficiencies arc shown for one of the detectors in Fig. 
18 for one of the Geometries (Geometry C). The PMT 
efficiency curves for the other detector and the other Ge- 
ometries were similar. These were incorporated in our 
Monte Carlo simulation code. 



L. MWPC Energy Calibration 

In addition to suppressing gamma-ray backgrounds 
and permitting the definition of a fiducial volume, Monte 
Carlo studies suggested that the measured energy de- 
position in the MWPC on the triggering scintillator 
side could be used to separate Type 2/3 backscatter- 
ing events. For example, in the event types schematic 
shown in Fig. 7, the depicted Type 2/3 events would 
otherwise appear identical simply in terms of detec- 
tor signals above threshold. However, in the depicted 
Type 2 event, the electron traverses the MWPC on the 
triggcring-scintillator side only once, whereas in the de- 
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FIG. 18. (Color online) Individual PMT trigger efficiencies ti 
for the four PMTs (indicated by the different colors) on one 
of the detectors, as extracted from an analysis of Geometry 
C data. 



picted Type 3 event, the electron traverses the MWPC 
on the triggering-scintillator side twice. Therefore, for 
the Type 2/3 events depicted there, the energy depo- 
sition in the MWPC on the triggcring-scintillator side 
would be expected to be greater for the Type 3 event 
based on path length considerations. Hereafter, we will 
refer to the MWPC located on the earlier (and the only) 
triggering-scintillator side as the "Primary MWPC" , and 
the opposite-side MWPC as the "Secondary MWPC" . 

An energy calibration of the MWPC response was per- 
formed according to the following procedure. First, we 
performed polynomial fits to Monte Carlo data of energy 
deposition in the MWPC for Type neutron /3-decay 
events (i.e., the calibration was based on Type /3-decay 
events). Denote the resulting fit function /mc(-E')- Sec- 
ond, the MWPC detector face was divided into 10 mm 
x 10 mm 2 square bins (with 88 of these bins providing 
coverage of the decay trap circular geometry). In each 
of these 88 bins, we then fitted the function fMc(E(x)), 
where x denotes the MWPC anode PADC readout chan- 
nel number, and E(x) was a function (taken to be poly- 
nomials) that coverted from channel number to energy. 
Lookup tables in binned (x, y) positions and anode chan- 
nel numbers were then constructed for each detector for 
the different Geometries (i.e., A, B, C, and D). 

The quality of the MWPC energy calibration is shown 
in Figs. 19 and 20. First, in Fig. 19 we compare a baseline 
calibrated (and background-subtracted) Primary MWPC 
energy spectrum for Type neutron /3-decay events from 
Geometry B with Monte Carlo calculations for this Ge- 
ometry. Recall that the calibration was based on polyno- 
mial fits to the Monte Carlo Type (Primary MWPC) 
spectra. Second, in Fig. 20 we compare calibrated (and 
background-subtracted) Primary and Secondary MWPC 
energy spectra for Type 1 and Type 2/3 neutron /3- 
decay events from Geometry B with Monte Carlo cal- 
culations. Note that the Type 1 Primary and Secondary 
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FIG. 19. Baseline calibrated Primary MWPC energy spectra 
for Type neutron /?-decay events from Geometry B (solid 
line) compared with Monte Carlo calculations (dotted line). 
The histograms are normalized to unity. 

and the Type 2/3 Secondary spectra provide a pseudo- 
independent check of the calibration, as the Monte Carlo 
calculations were based on fits to the Type MWPC Pri- 
mary spectrum (i.e., the spectrum on the earlier, or only, 
triggcring-scintillator side). Note that the calibrated re- 
sponse exhibited a strong position dependence, which can 
be seen in Fig. 21, where we show position maps of the 
two MWPCs' anode channel numbers that correspond to 
a particular fixed energy (for this plot, 4.14 keV, the rel- 
evance of which for the separation of Type 2 /3 events is 
discussed later in Section VI A 2). 

M. MWPC Position-Dependent Efficiency 

As discussed earlier in Section VE2, a fixed MWPC 
PADC channel number cut was used to separate gamma 
ray and charged particle events (i.e., this cut did not vary 
with the reconstructed (x,y) position). If the MWPC re- 
sponse was independent of position, a plot of the MWPC 
efficiency for the identification of charged particles as 
a function of the measured MWPC energy deposition 
would exhibit a step function at the cut energy (or 
PADC channel number cut). However, as just shown 
in Fig. 21, the MWPC response was strongly position- 
dependent 6 ; therefore, for some given energy deposition 
in the MWPC, the identification of an event as either a 
gamma ray event or a charged particle event was subject 
to the MWPCs' position-dependent response. 

The extracted position-dependent efficiency for the 
identification of charged particle events is illustrated in 



We believe an electronics issue (capacitances and/or induc- 
tances) on the MWPCs' circuit boards may be the source of 
their observed position-dependent response. 



Fig. 22, where we have plotted the fraction of events 
passing the standard gamma ray PADC channel num- 
ber cut as a function of the calibrated energy deposition 
in the two MWPCs for two examples of typical data sets 
obtained during Geometries B and D. As can be seen 
there, over a particular range of energy depositions the 
fraction of events passing this cut varies monotonically 
from 0.0 to 1.0, as a result of the MWPCs' position- 
dependent response. These measured MWPC efficiency 
curves were incorporated in our Monte Carlo simulation 
codes. Note that difference between the two MWPCs' 
efficiency curves was smallest (greatest) for Geometry B 
(Geometry D), with the differences for Geometries A and 
C in between those of Geometries B and D. 

The possible impact of this position-dependent effi- 
ciency on the identification of gamma ray events (as op- 
posed to a cut on the MWPCs calibrated energy re- 
sponse) and thus on the measured asymmetry is dis- 
cussed later in Section VIII M. 



N. Initial Energy Reconstruction 

The initial energy of the electron, hereafter denoted 
-E'recon, was reconstructed from the measured visible en- 
ergy in the scintillator based on the results of GEANT4 
Monte Carlo simulations for the relation between the 
measured visible energy in the scintillator, E V - 1S , and the 
actual initial energy of the emitted /3-decay electron. 
Parametrizations relating -Brecon to E v - m were constructed 
for the different event types (Type 0, Type 1, and Type 
2/3) and for the different Geometries (A, B, and C/D). 
These parametrizations were based on GEANT4 simula- 
tions of conversion electron source spectra, and were ex- 
tracted from fits of the predicted mean value for E vis 
to the true initial source energy, for the various sources 
employed in the experiment. 

Monte Carlo generated source electrons were separated 
into Type 0, Type 1, and Type 2/3 events, according to 
the same selection rules as applied in the data analy- 
sis. For each Geometry and event type, two different fits 
were constructed: "Fit 1" was based on the scintillator 
visible energy, E v i s , only; whereas "Fit 2" included both 
the scintillator visible energy, -E V is, and the calibrated 
MWPC energy, E'mwpc- The input variables to the fits 
were the mean values of: (1) the effective true source 
energy, E^°^ cc , which accounted for the difference be- 
tween the Monte-Carlo-gencrated initial /3-decay energy 
and the electron's subsequent energy loss in the 3.6-/im 
thick source enclosure foils; (2) the scintillator visible 
energy, E v - ls , summed over both scintillators, which ac- 
counted for possible sub-trigger-threshold energy deposi- 
tion in one of the scintillators (e.g., in Type 2/3 events); 
and (3) the MWPC energy deposition, -Emwpc, summed 
over both MWPCs. For Type and Type 1 events, the 
functional forms for the fits of E rccon (= El°™ cc ) to E vis 
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FIG. 20. Calibrated Primary (top panels) and Secondary (bottom panels) MWPC energy spectra for Type 1 and Type 2/3 
neutron /3-decay events from Geometry B (solid lines) compared with Monte Carlo calculations (dotted lines). The histograms 
are normalized to unity. 
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FIG. 21. (x, y) position maps of the two MWPCs' anode PADC channel numbers that correspond to a given calibrated energy 
(in this case, 4.14 keV). 



and Emwpc 



E r , 



and for Type 2/3 events 



■Ere 



eE vis 
i2 



£ E vis 



f Q E vis + AE, (Fit 1) 

fqEvis + AE + -Emwpc : (Fit 2) 



(21) 



f Q E vis + AE, (Fit 1) 

f Q E vis + AE + Emwpc (Fit 2) 

(22) 

For each Geometry, values were fitted for /q, which rep- 
resented a numerical scintillator quenching factor; AE, 
the energy loss in the decay trap end-cap foils, MWPC 
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MWPC Energy [keV] MWPC Energy [keV] 

FIG. 22. (Color online) Efficiency for the identification of charged particle events by the MWPC detectors (for a fixed PADC 
channel number cut) plotted as a function of the measured MWPC energy deposition. The non-step-function behavior of these 
efficiency curves is a result of the MWPCs' position-dependent response. Typical results for the East (red circles) and West 
(blue triangles) detectors from Geometries B and D are shown in the left and right panels, respectively. 




u 50 100 150 200 250 300 350 400 450 500 
E recon [keV] 

FIG. 23. (Color online) Reconstructions of E leccm values ac- 
cording to the Fit 1 (solid) and Fit 2 (dashed) parametriza- 
tions for 113 Sn source calibration Type events from Monte 
Carlo (thin lines) and experimental data (thick lines). The 
histograms are normalized to unity. 



windows, etc.; and e, a second-order parameter which 
was used to control the fits for the Type 2/3 events. 
The resulting fit parameters are given in Table V, the 
point of which serves to set the scale for the Monte Carlo- 
calculated AE energy loss for the different Geometries. 

Figure 23 compares applications of these Fit 1 and 
Fit 2 parametrizations to -Brecon reconstructions of 113 Sn 
source calibration Monte Carlo and experimental data 
for Type events. Good agreement between Fits 1 and 2 
is observed. In the experimental data, the RMS width of 
Fit 2 was slightly smaller (by ~ 2 keV), suggesting that 
the inclusion of the MWPC energy slightly improved the 
resolution. However, we ultimately chose to use Fit I in- 
stead of Fit 2 for our asymmetry analysis because inclu- 
sion of the MWPC energy subjected the value of -Brecon 
to the possibility of an overflow value for the MWPC 
anode PADC readout (see Fig. 10), which occurred at a 



calibrated MWPC energy of ~ 10 keV. Because we could 
not reliably construct a value for the MWPC energy in 
the event of an overflow readout, there was a potential 
bias in the application of Fit 2 to the data. Nevertheless, 
the possible impact of this small difference in the Fit 1 
and Fit 2 results is discussed later in Section VIII C. 



VI. ASYMMETRY ANALYSIS 

In this section we outline our extraction of the asym- 
metries from the calibrated data. We begin by defining 
our different Analysis Choices for the extraction of the 
asymmetries under which we included or excluded the 
various backscattering event types. Next, we describe 
our procedure for the extraction of the binned (in en- 
ergy) background-subtracted event rates (and their sta- 
tistical uncertainties) from the /3-decay and background 
runs, and proceed to an extraction of the asymmetries 
under the octet data-taking sequences via the super-ratio 
technique discussed earlier. We then show results from 
a number of basic data quality assessment checks, in- 
cluding comparisons of the reconstructed energy spectra 
with Monte Carlo results, assessments of the stability 
of the energy scale with time (as quantified via Kurie 
fits to the /3-decay endpoints), and assessments of the 
statistical properties of the asymmetries under the octet 
data-taking sequence. 

A. Analysis Choices 

1. Definitions and Selection Rules 

We extracted the asymmetries according to a number 
of different Analysis Choices under which we included 
or excluded the various backscattering event types and 
varied the selection rules for the assignment of the elec- 
tron's initial direction. The motivation for this study 
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TABLE V. Values of fitted parameters for the E Iccon parametrizations. See text for descriptions of parameters. Note that the 
uncertainties in the fitted parameters were highly correlated; for brevity, we do not report the error matrix here. 


Event Type, Fit 


Fit Parameter 


Geometry A 


Geometry B 


Geometry C/D 


Tvne fl Fit 1 


AE [keV] 
fQ 


50.7(2) 
1.0461(4) 


63.0(2) 
1.0459(4) 


33.3(1) 
1.0476(2) 


Tvnp 1 Fit 1 


AE [keVJ 
fQ 


129(1) 
1.031(2) 


149(1) 
1.040(3) 


73.5(6) 
1.035(1) 


Type 2/3, Fit 1 


AE [keVJ 
fQ 

e [keV" 1 ] 


145(8) 
1.327(56) 
-2.33(65) x 10" 4 


152(7) 
1.437(48) 
-3.60(58) x 10" 


56(1) 
1.420(13) 
4 -3.26(19) x 10" 4 


Type 0, Fit 2 


AE [kcV] 
fQ 


45.1(2) 
1.0486(4) 


57.4(2) 
1.0484(4) 


27.3(1) 
1.0510(2) 


Type 1, Fit 2 


AE [keVJ 
fQ 


115(1) 
1.031(2) 


136(1) 
1.041(3) 


57.3(9) 
1.039(2) 


Type 2/3, Fit 2 


AE [keVJ 
fQ 

e [keV" 1 ] 


121(6) 
1.276(40) 
-1.88(47) x 10" 4 


129(7) 
1.396(53) 
-3.24(64) x 10" 


41(1) 
1.274(15) 
4 -1.92(22) x 10~ 4 



TABLE VI. Definitions of the different Analysis Choices for the inclusion or exclusion of the various backscattering event types, 
and the selection rules for the assignment of the electron's initial direction. "Trigger Side" refers to the scintillator generating 
the trigger. 

Assignment of Initial Direction for Backscattering Event Types 



Analysis Choice 


Type 




Type 1 






Type 2/3 


1 


Trigger 


Side 


Earlier Trig 


ger 


Side 


Trigger Side 


2 


Trigger 


Side 


Earlier Trig 


ger 


Side 


Omit 


3 (default) 


Trigger 


Side 


Earlier Trig 


ger 


Side 


If Primary -Bmwpc > 4.14 keV Trigger Side, Else Other Side 


4 


Trigger 


Side 


Omit 






Omit 


5 


Trigger 


Side 


Earlier Trig 


ger 


Side 


Apply Likelihood Function Prob(_BMWPc) for Assignment of Direction 


6 


Omit 




Earlier Trig 


ger 


Side 


Omit 


7 


Omit 




Omit 






Trigger Side 


8 


Omit 




Omit 






If Primary £mwpc > 4.14 keV Trigger Side, Else Other Side 


9 


Omit 




Omit 






Apply Likelihood Function Prob(_B]viwPc) for Assignment of Direction 



was to compare the variation of the measured asymme- 
try as a function of the analysis choice with the variation 
predicted by the Monte Carlo as a robust benchmark of 
our Monte-Carlo-calculatcd corrections for backscatter- 
ing and the cos ^-dependence of the acceptance. 

The selection rules for our various Analysis Choices, 
numbered 1-9, are summarized in Table VI. Note that 
Analysis Choices 1-5 consider Type events, with differ- 
ent selections for the inclusion/exclusion of Type 1 and 
2/3 events, and further selection rules for the identifi- 
cation of the initial direction of incidence for Type 2/3 
events. Analysis Choices 6-9 were included in order to 
study the asymmetry from backscattering events for di- 
agnostic purposes. 

We employed Analysis Choice 3 as our default method 
for the extraction of our final results for Aq, because this 
choice provided for maximal use of detector information 
and minimized the magnitude of the systematic correc- 
tions for backscattering and the (cos 9) dependence of the 
acceptance. Nevertheless, we show several results from 
the other Analysis Choices below, as these results provide 
a powerful validation of our Monte Carlo calculations. 



2. Treatment and Separation of Type 2 and Type 3 Events 

Note that in Analysis Choices 1 and 7 we assigned the 
initial direction of incidence for Type 2/3 events sim- 
ply to the triggcring-scintillator side. However, we at- 
tempted in Analysis Choices 3, 5, 8, and 9 to separate 
Type 2 and 3 events according to selection rules on the 
energy response of the Primary MWPC (i.e., recall, the 
MWPC on the triggering scintillator side). To illustrate, 
the top panel of Fig. 24 shows Monte Carlo calculations of 
the energy deposition in the Primary MWPC for Type 2 
and Type 3 neutron /3-decay events separately (of course, 
such a direct separation was not possible in data analy- 
sis) . The distinct spectra for these events suggested such 
a separation could be performed via a probabilistic (like- 
lihood) approach, and the feasilibity of such was studied 
in Monte Carlo. Results from this study are shown in 
the bottom panel of Fig. 24, where calculations of the 
probability that an event identified in data analysis as a 
Type 2/3 event was actually a Type 2 event are plotted 
as a function of £mwpc- 

The £mwpc > 4.14 keV energy cut employed in Anal- 
ysis Choices 3 and 8 was chosen such that if a fixed Pri- 
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FIG. 24. Top panel: Monte Carlo calculations of the en- 
ergy deposition in the Primary MWPC for Type 2 and Type 
3 neutron /3-decay events. The histograms are normalized 
to unity. Bottom panel: Monte Carlo calculations of the 
probability an event identified in data analysis as a Type 
2/3 event was a Type 2 event as a function of the energy 
deposition in the Primary MWPC. The dashed lines indi- 
cate the threshold cut of £mwpc > 4.14 keV for the assign- 
ment of the initial direction to the triggering-scintillator side 
in Analysis Choices 3 and 8 for the separation of Type 2/3 
events. The solid line is a fit to the calculations of the form 
Prob(i?MWPc) = -poerf(pi_E M wpc - Pi) + Pa, where erf de- 
notes the error function and the pi are fit parameters. This 
likelihood function was employed for the separation of Type 
2/3 events in Analysis Choices 5 and 9. 



mary MWPC energy cut was used to assign the initial 
direction of incidence for events identified as Type 2/3 in 
data analysis to either the triggering-scintillator side or 
the non-triggcring-scintillator side, the calculated proba- 
bility for the incorrect assignment of the initial direction 
to the triggering scintillator side was less than 50% (i.e., 
the calculated probability that the event was actually a 
Type 2 event). This threshold cut is consistent with the 
conceptual expectation (see, e.g., Fig. 7) that an observed 
Type 2/3 event which was actually a Type 2 event (i.e., 
initially incident on the non-triggering-scintillator side) 
would deposit less energy in the Primary MWPC (one 
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; k Type 3 
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FIG. 25. (Color online) Comparison between data (red tri- 
angles) and Monte Carlo (filled squares) for Primary MWPC 
energy spectra for Geometry A Type 2 and Type 3 events, af- 
ter separation of Type 2/3 events in data analysis according 
to the likelihood function Prob(_EMWPc)- 



traversal) as compared to a Type 3 event (two traver- 
sal). 

In Analysis Choices 5 and 9 we then actually separated 
the Type 2/3 events according to the calculated likeli- 
hood function Prob(£ , MWPc)i such as that already shown 
in the bottom panel of Fig. 24. (Note that although the 
curve shown there is from a Monte Carlo calculation for 
the Geometry A configuration, the curves for the other 
Geometries were similar.) Fig. 25 then demonstrates the 
good agreement observed between the data and Monte 
Carlo for the Primary MWPC energy spectra for Geom- 
etry A Type 2 and Type 3 events, separated according 
to the likelihood function. 7 



B. Event Rates and Statistical Uncertainties 

For any given Analysis Choice, events were binned into 
25 keV Brecon energy bins from 0-1200 keV, and then 
assigned an initial direction of incidence according to the 
selection rules in Table VI. The measured event rate in 
each of these -EVecon bins was then computed for each 
detector according to 



N h 



n>i 



T 



(23) 



where A^in denotes the number of events passing cuts in 
that particular Brecon bin, and T denotes that detector's 
live time (defined previously in Section VD). 



7 Note that this good agreement was actually to be expected, be- 
cause the calibrated Primary MWPC energy spectra for Type 
2/3 events agreed well with Monte Carlo (as shown previously 
in Fig. 20), and the Type 2 and Type 3 separation follows from 
the Monte Carlo calculation of the likelihood function. 
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In energy bins containing N < 20 counts (e.g., in 
bins beyond the /3-decay endpoint for /3-decay runs, or in 
nearly all bins for the shorter background runs), where 
the assumption of Gaussian errors approximated by Pois- 
son uncertainty is no longer valid, we assigned a statis- 
tical uncertainty to the rate in these bins according to 
the following procedure whose key assumption was that 
the underlying background and neutron /3-decay spectral 
shapes did not change with time, even if the run-by-run 
rates varied. 

First, we generated high-statistics parent /3-decay and 
background spectra for each of the detectors (for a partic- 
ular Geometry and spin state) by combining many runs. 
Second, on an individual run-by-run basis, we then com- 
puted the measured rate for each detector within an en- 
ergy window, r™^j ow , nominally 275-625 keV (see Sec- 
tion VII F), and compared this measured rate with the 
rate for the parent spectrum, ^^dow' m ^ ms same en_ 
ergy window. Third, the ratio of these two rates then 
defined a scaling factor / = J^Xw/ r £St for each 
of the detectors, which we then used to compute on a 
run-by-run basis the expected rate for each energy bin, 
r bin = r bfn Cnt x /> f° r eacn detector. Finally, the statis- 
tical uncertainty we then assigned to the measured rate 
in each bin was then 



rates given previously in Eq. (10), 



„cxp 



(24) 



where T again denotes that particular detector's 
(blinded) live time. For bins with > 20 counts, we em- 
ployed the usual vN uncertainties. 

For each /3-decay and background run pair, we then 
subtracted on a bin-by-bin basis the measured back- 
ground rate from the measured /3-decay rate, with stan- 
dard propagation of the statistical errors calculated ac- 
cording to the above-described procedure for N < 20 or 
N > 20 counts. We emphasize that this procedure af- 
fected only the assignment of the statistical errors; the 
actual measured counts were still employed for the cal- 
culation of the run-by-run rates for each /3-decay and 
background run pair. 



R = 



+ - 



i - Vr 

l + y/R" 



(25) 



2. Quartet Asymmetries 

Conceptually, the asymmetry for a quartet run pairing 
(i.e., either A-type A1-A12 or B-type B1-B12 runs in 
Table II) can be computed from a "summed super ratio" 
of detector rates, 



Rl 



-(A1,A2) 



-(A10.A12) 



[r +(A4,A5) +r +(A7,A9) ] 



+ (A4,A5) +(A7,A9) 1 



-(A1,A2) 



-(A10,A12)-i 



(26) 

where the notation is such that the run pairs within (• • • ) 
parentheses indicate the /3-decay run and ambient back- 
ground run background-subtraction pairing. The expres- 
sion for B-type runs is similar. The merit of this ap- 
proach is that that linear background and detector effi- 
ciency drifts cancel in this definition of R' . (This cancel- 
lation is exact provided that the durations of the /3-decay 
and background runs do not vary and that any such lin- 
ear drifts are constant in time.) 

In practice, to properly account for the statistical un- 
certainty, within each quartet we computed the statis- 
tically weighted mean of the detector rates for each of 
the polarization states. The resulting statistical uncer- 
tainties in the mean rates were then propagated into the 
super ratio R' and then into the calculation of the asym- 
metry. For example, we defined the detector rate for 
an A-type quartet to be the statistically weighted mean 
of the r +( A1,A2 ' an d r +( A10 > A12 ) ra tes. The quartet-based 
asymmetry was then extracted from the super ratio R' 
of these quartet-averaged rates. 

Note that a "product super ratio", R", can also be 
defined in terms of geometric means of detector rates for 
spin states as 



C. Asymmetry Extraction 

We then extracted the experimental asymmetry on an 
energy-bin-by-energy-bin basis for data grouped into ei- 
ther individual spin-state run pairs, quartets, or octets 
(all of which were defined previously in Section IV A 1). 



1. Spin-State Pair Asymmetries 

For individual spin-state run pairings (i.e., A1-A6, A 7- 
A12, B1-B6, or B7-B12 run pairings in Table II), we 
calculated the experimental measured asymmetries ac- 
cording to the definition for the super-ratio of detector 



R" = 



-(A1,A2) -(AlO.A^h r +(A4,A5) +(A7,A9) 1 



+ (A4,A5) +(A7,A9)i 



-(A1,A2) -(A10,A12)i 



1/2 



(27) 

We extracted asymmetries via both the summed and 
product super-ratio approaches, and the central values 
from the two methods differed by only 0.07%. How- 
ever, the problem with the product super-ratio method 
is that the resulting statistical uncertainty in R" is domi- 
nated by the rate from the run with the largest statistical 
uncertainty (e.g., if one of the /3-decay runs is signifi- 
cantly shorter in duration than the others). In contrast, 
the summed super-ratio method employs inverse-squarc- 
uncertainty weighting. Therefore, in our final analysis we 
employed the summed super-ratio method. 
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FIG. 26. (Color online) Comparisons of visible energy spectra for the different event types extracted from data (black squares) 
and Monte Carlo (red triangles) for the Geometry C configuration. The same selection rules for the event types were applied 
to both data and Monte Carlo. See text for details. All error bars shown are statistical; if not visible, the errors are smaller 
than the marker size. 



3. Octet Asymmetries 

Octet-based asymmetries were calculated similarly to 
quartet-based asymmetries. Now, for example, we de- 
fined the detector rate rf for a complete octet to be the 
statistically weighted mean of the r +( A1 ' A2 ). r +( A10 > A12 ) ^ 
r +(B4,B5\ and r +( B7,B9 ) rates. As with the quartet- 
based asymmetry, the octet-based asymmetry was then 
extracted from a summed super-ratio R' of these octet- 
averaged rates, with propagation of the statistical un- 
certainties in the octet-averaged rates through the super 
ratio and the asymmetry. Again, the merit of the octet 
approach is that linear background drifts cancel (subject 
to the same caveats as for the quartet asymmetries). 



D. Comparison of Data and Monte Carlo Visible 
Energy Spectra 

1, Visible Energy Spectra 

We now show an example (for the entire Geometry C 
data set) of comparisons of the measured background- 
subtracted visible energy E v - ls spectra with results from 



Monte Carlo (GEANT4 except where noted) calculations 
(assuming the Particle Data Group value for A$) m Fig. 
26 for the different event types. 8 The fractional residu- 
als, (Data — Monte Carlo)/Data, are shown in Fig. 27. 
Although not shown explicitly in this article, we note 
that we achieved the same level of agreement between 
data and Monte Carlo shown in Fig. VN for all of the 
Geometries and both detectors. 

Both the measured and Monte Carlo data shown there 
were extracted from events triggering one of the particu- 
lar scintillators. Type 2/3 events were separated in both 
data and Monte Carlo according to the selection rules 
of Analysis Choice 3. The Type 0, Type 2, and Type 3 
i?vis spectra shown there are the spectra observed in the 
triggering scintillator, whereas the Type 1 spectra are 
summed over both of the scintillators. 

The spectra shown there are for one of the neutron spin 
states (in particular, AFP spin- flipper off). Note that 



8 Note that comparing the measured E v i s spectra with Monte 
Carlo is equivalent to comparing the measured E 

recon SpCC- 

tra with Monte Carlo, because the mappings from -E v ; s values 
to -Erccon values were via the parametrizations extracted from 
Monte Carlo calculations, discussed previously in Section VN. 
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FIG. 27. Fractional residuals for the comparisons between the data and Monte Carlo visible energy spectra shown in Fig. 26. 
As can be seen in Fig. 26, there are very few Type 2 and Type 3 events with visible energies greater than ~ 500 keV. 



the total measured background-subtracted rate during 
runs when the neutron spin was flipped with the AFP 
spin-flipper was ~ 30% less than that during runs for 
the unflipped spin state, due to UCN losses along the 
transport guides between the AFP spin-flip region and 
the SCS decay-trap volume after the 2-T equivalent UCN 
energy increase following the spin flip. 



2. Backscattering Strengths and Fractions 



TABLE VII. Measured event type fractions for each Geometry 
integrated over the entire visible energy range from 0-800 keV 
and averaged over the detectors and the spin states. Type 2/3 
events were separated according to Analysis Choice 3. 



Geometry 


Type 


Type 1 


Type 2 


Type 3 


A 


0.947 


0.032 


0.011 


0.011 


B 


0.952 


0.030 


0.010 


0.008 


C 


0.930 


0.047 


0.009 


0.014 


D 


0.939 


0.039 


0.008 


0.014 



The normalization of the measured and Monte Carlo 
(GEANT4) spectra shown in Fig. 26 was performed accord- 
ing to their respective integral of the Type spectrum 
over the complete visible energy range of 0-800 keV. Af- 
ter this relative normalization, the GEANT4 spectra were 
further internally normalized to account for a known 
deficit in the GEANT4 backscattering strength [69, 70]. 
This was accomplished by applying two scale factors, 
/bulk (for backscattering from the scintillator bulk ma- 
terial) and /thin (decay trap and MWPC thin windows), 
to the GEANT4 backscattering distributions such that the 
simulated backscattering strengths matched the mea- 
sured backscattering strengths. Specifically, we applied 
/bulk to Type 1 and Type 3 events, whereas we ap- 
plied /thin to Type 2 events, because the former (latter) 



correspond to backscattering from the bulk scintillator 
(MWPC windows, gaseous materials, etc.). The numer- 
ical values of these scale factors were /bulk = 1-3 for all 
of the Geometries, and / t hi n = 1-3 (1.6) for Geometries 
A and B (C and D). The Type 1, Type 2, and Type 3 
Monte Carlo spectra shown in Fig. 26 are scaled by these 
scale factors. 

For completeness, the measured event type fractions 
integrated over the entire visible energy range from 800 
keV for the different Geometries are tabulated in Table 
VII. 
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FIG. 28. Example of a Kurie fit to the (background- 
subtracted) reconstructed energy E lecon spectrum for a ~ 1- 
hour long /3-decay run. The errors shown are statistical. 



E. Fitted Endpoint Distributions 

As a measure of the stability of the energy calibration 
with time, background-subtracted E Tecon spectra for each 
/3-decay run were converted to Kurie plots and fitted over 
the range of 400-700 keV. An example of such a fit to 
the -Ei-ocon spectrum from a ~ 1-hour long /3-decay run 
is shown in Fig. 28. There, we fitted the measured rate, 
binned in energy, dW/dE e , to the function 




cx(To-T e ), (28) 



where T denotes the fitted endpoint (kinetic) energy, 
p e = y/Tj? + 2T e m e is the magnitude of the electron 
momentum, T e (= E Tecon ) is the kinetic energy, and 
dE e /dp e = p e /E e = p e / (T e + m e ). 

Fits to Kurie plots of Monte Carlo spectra for the two 
spin states (in order to account for the /3 cos 6 factor 
in the angular distribution) which included the energy- 
dependent recoil order effects and the finite scintillator 
energy resolution yielded fitted values for the endpoint 
of ~ 787 keV (with a ~ ±0.4 keV difference for the two 
spin states). 

Distributions of fitted endpoints for the two detec- 
tors (employing "Fit 1" for the E lecon reconstruction), 
summed over all Geometries and spin states and weighted 
by the inverse square of their statistical uncertainties are 
shown in Fig. 29. For both detectors, typical run-to- 
run fluctuations (as characterized by the RMS) were less 
than ~ 13 keV. However, there was a systematic ~ 10- 
14 keV difference between the mean fitted endpoints and 
the Monte Carlo result of ~ 787 keV, which is addressed 
later in Section VIII C. 
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FIG. 29. Distributions of fitted endpoints for the two de- 
tectors, averaged over Geometries and spin states. The his- 
togram contents were weighted by the inverse square of their 
uncertainties. 



F. Statistical Properties of the Asymmetries 

As a demonstration of the utility of our octet data- 
taking procedure, we extracted values for the measured 
blinded asymmetries, j4 meas , under Analysis Choice 3 in- 
tegrated over a 225-675 keV energy window for the three 
different run groupings discussed previously: spin-state 
pairs, quartets, and octets. Incomplete octets were not 
discarded; instead, we retained individual spin-state pairs 
or quartets within incomplete octets in our analysis. 

Sample results from an analysis of the Geometry B 
data set are shown in Fig. 30. As can be seen there, 
the x 2 /^ value for the quartet analysis is improved over 
that for the spin-state pair analysis. Further, there is a 
slight shift in the central value of the asymmetry between 
the spin-state pair analysis and the quartet analysis. If 
there were linear drifts in the backgrounds and/or detec- 
tor efficiencies, this is the expected result, as any such 
linear drifts would bias the spin-state pair analysis, but 
would cancel in the quartet analysis. In comparing all 
of the Geometries, the x 2 /V values for the octet analy- 
ses of Geometries A, B, C, and D were 30.1/29, 12.8/22, 
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FIG. 30. Extracted values for the (blinded) measured asym- 
metries, ylmcas, from the Geometry B data set obtained un- 
der Analysis Choice 3 for spin-state run pairings (top panel), 
quartet run pairings (middle panel), and octet run pairings 
(bottom panel). 



20.1/13, and 6.5/7, respectively (i.e., the relatively small 
X 2 jv value for Geometry B was not representative of the 
entire data set). 

All of the asymmetry results presented hereafter were 
obtained under the octet analysis. 



VII. SYSTEMATIC CORRECTIONS FOR A 
EXTRACTION 

In this section we discuss the systematic corrections 
for backscattering, the cos ^-dependent acceptance (here- 
after also termed the "angle effect"), and the Standard 
Model recoil-order and radiative corrections that were 
applied to the measured asymmetries in order to ex- 
tract the desired /3-asymmetry parameter Aq. The Monte 
Carlo corrections for the backscattering and angle ef- 
fect corrections presented here were based on our GEANT4 
Monte Carlo simulation code. However, as noted ear- 
lier, we also developed a PENELOPE simulation code and, 



in general, we obtained good agreement between the 
GEANT4- and PENELOPE-calculated corrections. The small 
differences between the results from these two simulation 
programs are discussed in the context of the systematic 
uncertainty for the backscattering and cos ^-dependent 
acceptance corrections later in Section VIII L. 



A. Monte Carlo Benchmark: Scaling of 
Asymmetries with Analysis Choice 

First, benchmark of the validity of our Monte 

Carlo treatment of the angle and backscattering effects, 
we demonstrate that the scaling of our (still-blinded) 
measured asymmetries, j4 meas , with the Analysis Choice 
(calculated according to the selection rules in Table VI) 
is consistent with our Monte Carlo predictions. An ex- 
ample of this is shown in Fig. 31, where we have plotted 
the (blinded) measured asymmetries A meas and both the 
GEANT4 and PENELOPE Monte Carlo predictions for A meas , 
integrated over an energy window of 225-675 keV, as a 
function of the Analysis Choice for each of the Geome- 
tries. Note that the Monte Carlo predictions shown there 
were arbitrarily scaled (for effectively an arbitrary Aq) in 
order to match the data at Analysis Choice 3 (given that 
the measured asymmetries shown there were blinded); 
the point of this exercise was to demonstrate that the 
scaling of the Monte Carlo-calculated asymmetries with 
the Analysis Choice matched that of the measured asym- 
metries. [An absolute comparison would have required 
scaling each by the (unknown) unblindcd values of Aq 
for each Geometry] 

This agreement between the scaling of our measured 
and Monte-Carlo-calculated results for ^4 mea s with the 
Analysis Choice (and, thus, upon the inclusion or ex- 
clusion of the different backscattering event types, with 
different selection rules therein for identification of the 
backscattering event types) provides a powerful demon- 
stration of the validity of our Monte Carlo-calculated cor- 
rections for backscattering and the cos ^-dependence of 
the acceptance, and also validates our application of the 
scale factors /bulk and /thin to the GEANT4 backscattering 
distributions, discussed previously in Section VI D 2. 

Hereafter, all of the results shown for the asymmetries 
and the Monte Carlo-calculated systematic corrections to 
the asymmetries were obtained under the default Analy- 
sis Choice 3, the motivation for which was noted previ- 
ously in Section VI A 1 . 



B. Asymmetry Unblinding 

All of the asymmetry results shown hereafter consti- 
tute our final (unblinded) results. We note that during 
our actual analysis of the blinded data, we did not un- 
blind our asymmetries until after all of the (already dis- 
cussed) calibrations and cuts, the systematic corrections 
and uncertainties now being discussed, and the analysis 
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FIG. 31. Comparison of the scaling of the (blinded) measured asymmetries extracted from the data, ylmeas, and the Monte 
Carlo-calculated asymmetries (from both GEANT4 and PENELOPE) with the Analysis Choice for each of the Geometries. The 
asymmetries shown here were integrated over an energy window of 225-675 keV. The results from the Monte Carlo calculations 
were scaled for each Geometry in order to match the data at Analysis Choice 3. All errors are statistical, with the GEANT4 
statistical errors smaller than the symbol size. 



energy window (discussed later in Section VII F) were fi- 
nalized. At that time, the asymmetries were unblmdcd 
by removing the scale factors applied to the detector rates 
which were used to blind the data (as discussed previ- 
ously in Section VC). 

We emphasize that our final results for Aq presented 
later were those obtained at the time the data were un- 
blinded; no further data analysis was conducted after the 
unblinding. 



C. Results for Measured Asymmetries 

Our resulting measured asymmetries, A meas (E lecon ) , 
are shown in Fig. 32 as a function of E recon in 25 keV bins 
over the entire detectable energy range, 50-800 keV, for 
each of the Geometries. There, we also show for each Ge- 
ometry the measured background spectrum and the re- 
sulting bin-by-bin background-subtracted /3-decay spec- 
trum. 



D. Overview of Method for Systematic Corrections 
to Extract Ao 

For a given energy bin in "true initial energy" , the true 
energy-dependent physics asymmetry under the Stan- 
dard Model would be P„Ao(/3trueCOs6')(l + /ro), where 
P n is the neutron polarization, (/?true cos 9) denotes the 
average value of the product of the true electron velocity 
in units of c with the true electron pitch angle cos 9 for 
that energy bin, and /ro denotes the (energy-dependent) 
recoil-order physics correction to Aq. For now, we have 
omitted the small O(0.1%) radiative correction to the 
asymmetry discussed in the Introduction to this article. 

Therefore, it is clear that an extraction of the desired 
/3-asymmctry parameter Aq from the (energy-dependent) 
measured asymmetry A meas (£' recon ) in each i^econ bin re- 
quires a correct reconstruction of the mean (initial) en- 
ergy in each energy bin, a correction to the measured 
asymmetry for missed backscattering, and a calculation 
of the mean value of (/3 r ocon cos 9) in that energy bin. For 
now, we will assume that in each particular energy bin 

(Brecon) = (Etrue), Or (Brecon) = (Aruc); later, we will 

explore the implication for a systematic uncertainty to 
the asymmetry resulting from a possible -E^ccon 7^ E trno 
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FIG. 32. (Color online) Results for the measured background E lacon energy spectra (open circles) and the background-subtracted 
neutron /3-decay E lacon energy spectra (filled circles) for each of the Geometries. The spectra shown here were summed over 
both detectors and then averaged over the two neutron spin states. The resulting measured asymmetries A mcas (filled red 
squares; see vertical scale on the right) are shown as a function of -E rC con for each of the Geometries. These are compared with 
the Monte Carlo calculated values for the measured asymmetries (solid red lines). The energy dependence of the measured 
asymmetries is primarily due to the /? cos dependence of the measured asymmetries [see Eq. (8)] ; there are also (/? cos 9) 
acceptance and backscattering effects (e.g., compare Geometry B). All errors are statistical. 



error, such as from an energy calibration error. Note 
that even under the assumption that E rccon = E tTUC (i.e., 
"perfect" calibration), an acceptance correction must still 
be applied for the average value of (/3 cos 9) in each en- 
ergy bin, because the electron energy loss (in the decay 
trap end-cap foils, MWPC windows, etc.) is strongly 
angle dependent (hence the designation "angle effect"), 
with the energy loss a monotonically increasing function 
of the pitch angle. Therefore, each E Tecon bin includes 
a distribution of events in initial true energy and pitch 
angle. 



We used our Monte Carlo simulation code to compute 
these required corrections for missed backscattering and 
the ((3 cos 9) angle effect. After extraction of the mea- 
sured asymmetry A meas (-Brecon) in each E TCCOn energy 
bin, we then applied the following corrections in a sequen- 
tial manner in order to extract Aq from A meajs (E Tecon ): 



1. Asymmetry Energy Dependence and /3-Decay Spectra 

(1) First, we made a first-order bin- by-bin correction 
for the energy-dependence of the measured asymmetries, 
Ancas (-Brecon), by calculating an asymmetry Ax(E Tecoa ), 
defined by 

A\ (Brecon) — ^mcas/ ( ^/^recon) 5 (^^) 

where /3 r econ was calculated for the central value of each 
Brecon bin for each event class. Note that this first-order 
correction simply assumed a uniform value (or distribu- 
tion) of /3 r ccon within each energy bin, and also a sym- 
metric value for (cos 9) of 1/2 in each energy bin. 

2. Backscattering Corrections 

(2) Second, we applied a bin-by-bin correction for 
missed (or misidentified) backscattering events, with the 
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result of this correction ^(-Brecon) denned to be 

ccon ccon 

)(1 + A 2 ), (30) 

where A 2 = A2,o + A 2 ,i + A2 ; 2 represents the total 
backscattcring correction in each E Tecon bin, which we de- 
fine to be the sum of the individual corrections for events 
misidentified as Type 0, Type 1, and Type 2/3 events, 
respectively. The fractional correction factors A2,j were 
extracted from Monte Carlo. 



3. (j3 cos 9) Acceptance Corrections 

(3) Third, for each event class we applied a bin-by-bin 
correction for the (/3cos#) angle effect, with the result of 
this correction A3(E Iccon ) defined to be 

^(-Brecon) = ^4-2 (-Brecon) X -/3 r ccon/ (ftruc COS 9) 

= A 2 (E Iccon )(l + A 3 ). (31) 

Again, the fractional correction factors A3 were extracted 
from Monte Carlo. 

4- Standard Model Corrections 

(4) After application of the above three corrections, the 
resulting asymmetry ^i3(-B rccon ) is proportional to the 
product of P n , Aq, and the (energy-dependent) recoil- 
order and radiative corrections (which are calculable un- 
der the Standard Model). As a final step, we extracted 
a value for A in each E lecon bin according to 

ccon )/Pn = 

A 3 (E 

nX'Oll )[l + A RO (£ r ccon 

)][1 + A rad (-Brecon)], 

(32) 

where Aro (-Brecon) and A ra d (Brecon) denote the energy- 
dependent recoil-order and radiative corrections, respec- 
tively, for that particular -E rocon bin. In the absence of 
new physics, there should be no residual energy depen- 
dence to the values of the asymmetries Ao(E Tecon ) ex- 
tracted from all of the -B recon bins after application of 
these Standard Model corrections; therefore, the final 
energy-averaged value for A will be the statistically- 
weighted average of the bin- by-bin j4o(-E rccon ) values, 

ccon 

)). (33) 

We now describe the results from our Monte Carlo cal- 
culations of the A 2l , and A3 correction factors in more 
detail. 



E. Results for Monte Carlo Corrections 

Monte Carlo calculations of the A 2 ,i and A3 correction 
factors were carried out for each of the Geometries, thus 



accounting for their different foil thicknesses, their differ- 
ent measured MWPC efficiencies (recall the discussion of 
such in Section VM), and their different two-fold PMT 
trigger threshold functions (Section VK). To account for 
the previously-discussed (Section VI D 2) known deficit 
in the GEANT4 backscattering strengths, backscattering 
events in the Monte Carlo were re-weighted, with a 
scaling factor of /bulk applied to Type 1 and Type 3 
events, and /thin applied to Type 2 events and those 
events misidentified as Type events (i.e., the "missed" 
backscattering events discussed in Section VA). How- 
ever, note that the A3 correction factors were evaluated 
in Monte Carlo without the /bulk and /thin re-weighting 
factors for the backscattering events. After application 
of the A 2j j backscattcring corrections, the remaining A3 
correction factor accounts for the distortion of the de- 
tected (Brecon cos 6) acceptance from the "true" (^/3true) 
distribution. Therefore, employing these /bulk and /thin 
re- weighting factors for the evaluation of the A3 correc- 
tion factor in Monte Carlo would have biased the angular 
distribution of the simulated events, as the probablity for 
backscattcring increases with angle. 

The results of these calculations are plotted in Figs. 
33 and 34 as a function of E Tecon for each of the Geome- 
tries. Figure 33 shows the calculated results for the A 2 
and A3 corrections, and the size of their combined cor- 
rection, (1 + A 2 )(1 + A 3 ). For completeness, Fig. 34 then 
shows the sources of the individual A 2 ^ contributions to 
A 2 . All of the corrections shown here were calculated 
for the default Analysis Choice 3. The systematic uncer- 
tainties in these calculations are discussed in detail later 
in Section VIII L. Note that the results for A2 and A3 
are shown over an E recon energy range of 150-750 keV. 
Below 150 keV, the corrections become quite large, as 
the acceptance is highly suppressed for large pitch angle 
events and the probability for backscattering increases 
with decreasing energy. 

Our sign convention for these corrections is such that 
if the correction factor Aj > 0, the resulting asymmetry 
calculated according to Ai(l + Aj) is larger in magnitude 
(i.e., a more negative value for the asymmetry). As can 
be seen in Fig. 33, the A2 correction factors are > for 
all energies and Geometries, whereas the A3 correction 
factors are < 0. The conceptual physical explanation for 
this is as follows. 

With regard to a physical explanation for A2, the A 2 ^ 
correction factors correct the measured asymmetries for 
events misidentified in data analysis as a Type i event. 
For example, an electron incident initially on one of the 
detectors could backscatter from either the decay-trap 
end-cap foil or the MWPC entrance window and then be 
detected in the opposite side's MWPC and scintillator 
(e.g., the "Missed" event illustrated in Fig. 7); however, 
such an event would subsequently be misidentified in data 
analysis as a Type event (with an incorrect initial direc- 
tion). Because this type of event would dilute the magni- 
tude of the measured asymmetry for Type events, the 
Monte-Carlo-calculatcd A2.0 correction factor would be 
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FIG. 33. Results from GEANT4 Monte Carlo calculations of the energy dependence of the relative sizes of the A2 backscattering 
(thin solid line) and A3 angle effect (thick dashed line) systematic corrections. The energy dependence of the relative size of the 
combined correction (1 + A2)(l + A3) is then shown as the thick solid line. The systematic uncertainties in these calculations 
are discussed in detail later in Section VIII L. 



> (i.e., so as to increase the magnitude of the asymme- 
try, per our sign convention) to compensate for the dilu- 
tion. Further, as the probability for backscattering from 
plastic scintillator decreases with energy over the energy 
range of neutron /3-decay [70] (and also for backscatter- 
ing from thin Mylar films, as calculated in GEANT4 and 
verified via analytic integration of the differential cross 
section for the Mott scattering of electrons from atomic 
electrons), the A2,o correction factor decreases with en- 
ergy. Thus, conceptually, the /S.2,1 correction factors are 
expected to be > for all energies and to decrease in 
magnitude with energy, which is consistent with the re- 
sults from our Monte Carlo calculations shown in Fig. 
33. 

With regard to a physical explanation for A3 , there are 
two potential sources of bias to the (/3 cos 9) acceptance. 
First, there is a bias to the (cos 9) value, as the acceptance 
for large pitch angles is suppressed (e.g., from "Lost" 
events such as those shown in Fig. 7), with the value of 
(cos 9) > 1/2 for detected events (as the acceptance is 
biased towards small pitch angle events). Thus, in gen- 



eral, A3 is expected conceptually to be negative, given 
our definition of (1 + A3) = (^f3 rcC on) / (Pu-uc cos 9) in Eq. 
(31). Second, because the energy loss in the decay-trap 
end-cap foils and the MWPC windows is strongly angle- 
dependent, the E v i s visible energies in each visible energy 
bin, which then map to E rccon values via the (mostly 
linear) parametrizations discussed in Section VN, actu- 
ally correspond to a distribution of initial true energies. 
Therefore, there is a bias to the assumption that the value 
of (/Jtruc) = (Brecon) in each E recon bin, which is corrected 
for in our Monte Carlo calculations of A3 . The variation 
of A 3 with energy depends on the details and shape of 
the initial energy distribution. 

The differences between the values of the A2 and A3 
correction factors for Geometry C and Geometry D (de- 
spite their identical foil thicknesses) is primarily because 
the measured MWPC thresholds for Geometry D were 
higher than those for Geometry C. 

Finally, it is interesting to note that because A2 and 
A3 are of opposite signs, there is a "zero crossing" in the 
total correction factor (1 + A 2 )(l + A 3 ) sa 1 + A 2 + A 3 , 
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FIG. 34. Results from GEANT4 Monte Carlo calculations of the relative sizes of the individual A2,o (thin solid line), A2,i (thick 
solid line), and A2,2 (thick dashed line) contributions to the total A2 backscattering systematic correction. The systematic 
uncertainties in these calculations are discussed in detail later in Section VIII L. [The A2.1 backscattering correction is non-zero 
because some small fraction of events will initially backscatter from the decay-trap end-cap foils, and then undergo Type 1 
backscattering (however, the reconstructed initial direction will then be incorrect). The probability for this type of event 
increases with decreasing energy, and thus the acceptance for this type of events was suppressed in Geometries A and B (i.e., 
thicker MWPC windows).] 
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FIG. 35. Energy dependence of the Standard Model recoil- 
order correction (calculated according to the formalism of [21, 
22]) for an extraction of Aq from the measured asymmetry A3. 



which can be seen in Fig. 33. 

As discussed in the upcoming Section VIII L, the sys- 
tematic uncertainties in these corrections were taken to 
be a relative fraction of the magnitude of these correction 
factors. 



F. Analysis Energy Window 

Our final results for Aq presented later in Section IX 
were obtained over an analysis energy window of 275-625 
keV. This optimized energy window minimized the total 
absolute error in Aq resulting from the quadrature sum of 
the statistical error and the energy- dependent systematic 
errors in the above-described Monte Carlo calculations of 
the A2 and A3 backscattering and angle effects correc- 
tions (which varied with energy as just shown in Figs. 33 
and 34). 
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FIG. 36. Results for the asymmetries Ai (open circles), A2 
(open squares), A3 (filled circles), and Ao (filled squares) for 
each of the Geometries, integrated over the analysis energy 
window of 275-625 keV. All errors are statistical. 



G. Recoil-Order and Radiative Corrections 

Recoil-order corrections to the asymmetry were calcu- 
lated within the context of the Standard Model according 
to the formalism of [21, 22]. The numerical results from 
these two parametrizations agree to better than 2 x 1CP 5 . 
These parametrizations for the asymmetry A, Eq. (6), 
were then folded over the /3-decay energy spectrum (in- 
cluding contributions from the Fermi function). In gen- 
eral, these recoil-order corrections increase the magni- 
tude of the measured asymmetry over that of Ao; thus, 
an extraction of Ao from the measured asymmetry re- 
quires (per our sign convention) a negative correction 
(i.e., a decrease in the magnitude of the asymmetry). 
The energy-dependence of this correction is shown in Fig. 
35. For our analysis energy window of 275-625 keV, the 
integrated recoil-order correction (assuming the Particle 
Data Group average value for A) was (—1.79 ± 0.03)% 
to Ao, where the ±0.03% uncertainty corresponds to the 
Particle Data Group's statistical uncertainty in A. 

The value for the radiative correction was taken from 
the calculations of [25] , who presented results for the ab- 
solute (as opposed to relative) value of the radiative cor- 
rection to the asymmetry. These calculations were per- 
formed using the average value for A available at that 
time, and were reported at six discrete values of electron 
energy. In general, these radiative corrections increase 
the magnitude of the measured asymmetry by ~ 0.1% 
over the range of our 275-625 keV analysis energy win- 
dow. However, at the time of this analysis, we were not 
aware of the functional form for the energy-dependence 
of these radiative corrections which is presented in [24] . 
Therefore, we applied a (0.10 ± 0.05)% correction to our 
asymmetries (again, per our sign convention), where our 
estimated ±0.05% uncertainity accounts for differences 
between the value of A available at the time of the calcu- 
lation (1992, [25]) and its present value, and for our in- 



complete knowledge at the time of the energy-dependence 
for the radiative correction. 

Per the discussion in [15], application of this radiative 
correction to our reported value for Ao then permits ex- 
traction of a value for which can be compared with 
the expression relating Gf, t„, <m, and V u d, Eq. (3), 
in which the (1 + RC) electroweak radiative corrections 
have been factorized in the same way for both the vector 
and axial- vector interactions [15]. 



H. Comparison of Geometry-by-Geometry Results 

The values for the A2 and A3 corrections in Geome- 
try over our analysis energy window of 275-625 keV are 
tabulated in Table VIII for all of the Geometries. Figure 
36 then compares the resulting Geometry-by-Geomctry 
asymmetries A\ , A2 , A3 , and Ao integrated over the anal- 
ysis energy window. The final values for Ao from each 
of the Geometries are all seen to be statistically consis- 
tent. Geometry-by-Geometry results for the Ao values 
extracted bin-by-bin are shown in Fig. 37. 

We emphasize that the agreement between our re- 
sults from the different Geometries provides confidence 
in our Monte Carlo treatment of the backscattering and 
(/3 cos 6) acceptance effects. In particular, the result from 
Geometry B (thickest decay trap end-cap and MWPC 
foils) agrees well with those from the other Geometries, 
in spite of its relatively larger A2 and A3 corrections 
(compare the differences in Ai, A2, and A3 for Geome- 
try B with those for the other Geometries). 

Finally, for completeness, our background-subtracted 
neutron /3-decay rates integrated over the complete en- 
ergy range of 0-800 keV after all analysis cuts are listed 
in Table IX for all the Geometries and for both spin 
states. As previously discussed in Section VD, the data 
quality analysis cuts related to DAQ electronics prob- 
lems removed a significant fraction of the Geometry A 
and Geometry B events (up to ~ 30%). Another sig- 
nificant analysis cut included the 45 mm radius fiducial 
cut discussed in Section VF2, which removed ~ 25% of 
the events. Finally, under the 275-625 keV analysis en- 
ergy window, the rates from 0-800 keV listed in Table 
IX would be reduced by another factor of ~ 40%. 



VIII. SYSTEMATIC UNCERTAINTIES 

A. Summary 

Our systematic corrections and uncertainties are sum- 
marized in Table X, where we have categorized the ef- 
fects as either Geometry-Dependent (i.e., effects which 
varied with the decay trap end-cap foil and MWPC 
window thicknesses, measured detector thresholds, etc.), 
or Geometry-Independent (e.g., UCN polarization, dead 
time effects, etc.) In the rest of this Section we discuss 
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TABLE VIII. Values of the A2 and A3 correction factors over the analysis energy window of 275-625 keV for each of the 
Geometries. 



Geometry A Geometry B Geometry C Geometry D 



'rccon [kfiV] 


A 2 [%] 


A 3 [%] 


A 2 [%} 


A 3 [%] 


A 2 [%] 


A 3 [%] 


A 2 [%} 


A 3 [%] 


275-300 


2.28 


-1.17 


7.12 


-1.28 


1.57 


-0.43 


1.55 


-0.25 


300-325 


1.89 


-1.58 


6.32 


-1.52 


1.36 


-0.49 


1.41 


-0.24 


325-350 


1.63 


-1.38 


5.60 


-1.71 


1.30 


-0.40 


1.32 


-0.24 


350-375 


1.54 


-1.43 


5.12 


-1.94 


1.18 


-0.28 


1.19 


-0.25 


375-400 


1.33 


-1.45 


4.66 


-2.38 


1.04 


-0.50 


1.08 


-0.25 


400-425 


1.25 


-1.67 


4.19 


-2.61 


1.05 


-0.43 


1.00 


-0.28 


425-450 


1.16 


-1.66 


3.82 


-3.01 


0.91 


-0.53 


0.94 


-0.29 


450-475 


1.04 


-1.96 


3.51 


-3.56 


0.89 


-0.76 


0.90 


-0.29 


475-500 


1.04 


-1.91 


3.21 


-4.02 


0.87 


-0.63 


0.88 


-0.35 


500-525 


0.97 


-2.16 


2.96 


-4.39 


0.85 


-0.63 


0.80 


-0.42 


525-550 


0.92 


-2.51 


2.70 


-4.99 


0.77 


-1.00 


0.79 


-0.52 


550-575 


0.85 


-2.44 


2.46 


-5.57 


0.83 


-1.10 


0.83 


-0.65 


575-600 


0.84 


-3.03 


2.25 


-6.42 


0.73 


-1.20 


0.76 


-0.75 


600-625 


0.84 


-3.02 


1.99 


-6.86 


0.70 


-1.13 


0.80 


-1.13 



TABLE IX. Background-subtracted neutron /3-decay rates in- 
tegrated over the complete energy range of 0-800 keV after 
all analysis cuts. 







East Detector 


West Detector 


Geometry 


Spin State 


Rate [s _1 ] 


Rate [s _1 ] 


A 




4.79 


4.90 


A 


+ 


3.82 


3.25 


B 




6.50 


6.87 


B 


+ 


5.43 


4.78 


C 




3.29 


3.42 


C 


+ 


2.81 


2.44 


D 




6.45 


6.61 


D 


+ 


4.80 


4.12 



each of these systematic effects (in the order in which 
they appear in Table X) in more detail. 



B. Geometry-Independent: Dead Time 

Nearly all dead time effects cancel in the super-ratio 
technique. Indeed, in order for there to be any bias to 
the asymmetry resulting from dead time effects in the 
background-subtracted /3-decay rates, there must be a 
difference in the two detectors' dead times, and there 
must be a difference in a particular detector's dead time 
for the two neutron spin states. Thus, these effects arc 
expected to be quite small. 

Nevertheless, as previously noted in Section III H, the 
dead time of the DAQ system was monitored by counting, 
in scalers, the total number of detector two-fold PMT 
triggers, including those that were vetoed by the DAQ 
"busy logic" during the ~ 12 fis readout gates for the 
PADC modules. However, to avoid spurious (and cor- 
related) trigger chains from scintillator afterpulses (as 
noted earlier in Section III H) distorting the determina- 
tion of the dead time, the dead time was determined only 



TABLE X. Summary of systematic corrections and uncer- 
tainties. All numbers quoted are fractional [%] relative to Ao- 
Upper Table: Geometry-Independent systematic uncertain- 
ties. No systematic corrections were applied for these effects 
(with the exception of the radiative corrections, already dis- 
cussed in Section VII G). Lower Table: Geometry-Dependent 
systematic corrections and uncertainties. The quoted value 
denotes the systematic correction, with the error the system- 
atic uncertainty. As discussed in the text, A 2 represents the 
correction for backscattering, and A3 the correction for the 
angle effect, emwpc denotes the systematic uncertainty asso- 
ciated with the MWPC efficiency. 



Geometry-Independent Effect 



Uncertainty 



Dead Time 

Energy Reconstruction 

Fiducial Cut and Coordinate Systems 

Gain Fluctuations 

Live Time 

Magnetic Field Nonuniformity 
Muon Veto Efficiency 
Neutron-Generated Backgrounds 
Polarization 
Radiative Corrections 
Rate-Dependent Gain Shifts 



±0.01 
±0.47 
±0.24 
±0.20 

±0.24 
+0.20 
-0.00 
±0.30 
±0.02 

+0.52 
-0.00 

±0.05 
±0.08 



Geometry-Dependent Effects 



A 



B 



C 



D 



A 2 1.34 ±0.40 
A 3 -1.81 ±0.45 
emwpc 0.00 ± 0.02 



4.32 ± 1.30 
-3.22 ±0.81 
0.00 ±0.01 



1.07 ±0.32 
-0.60 ±0.15 
0.00 ±0.16 



1.08 ±0.32 
-0.36 ±0.09 
0.00 ±0.50 



from the scaler counts of detector triggers that occurred 
during triggers from the opposite-side detector or from 
other experimental triggers, such as the UCN monitors. 
The dead time, as extracted from the correlation between 
the DAQ trigger rate and the fraction of these "missed 
triggers" , was found to be ~ 13.5 /xs, which is consistent 
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FIG. 37. (Color online) Geometry-by-Geometry results for Ao values extracted bin-by-bin. The solid red lines indicate the 
analysis energy window of 275-625 keV. The quoted errors are statistical. 



with the nominal ~ 12 fis system dead time (associated 
with the gate for the PADC readout). Further, the dif- 
ference in the fraction of "missed triggers" for (up to) a 
20 s _1 trigger rate difference between the two spin states 
is no larger than ~ 0.03%. Considered together, any pos- 
sible bias to the asymmetry was no greater than 0.01%, 
which is the error we quote in Table X. 

Alternatively, another possible way dead time effects 
could bias the asymmetry is in the background subtrac- 
tion procedure, resulting from differences in the DAQ 
total trigger rates during /3-decay and background runs. 
However, these effects tend to cancel in the super-ratio, 
as the four background-subtracted /3-decay rates appear- 
ing in the super-ratio would be expected to be biased 
in the same direction. Further, the effect is minimized 
as the signal-to-background ratio increases. Under the 
conservative assumption of a 200 s _1 DAQ trigger rate 
difference (e.g., from differences in the scintillator trigger 
rates, UCN monitor trigger rates, etc.) for /?-decay ver- 
sus background runs, and a signal-to-background ratio 
greater than 5, any such systematic bias to the asymme- 
try from dead time effects is <C 0.01%. 



C. Geometry-Independent: Energy Reconstruction 

Figure 17 showed the error envelope for the uncertainty 
in the visible energy calibration. To estimate the system- 
atic error associated with possible errors in our energy 
calibration, we generated a large number (200 per Ge- 
ometry) of random error curves that were constrained to 
fit within the limits of this error envelope. We then ex- 
tracted from these error curves their contributions to an 
error in the asymmetry, resulting from an incorrect re- 
construction of the electron energy, and hence ((3 cos 6) . 
From these calculations we concluded that the maximum 
(i.e., worst case) error, resulting from the case where the 
error curves for the two detectors are identical, is a frac- 
tional 0.47% uncertainty in the asymmetry for the anal- 
ysis energy window of 275-625 keV. As a conservative 
estimate of the systematic uncertainty associated with 
our energy calibration, we then assign this worst-case er- 
ror of 0.47% to be the systematic uncertainty associated 
with possible errors in our energy calibration. 

As discussed earlier in Section V N, the default E lecon 
paramctrization we employed was based on a fit to the 
scintillator visible energy E v i s only; by constrast, an al- 
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tcrnative E recon fit included both E v { s and the calibrated 
MWPC energy -Emwpo To study the sensitivity of the 
reconstructed asymmetry to these two different -Brecon 
parametrizations, we extracted values for the energy- 
corrected asymmetry A\ for these two different fits. The 
difference between these two methods, averaged over the 
entire data set, was 0.2%. This is small relative to the 
0.47% systematic uncertainty associated with the energy 
calibration, and we also noted in Section VN that this 
alternative E Iccon paramctrization based on both E v - ls 
and £mwpc is subject to (uncorrectable) overflow of the 
MWPC anode readout. 

Another source of a systematic error resulting from the 
energy reconstruction as discussed in detail in Section 
VIE (and shown in Fig. 29) was the observed systematic 
~ 10-14 keV difference between the fitted endpoints and 
the Monte Carlo prediction. We investigated the sys- 
tematic uncertainty to the extracted asymmetry due to 
this systematic difference by extracting a "stretching fac- 
tor", / = T^Mc/To.fit, where T ,mc denotes the Monte 
Carlo predicted endpoint and To imcas the fitted endpoint, 
for each run. The data were then re-analyzed by apply- 
ing on an event-by-event basis this "stretching factor" 
to the reconstructed energy E recon , thus effectively forc- 
ing the fitted endpoints to match the Monte Carlo pre- 
dicted endpoints. The asymmetries extracted from the 
"strcched" data differed by < 0.07% from the (original) 
"unstretched" data which, again, is much less than the 
0.47% error associated with the energy calibration. 

Finally, to account for a slight mismatch (~ 2 keV) 
between the Monte Carlo and measured energy spectra 
(this is visible in the final E Tecon spectrum later in Fig. 
40) we fitted the Monte Carlo visible energy spectra to 
the measured visible energy spectra, and then extracted 
values for the asymmetry assuming these modified values 
for the visible energy. The bias to the asymmetry was 
0.13% averaged over all four Geometries which, again, is 
much less than the 0.47% the energy calibration uncer- 
tainty. 



D. Geometry-Independent: Fiducial Cut and 
Coordinate Systems 

As discussed in Section VF 2, we required backscatter- 
ing events to satisfy a default vertex cut of \xe~ xw\ < 25 
mm. We studied the impact of this cut on the asymme- 
try by varying this cut from 10 mm to 40 mm; the effect 
on the asymmetry was < 0.1%, indicating a negligible 
systematic effect. 

As was also discussed there, our fiducial cut required 
the position (radius) of the event on the primary trigger- 
ing scintillator side to satisfy rigger < 45 mm. To exam- 
ine whether there was any position bias, we extracted the 
asymmetry in successive annuli via cuts on r 2 in six differ- 
ent annular bins, ranging from [0,400] mm 2 to [2025,2500] 
mm 2 . The asymmetries in all of these annuli were in sta- 
tistical agreement, with no statistical evidence for any 



systematic difference with position. 

Recall also in Section VF2 we noted the possibility 
for the definition of four different coordinate systems. 
To determine whether there was any bias resulting from 
the choice of the coordinate system (for example, a con- 
sideration could be whether there were any systematic 
variations in the backscattering fractions in the vicin- 
ity of the fiducial cut), we studied the variation of the 
asymmetry with the choice of coordinate system, and for 
fiducial volume radius cuts of 45 mm and 50 mm. The 
RMS spread in the asymmetries for the different coor- 
dinate system choices was 0.24% for the 45 mm radius 
cut and 0.21% for the 50 mm radius cut. Although the 
RMS spread for the 50 mm radius cut was actually some- 
what smaller (suggesting that employing a larger fiducial 
volume would have introduced no bias to the asymme- 
try), we nevertheless chose the 45 mm radius cut as our 
(conservative) definition of the fiducial volume, and thus 
assigned a 0.24% systematic uncertainty to the definition 
of the fiducial volume. 



E. Geometry-Independent: Gain Fluctuations 

As noted in Section IIIF, the PMT gains were 
monitored on a run-to-run basis using the minimium- 
ionizing peak from cosmic-ray muon events. Neverthe- 
less, any residual uncompensated run-to-run gain fluc- 
tuations could bias the asymmetry on a run-to-run ba- 
sis; however, any such short-term run-to-run fluctuations 
will average away according to the usual 1 / \/N Tun statis- 
tics assuming the long-term gain corrections are accurate. 
We estimated the level of any such run-to-run residual 
gain errors by extracting the level of fluctuations in the 
run-to-run fitted values for the /3-decay spectrum end- 
point. These were typically of order ~ ±1.2% in each 
detector, with the gain fluctuations in the two detectors 
only slightly correlated relative to each other. [Correlated 
gain fluctuations arc significantly more problematic than 
are anti-correlated gain fluctuations.] Conservatively as- 
suming the worst-case sensitivity for gain fluctuations in 
one of the Geometries to be representative of the entire 
data set, we quote a systematic uncertainty of 0.20% for 
uncompensated gain fluctuations. 



F. Geometry-Independent: Live Time 

The detector rates (and, hence, asymmetries) were ul- 
timately calculated from the number of events passing 
the analysis cuts normalized to the detectors' respective 
live times; the concept of the detector live time was dis- 
cussed in detail previously in Section VD. As discussed 
there, we defined a run's live time to be the fraction of 
that run surviving all of the global data quality cuts. 
However, as we noted there, it was necessary to apply 
a correction for the Geometry B live times due to the 
large fraction (up to ~ 30%) of events suffering from an 
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event-by-event TDC event counter problem. The correc- 
tion factors for each run were determined using events 
identified as gamma rays, which were statistically inde- 
pendent of the neutron /3-decay events and also provided 
higher statistics (event rates up to 100 s _1 in each de- 
tector) than the neutron /3-decay events themselves for 
the calculation of the correction factors. The resulting 
correction factors, defined to be the ratio of the num- 
ber of gamma ray events surviving the event-by-event 
TDC event counter cut to the total number of gamma 
ray events, were then computed on a run- by-run basis 
for each detector. 

Only the Geometry B live times were corrected accord- 
ing to this procedure. Nevertheless, to assess the system- 
atic error associated with our definition of and calculation 
of the live time, we extracted values for the asymmetries 
for all four Geometries with and without application of 
these live time correction factors (the correction factors 
for Geometries A, C, and D were small, with the values 
for the asymmetries differing by < 0.1% under the two 
scenarios). Averaged over all four Geometries, the dif- 
ference between the asymmetries extracted under these 
two different scenarios was 0.24%, which is the value for 
the systematic uncertainty associated with this effect we 
quote in Tabic X. 



-o o% systematic uncertainty to this effect. We also note 
that our Monte Carlo calculations found that the time 
for a trapped electron to scatter from residual gas for a 
vacuum pressure of 10~ 5 Torr is ~ 4 ms, with a small 
distortion to their energy distribution of ~ — 8 keV. 

Note that our Monte Carlo results for the fraction of 
events trapped by the field dip and their average value 
of (/3 cos 9) are consistent with the following simple esti- 
mates. As discussed earlier in Section III El, electrons 
emitted with some momentum po = {p\ o+P 2 o) 1 ^ 2 ' w ^h 
P±,o (p\\.o) the initial transverse (longitudinal) momen- 
tum component, in some local field Bo will be reflected 
from higher field regions B if B > B clit = {p\jp\ )B 
(thus, only the pitch angle 9 of the emitted electron is 
relevant, not the magnitude of the momentum). Taking 
B a = 0.9925 T and B = 0.9955 T for the measured 2009 
field profile (here, B is taken to be the average of the lo- 
cal maxima at z = —100 cm and +50 cm) shown in Fig. 
4, one finds electrons with pitch angles 9 > 9 a n = 86.9° 
will be trapped in the field dip region. 

Approximating the initial angular distribution of emit- 
ted electrons as isotropic (reasonable, given that the /3- 
asymmetry is an O(10%) effect), one finds that the frac- 
tion of electrons emitted in the local field dip region 
B a = 0.9925 T which will be trapped is 



G. Geometry-Independent: Magnetic Field 
Nonuniformity 

Our Monte Carlo calculations of the corrections for 
backscattcring and the (/3 cos 9) acceptance discussed in 
Section VII D assumed a uniform magnetic field in the de- 
cay trap region. We studied the impact of the actual mea- 
sured nonuniformity in the spectrometer magnetic field 
shown previously in Fig. 4 in Monte Carlo. Qualitatively, 
the impact of the ~ 30 Gauss "field dip" in the central 
decay-trap region is such that electrons from decays oc- 
curing in this "field dip" region are cither reflected (anal- 
ogous to backscattering) or are trapped (for large pitch 
angles). We studied these effects in our GEANT4 Monte 
Carlo simulation program by implementing the magnetic 
field profile shown in Fig. 4 directly in the simulation. 
Neutron /3-decay events were then generated uniformly 
along the length of the decay trap. 

In the Monte Carlo, ~ 0.3% of the events incident ini- 
tially on one of the two detectors were reflected from the 
field dip, with an average (/3cos#) of ~ 0.02. Because 
this small fraction of events carries little (/3 cos 9) "an- 
alyzing power" , the resulting bias to the asymmetry is 
negligible. The fraction of electrons trapped by the field 
dip was ~ 2.6%, again with an average ((3 cos 9) of ~ 0.02. 
The remaining 97.1% of the events were not impacted by 
the field dip. Assuming that the electrons trapped by the 
field dip eventually scatter from residual gas molecules, 
the impact is a dilution to the asymmetry. The calcu- 
lated dilution to the asymmetry was —0.2%. In lieu of 
applying a correction to the asymmetry, we assigned a 



1 nit/ 2 ,2tt 

/trap = -r- 2 x / s'm 9 d9 / d(j) = cos6> crit 

4?F Jfcrit JO 

= 0.054. (34) 

Then, assuming a uniform distribution of events along 
the 300-cm long decay trap, the fraction of events emitted 
in the ~ 150-cm long field dip region is ~ 0.5, implying 
the total fraction of events generated over the length of 
the decay trap which will be trapped in the field dip 
region is ~ 0.027, which is consistent with the Monte 
Carlo result of 2.6%. For a nominal value of (3 ~ 0.75, 
{/3cos9) « /3cos[(6» crit + 7r/2)/2] = 0.02 for these trapped 
events, again, consistent with the Monte Carlo result. 



H. Geometry-Independent: Muon Veto Efficiency 

We estimated the effect of a possible systematic uncer- 
tainty resulting from fluctuations in the muon-veto effi- 
ciency by extracting values for the asymmetries with and 
without application of the muon-veto detector cuts. Av- 
eraged over Geometries, the variations in the asymmetry 
were at the 0.3% level. 

We note that the assignment of this 0.3% uncertainty 
is quite conservative. A linear drift in the muon veto cut 
efficiency would be equivalent to a linear drift in the back- 
grounds, and as discussed in Section VIC, linear back- 
ground drifts cancel under the octet-based super-ratio 
asymmetry structure. 
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I. Geometry-Independent: Neutron-Generated 
Backgrounds 

As already dicussed in detail, ambient backgrounds 
were measured and subtracted on a run-by-run basis. 
However, a possible source of irreducible backgrounds 
was neutron capture on materials near the electron de- 
tectors, generating prompt gamma rays with energies up 
to 7.9 MeV, 7.1 MeV, 6.8 MeV, 4.9 MeV, or 8.2 MeV for 
capture on 63 Cu, 65 Cu, 9 Be, 12 C, or 13 C, respectively, 
the elements of which the decay trap and end-cap foils 
were primarily composed. Such backgrounds cannot, of 
course, be subtracted. 

This background was expected to be significantly sup- 
pressed in the UCNA experiment as compared to previ- 
ous cold neutron beam experiments because, as discussed 
earlier in Section III A, the fraction of neutrons present 
in the apparatus which contribute to the decay rate is or- 
ders of magnitude larger in the UCNA Experiment than 
in previous cold neutron beam experiments, and also be- 
cause of the small probability for capture and upscatter 
by UCN stored in the decay trap. 

We carried out three different approaches to our as- 
sessment of the contamination level from any such back- 
grounds. The idea of our first approach is as follows. If a 
gamma ray emitted from a neutron capture subsequently 
interacted with the scintillator, the MWPC should not 
have recorded any energy deposition if the gamma ray 
forward Compton scattered in the scintillator. Further, 
as calculated in simulations, there is a factor of 10-20 
suppression in the fraction of gamma ray events incident 
on the electron detectors triggering both the scintilla- 
tor and MWPC as compared to those triggering only the 
scintillator. Therefore, any such neutron-generated back- 
grounds should appear as non-zero residuals in a compar- 
ison of background-subtracted scintillator spectra [i.e., 
(/3-decay run — background run) spectra] formed with 
and without application of a MWPC-scintillator coinci- 
dence cut. In particular, an excess would be expected in 
the background-subtracted spectrum constructed with- 
out the requirement of a MWPC-scintillator coincidence 
cut as compared to the background-subtracted spectrum 
obtained with the requirement of a MWPC-scintillator 
coincidence cut. 

Now consider the following model. Under application 
of a MWPC coincidence cut, let 

-Rcut = £ (Scut ~ B cut ) , (35) 

where i? C ut denotes the resulting background-subtracted 
scintillator event rate, S C ut and B cnt denote, respectively, 
the underlying /3-decay + background and background 
event rates, respectively, and e denotes the MWPC cut 
efficiency. We then write a similar expression for the 
background-subtracted scintillator event rates obtained 
without application of a MWPC cut as 

-Rno cut = (Scut — -Bcut) + (5*7 — By) , (36) 



where now S 7 and By denote the signal and background 
rates during the /3-decay and background runs from 
gamma ray events which would otherwise fail the MWPC 
cut. Note that in the absence of any neutron- generated 
gamma rays, the statistical averages of S 7 and By should 
be identical. 

The difference between R no cu t and i? C ut is then 

AR = R no cut — i?cut 

= (S 7 - B 7 ) + (Scut - Bcut) (1 - e) 

= Ty + (1 - e)r„, (37) 

where T 7 and T n denote, for gamma ray and neutron 
/3-decay events, respectively, the difference between the 
background-subtracted scintillator rates with and with- 
out application of a MWPC coincidence cut. Thus, this 
model then requires an estimate for the MWPC cut effi- 
ciency. 

We extracted a value for our MWPC cut efficiency by 
examining the MWPC anode spectrum for 113 Sn source 
calibration data. After placing a FWHM cut on the scin- 
tillator visible energy spectrum, we then fitted the re- 
sulting MWPC spectrum (such as shown, for example, 
in Fig. 10) to a Landau distribution, and its pedestal to 
a Gaussian. We then calculated the fraction of the Lan- 
dau distribution falling below the cut line, which should 
provide an estimate of the MWPC cut efficiency. These 
values were found to be 99.93(3)% and 99.95(3)% for the 
East and West detector, respectively. We do note that 
this provides for an estimate of the MWPC efficiency only 
at the 113 Sn energy. 

The residual AR (no MWPC cut - MWPC cut) rates 
integrated over the analysis energy window (and after all 
analysis cuts) ranged from 1.0 x 10 -3 s _1 to 6 x 10 -3 s _1 
for the two detectors and two spin states for all four Ge- 
ometries, representing ~ 10~ 3 of the /3-decay rates. Af- 
ter accounting for the factor of 10-20 supression for the 
fraction of events which would trigger both the scintilla- 
tor and MWPC, the estimated contamination fractions 
for the actual /3-decay analysis employing the MWPC- 
scintillator coincidence cut are then on the order of 10~ 4 . 
Propagation of the measured contamination fractions (for 
each detector and spin state) through the super ratio 
then led to a systematic bias to the asymmetry of order 
- 0.02%. 

In our second approach, we extrapolated the residual 
background (i.e., after background subtraction) above 
the /3-decay endpoint into the signal region. Above 
the endpoint, the residual background-subtracted rates 
in the 25 keV E iecon bins were typically of order 10~ 4 
s _1 or less. Under the assumption that the neutron- 
generated background is independent of energy (e.g., as 
was employed in the analysis of [46]), we then extrap- 
olated these above-the-endpoint rates into the analysis 
energy window. The resulting contamination fractions, 
~ 10~ 4 -10 -3 , were similar to the analysis in our first 
approach comparing background-subtracted scintillator 
spectra obtained with and without a MWPC cut. These 
contamination fractions were again propagated through 
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FIG. 38. Schematic diagram of the "beta-blocker" measure- 
ment of neutron-generated backgrounds. The dashed lines are 
meant to indicate the solid angle for production of Compton- 
scattered electrons in the acrylic at the two different positions. 



J. Geometry-Independent: Polarization 

The UCN polarization systematic was discussed earlier 
in Section IV C. Because the measured depolarization 
was consistent with zero at the la level (i.e., P > 0.9948), 
we did not apply a correction for the polarization, and 
instead quote a one-sided systematic uncertainty in Aq 



of 



+0.52% 
-0.00% 



resulting from the constraint P > 0.9948. 



K. Geometry-Independent: Rate-Dependent Gain 
Shifts 



our super-ratio asymmetry analysis, and the systematic 
bias was again found to be of order ~ 0.02%. 

Finally, in our third approach, we carried out a "beta- 
blocker" measurement in which a 6.35-mm thick piece 
of acrylic was placed between the decay trap and the 
MWPC in the field-expansion region of the spectrom- 
eter, as indicated schematically in Fig. 38. The idea 
for this measurement was two-fold: (1) the acrylic was 
sufficiently thick to stop the endpoint /3-decay elec- 
trons, thereby "blocking" the /3-dccay signal of inter- 
est; and (2) the acrylic then served as a "source" of 
Compton-scattered electrons, resulting from interactions 
of neutron-generated gamma rays with the acrylic. 

Measurements were conducted with this piece of 
acrylic at two different positions, A and B, as shown in 
Fig. 38, in front of one of the detectors. The motiva- 
tion for doing so was that a comparison of the results 
from positions A and B should, in principle, permit a 
decomposition of the measured detector signal into con- 
tributions from Compton-scattered electrons originating 
in this acrylic piece (the solid angle for which was clearly 
smaller in position B as compared to position A), and 
direct neutron-generated gamma-ray interactions in the 
plastic scintillator detector (which should not have varied 
with the position of the acrylic piece) . 

Note that this measurement is subject to some model 
dependence, including an assumption for the source posi- 
tions along the decay trap of the neutron-generated back- 
grounds (which determines the ratio of the solid angles 
for positions A and B). Our resulting estimates for the 
contamination fraction, as extracted from our measure- 
ments of the residual (background-subtracted) rates with 
the acrylic piece located at both positions A and B (under 
the assumption that the ratio of the A and B solid angles 
for production of Compton-scattered electrons was 20:1), 
were of order 10 -3 . Then accounting for the factor of 10- 
20 suppression for gamma ray events triggering both the 
scintillator and MWPC (in the actual geometry) leads to 
an estimate for the contamination fraction on the order 
of 10~ 4 , consistent with the other two approaches. 



A potential systematic effect would arise from any 
spin-state-correlated systematic gain shifts, such as from 
rate-dependent gain shifts. However, any such gain shifts 
shared by both detectors cancel to first order in the super 
ratio asymmetry and are thus expected to be small. As 
noted earlier in Section IIIB, during the operation of the 
experiment, the total DAQ recorded data rates (i.e., the 
"online" /3-dccay rates integrated over all energies with 
few cuts) during /3-decay runs for the non-flipped spin 
state were typically ~ 10 s _1 greater than those recorded 
during measurements of the flipped spin state. By taking 
data with calibration sources with different activities, we 
were able to bound any such rate-dependent gain shifts 
to then be < 0.02%/(~ 10 s _1 ), which corresponds to a 
systematic uncertainty in the asymmetry of < 0.08%. 



L. Geometry-Dependent: Backscattering A2 and 
Angle Effect A3 Corrections 

Our working assumption was that application of the 
/bulk and /thin scale factors to the GEANT4 backscatter- 
ing distributions calibrated our Monte Carlo calculations 
of our A2 backscattering corrections to the asymme- 
try. To estimate the uncertainty in these now-calibrated 
corrections, we compared the GEANT4 results for the 
A2 backscattering correction with the PENELOPE results. 
[Note that the PENELOPE calculations required /bulk and 
/thin scale factors of 0.9-1.1 and 1.3, respectively, some- 
what smaller than those required by GEANT4.] These 
agreed to better than 22% for Geometry A, and better 
than 6% for Geometries B, C and D. We also note that 
the RMS of the A2 corrected asymmetries for Analysis 
Choices 1-5 was 0.10%, 0.13%, 0.30%, and 0.27% for Ge- 
ometries A, B, C, and D, respectively, providing a pow- 
erful check of the robustness of the calculation of the A2 
correction (indeed, consistent with the robustness of the 
agreement between the measured and simulated asym- 
metries for the various Analysis Choices demonstrated 
previously in Fig. 31). We have taken a conservative 
approach to our estimate of the systematic uncertainty 
in our backscattering corrections, and quote a 30% rela- 
tive uncertainty in the A2 backscattering correction (and, 
thus, in the asymmetry) for all of the Geometries. 
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FIG. 39. Results from GEANT4 Monte Carlo calculations of the relative size of the combined A2 backscattering and A3 angle effect 
fractional systematic correction (thick solid lines) and the total fractional systematic uncertainty assigned to these corrections 
(gray bands). 



To estimate the systematic uncertainty in the A3 angle 
effect correction, we varied the thickness of the decay trap 
end-cap foil thicknesses in the Monte Carlo. Assuming 
0.5 11m. to be a reasonable uncertainty in the foil thick- 
ness, the relative uncertainty in A3 was found to be no 
larger than 25%. Further, we note that the GEANT4 and 
PENELOPE results for A 3 agreed to better than ~ 25% for 
all of the Geometries. Therefore, we again quote a con- 
servative 25% relative uncertainty in the A3 angle effects 
correction for all of the Geometries. 



Figure 39 shows the resulting energy-dependent error 
bands for the combined A2 and A3 corrections for each 
of the Geometries. The analysis energy window of 275- 
625 keV referenced earlier in Section VII F was chosen 
such that the total error in Aq resulting from integration 
of these error bands over some energy window combined 
with the statistical error within that energy window was 
a global minimum. 



M. Geometry-Dependent: MWPC Efficiency 

As discussed earlier in Section VM, the "standard" 
MWPC cut for the separation of gamma rays and charged 
particles was a cut on a fixed PADC channel number. 
However, as shown there, the MWPC exhibited a strong 
position-dependent response. Therefore, for some partic- 
ular energy deposition in the MWPC, employing such a 
standard cut resulted in a position-dependent efficiency 
for the identification of an event as either a gamma ray 
or charged particle event. 

We investigated the impact of this position-dependent 
efficiency on the asymmetry by comparing results for the 
asymmetry extracted from analyses employing the stan- 
dard PADC channel number cut with those obtained with 
fixed MWPC energy cuts ranging from 0.1-0.7 keV. The 
error bounds from this analysis appear as the MWPC 
efficiency systematic uncertainties in Table X. Note that 
this effect is largest for Geometry D, for which the dif- 
ferences between the two MWPCs' efficiency curves was 
greatest, as shown previously in Fig. 22. 
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FIG. 40. (Color online) Upper plot: Final results for the 
measured background -E rcc0 n energy spectrum (open circles) 
and the background-subtracted neutron /3-decay Jf? rC con en- 
ergy spectrum (filled circles) summed over both detectors, 
averaged over the two neutron spin states, and then averaged 
over all four Geometries. The Monte Carlo prediction for the 
Brecon spectrum is shown as the solid line. Lower plot: Final 
results for Aq extracted bin-by-bin and then averaged over 
all four Geometries. The solid red line indicates the analysis 
energy window of 275-625 keV. The quoted error is statistical. 



IX. SUMMARY OF FINAL RESULTS 

A. Energy Spectra and Binned Ao Results 

Our final results are shown in Fig. 40, where we com- 
pare the background-subtracted /3-decay -Brecon spec- 
trum, summed over both detectors, averaged over the 
two neutron spin states, and averaged over all four of 
the Geometries, with the Monte Carlo predicted -Brecon 
spectrum. There, we also show our Geometry-averaged 
energy-binned values for the /3-asymmetry Aq, and the 
final statistical result of Aq = -0.11966 ± 0.00089. Note 
that the central value for Aq was insensitive to the choice 
of analysis energy window, with the variation less than 
15% of the statistical uncertainty for other windows be- 



tween 150 and 750 keV. 



B. Final Combined Result 

Each of the four Geometries yielded a data set with a 
statistical error and a systematic error. The statistical 
errors for each of these Geometries were, of course, in- 
dependent. However, three of our dominant systematic 
uncertainties (energy reconstruction, backscattering, and 
angle effects) were all correlated. For example, a mistake 
in the decay trap end-cap window thickness would have 
biased the angle effects correction for all of the data sets. 
Therefore, we assume all three of these systematic un- 
certainties are 100% correlated across the four different 
Geometries. 

Under this assumption, we then combined the results 
from the four Geometries according to the following pro- 
cedure which incorporates correlations properly in the 
construction of a global % 2 [74]. Each individual mea- 
surement i gives a constraint of 



(38) 



k=i 



where A and &i denote, respectively, the central value 
and statistical uncertainty for the i th measurement (i.e., 
Geometry), and the cr(A^) denotes the correlated uncer- 
tainty due to the three systematic effects just discussed. 
We then constructed a x 2 as 



££(4 -A )(V-%(Al-A ), 



(39) 



where i and j are the indices of the measurements, 
and V is the covariance matrix with Vij — erf Si j + 
J2l =1 cr(^i)cr(Al). As was shown in [74], the \ 2 con- 
structed this way satisfies a standard x 2 distribution, in 
which the value of Aq follows from minimization of this 
X 2 , and the one-sigma uncertainty is obtained from the 
usual condition \ 2 — Xmin = 1- 

Using standard minimization techniques, 
combined result for Aq we obtained is 



the final 



A 



-0.11966 ± 0.00089 t^oouo) 
[ X 2 /v = 2.4/3 (Prob = 0.49)] 



(40) 



where the first (second) error represents the statistical 
(systematic) error. From this, we extract the following 
value for A = QaI 9v under the Standard Model, 



A = — 

9v 



4.27590 ±0.00239 



+0.00331 
-0.00377- 



(41) 



X. SUMMARY AND CONCLUSIONS 

In this article we have presented a comprehensive 
and detailed description of the first precision result [39] 
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FIG. 41. (Color online) Ideogram of values for gA/gv ex- 
tracted from measurements of the neutron /3-asymmetry pa- 
rameter Aq (this work [39], filled red square; all other mea- 
surements of Ao [40, 42, 44, 47], open red squares), a simul- 
taneous measurement of Ao and Bo (black triangle, [75]), the 
current average value of gA/gv — — 1.2726±0.0024 from these 
measurements (gray band; note that the current 2012 Particle 
Data Group average value for gA/gv has not yet incorporated 
the result of [47]), measurements of the neutron lifetime r n 
(blue circles, [76-82]) assuming the superallowed + — > + 
value for V u d [35], and the current 2012 Particle Data Group 
average value for the neutron lifetime of 880.1 ±1.1 s [34] 
(between dashed lines). 

from the UCNA Experiment, an experiment designed to 
perform the first-ever measurement of the neutron /3- 
asymmctry parameter Aq with polarized ultracold neu- 
trons. As demonstrated here, the use of UCN in a 
neutron /3-asymmetry experiment controls key neutron- 
related systematic corrections and uncertainties, includ- 
ing the neutron polarization and neutron-generated back- 
grounds. Our result for the neutron polarization was 
shown to only be statistics limited, and our neutron- 
generated backgrounds were negligible to the level of 
< 0.1% precision. All of our results reported here are con- 
sistent with our previously published proof-of-principle 
results obtained during data-taking runs in 2007 [38]. 

To evaluate the immediate impact of this work, we first 
note, as shown earlier in Fig. 1, that our measurement 
agrees well with the most recent (and most precise) pub- 
lished result for the /3-asymmetry Aq from the PERKEO 
II experiment [45-47], but is in poorer agreement with 
the three other results [40-44] employed by the Particle 
Data Group in their averaging procedure. Because of the 
differences between the PERKEO II and UCNA experi- 
mental techniques, we believe this is a significant result. 
Second, the value for gA/gv that we extract from our 
measurement is 1.27590 to!ou445- We compare our value 
for gA/gv with results from a global fit under the Stan- 
dard Model for gA/gv values extracted from results for 



the /3-asymmetry Aq [39, 40, 42, 44, 47], a simultaneous 
measurement of Aq and the neutrino asymmetry Bq [75] , 
and from individual measurements of the neutron life- 
time [76-82] and the current Particle Data Group average 
value for the lifetime (880.1 ± 1.1 s [34]; recently updated 
for the corrected result of [82] which supersedes the orig- 
inal result of [83]), where the lifetime results assume the 
superallowed 0+ -> 0+ value for V ud of 0.97425 ± 0.00022 
[35] in Eq. (3). The results of our global fit are dis- 
played as an ideogram in Fig. 41, where it can be seen 
that under the Standard Model the PERKEO II [47] and 
UCNA [39] /3-asymmetry experiments are in agreement 
with the three most recent (or updated) results for the 
neutron lifetime reported from experiments with stored 
UCN [76, 77, 82], but in poorer agreement with the other 
results for the lifetime employed by the Particle Data 
Group in their averaging procedure obtained in exper- 
iments utilizing stored UCN [78, 79] and cold neutron 
in-beam [80, 81] techniques. Indeed, from our result for 
Aq alone of —0.11966 1q!ooi66' we extract, according to 
Eq. (3), a value for the neutron lifetime of 

r n = 879.0 tfl s, (UCNA A and 0+ -> + V ud ) (42) 

in agreement with the measured values from the three 
most recent results reported from experiments using 
stored UCN [76, 77, 82]. 

Thus, we conclude that our /3-asymmetry measurement 
already provides significant impact to a self-consistent 
evaluation of the landscape of neutron /3-decay observ- 
ables and the superallowed + — > + V u d data set. The 
fact that the most recent values for the /3-asymmetry and 
the neutron lifetime currently exhibit statistically signifi- 
cant deviations from their respective world averages pre- 
pared by the Particle Data Group, but are seen to be in 
agreement with each other (as shown in Fig. 41), moti- 
vates further refinement of the UCNA technique, with its 
novel approach to key neutron-related systematic errors, 
in order to conduct a more precise evaluation of neutron 
/3-decay observablcs. Indeed, recently demonstrated im- 
provements to our UCN source and refinements to the 
energy calibration and gain monitoring systems will per- 
mit the future collection of a data set with significantly 
improved statistical and systematic uncertainties. 
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